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ABSTRACT 

An adaptive numeric predictor-corrector guidance algorithm is developed for atmos- 
pheric entry vehicles which utilize lift to achieve maximum footprint capability. Applicability 
of the guidance design to vehicles with a wjde range of performance capabilities is desired 
so as to reduce the need for algorithm redesign with each new vehicle. Adaptability is 
desired to minimize mission-specific analysis and planning. The guidance algorithm moti- 
vation and design are presented. 

Performance is assessed for application of the algorithm to the NASA Entry Research 
Vehicle (ERV). The dispersions the guidance must be designed to handle are presented. 
The achievable operational footprint for expected worst-case dispersions is presented. The 
algorithm performs excellently for the expected dispersions and captures most of the achiev- 
able footprint. 
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= acceleration vector 
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= total inertial acceleration vector 
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— velocity vector term in the integration algorithm 
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= longitude 

= coefficient of viscosity for air 
= Earth gravitational constant 
= Earth's rotation vector 

= natural frequency of heat rate control response 
= bank angle 
= universal gas constant 
= atmospheric density 

— standard deviation 

= time constant of filter 

— damping ratio of heat rate control response 


SUBSCRIPTS 


aero = aerodynamic 

c — geocentric 

cmd = command 

d — desired 

des = desired 

drag = drag term 


18 



e 

El 

f 

g 

imu 

inplane 

lat 

lift 

lim 

LID 

max 

min 

nom 

per pen 

pole 

R 

si 

std 

t 


= error 

= entry interface 
= final 
= geodetic 

= inertial measurement unit 

— projection of target unit vector into plane formed by the 
position vector and the relative velocity vector 

= direction perpendicular to the plane formed by the position vector 
and the relative velocity vector 
= lift term 

= limiting value or boundary 
= lift-to-drag ratio 
= maximum 
= minimum 
= nominal 

= direction perpendicular to the plane formed by the position vector 
and the relative velocity vector 
= direction of the north pole 
= direction of the position vector 
= sea level 
= standard value 

— target aim point 


P 


— atmospheric density 


19 



SUPERSCRIPTS 


EF 

imu 

imu past 

A 


coordinatized in Earth-fixed coordinates 

value measured by inertial measurement unit on current cycle 
value measured by inertial measurement unit on past cycle 
estimated or measured 
value from previous guidance cycle 


20 



1.0 INTRODUCTION 


Routine access to space and the maintenance of a Space Station will increasingly 
require greater flexibility in mission planning and the requirement for lower system mainte- 
nance costs. The launch and recovery phases of space flight have historically been the 
most demanding phases of space flight and therefore require the most development effort 
and investment. Mission flexibility requires more frequent launch and deorbit opportunities. 
For the case of re-entry vehicles, deorbit opportunities are defined by the ranging capability 
of the vehicle. A high L/D vehicle increases the available deorbit opportunities increasing 
mission flexibility. High L/D vehicles also are of interest for over-flight missions for the pur- 
pose of reconnaissance. 

Entry guidance algorithms developed to date have been highly vehicle-specific and 
required great development and maintenance efforts over the life of the vehicle. These algo- 
rithms were not applicable to other vehicles without extensive modification. 

This study seeks to design an adaptive entry guidance algorithm that maximizes the 
usable footprint by making full use of the available vehicle capability. This algorithm should 
also be easy to maintain throughout the vehicle definition phase and operational life. Mini- 
mizing the number of mission-dependent input parameters (l-loads) is desirable. The algo- 
rithm should also be easily transported to other vehicles to minimize development cost. 
Transportability is accomplished by minimizing vehicle-specific features of the algorithm. 
Explicit heat rate control should be provided to allow full use of the entry corridor up to the 
heat rate limits. 
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This study seeks to design such an algorithm. A candidate entry guidance algorithm is 
defined for the NASA Entry Research Vehicle (ERV), but is easily adapted to other vehicles 
with minimal modification. The proposed algorithm attains almost complete coverage of the 
achievable footprint, while employing a simple one-phase entry algorithm with explicit heat 
rate control. Vehicle-specific features and l-loads are minimized, reducing algorithm devel- 
opment and maintenance costs. 

The ERV [1] is a proposed high-performance entry vehicle designed as a test bed for 
future technology development in the areas of: 

1. Maneuvering entry/synergetic plane change 

2. Atmospheric uncertainties 

3. Advanced thermal protection systems 

4. Aerodynamic/aeroheating prediction 

5. Adaptive guidance and navigation 

6. Load-bearing thermostructures 

The ERV is designed for deployment from the Space Shuttle, after which the ERV enters the 
atmosphere for demonstration of the synergetic plane change, over-flight, and entry mis- 
sions. Figure 1 on page 78 shows a three-view drawing of the ERV and the surface areas of 
the aerodynamic control surfaces. Also seen is the size of the ERV in relation to the diam- 
eter of the Shuttle payload bay in which the ERV musrfit 
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2.0 MOTIVATION 


2.1 INTRODUCTION 


The goal of any entry guidance algorithm is to successfully guide the vehicle to the 
desired final state for the largest range of dispersions possible without violating any vehicle 
constraints while also maximizing the achievable footprint. It is also desirable to minimize 
the mission and vehicle-specific aspects of the guidance algorithm so as to minimize pre- 
mission analysis and planning. Transportability of the algorithm from one vehicle to another 
significantly reduces guidance algorithm development effort and cost. 

To maximize the footprint attainable, the guidance algorithm must follow the optimal 
path to any particular point in the footprint. The algorithms developed to date for such vehi- 
cles as the Apollo capsule [2] and the Space Shuttle [3] have attempted to do this by fitting 
the optimal trajectory with phases that follow important parameters (reference profiles) over 
some range of conditions. These guidance algorithms were required to be computationally 
efficient because of the limited on-board computer resources available. Analytic 
expressions for the reference profiles allowed for low execution time and tailoring of the tra- 
jectory for vehicle-specific constraints. For example, trajectories for these vehicles had to 
be shaped to reduce and control the maximum heat rate experienced below that allowed for 
the available thermal protection system materials. 
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The Space Shuttle entry guidance system employs three major modes with seven phas- 
es: 

1. Entry 

a. Pre-entry 

b. Temperature control 

c. Equilibrium glide 

d. Constant drag 

e. Transition 

2. Terminal Area Energy Management 

3. Approach and Landing 

Except for the pre-entry phase which is open-loop, each phase is described by an analytic 
expression relating the desired drag and altitude rate (the measured feedback terms used) 
to the desired profile. Because the algorithms are tailored for a particular vehicle and the 
reference profiles do not follow the optimal profile to all points in the footprint, guidance 
algorithms developed to date can not be easily adapted to other vehicles or provide full cov- 
erage of the. theoretically achievable footprint. 

The next generation of entry vehicles will not be so constrained due to advances in ther- 
mal protection system materials and computer technology. For example, flight computers 
are now capable of supercomputer speeds on the order of 40 million instructions per second 
utilizing parallel processing architecture [4]. A different approach to guidance that attempts 
to follow an optimal profile to maximize footprint capability is therefore possible. 

The proposed approach is a predictor-corrector algorithm that numerically predicts the 
final state for a particular control variable history and then corrects the control variable his- 
tory to satisfy the specified final state constraints. This approach, proposed previously for 


24 



various guidance problems, has most often been impractical because of the long trajectories 
that must be predicted and the slow computer speeds. 

Such an approach has been employed for the Space Shuttle Powered Explicit Guidance 
(PEG) [5] used for second stage ascent and orbit insertion burns where the trajectory is 
short enough to be predicted with the available computer resources. The Shuttle algorithm 
numerically predicts the gravitational effects during the powered flight phase with a 10 step 
integration of the 500 second trajectory. 

A predictor-corrector has also been proposed for Aerobraking Orbital Transfer Vehicles 
(AOTV) [6] which would utilize more advanced computers. This algorithm numerically inte- 
grates the equations of motion along a skimming trajectory through the upper atmosphere 
that is approximately 500 seconds long and requires about 100 integration steps. 

The trajectories flown by the ERV or any high L/D entry vehicle are typically from 30 to 
100 minutes long from entry interface (400K feet) to landing, so the computational demand 
for such an algorithm is very great early in the entry when the time to landing is long. How- 
ever, because the entry is long and the vehicle has excess ranging capability for all but a 
small region along the edge of the footprint, the accuracy of the early predictions need not 
be as high as for the later predictions. Hence, large time steps can be used early in the pre- 
dictor algorithm. Later, when the vehicle nears the landing site, the time remaining is short, 
and hence, the prediction is short. This allows the predictor-corrector to be executed more 
often near landing just like the current analytic algorithms. Throughout the entry, vehicles 
using an analytic guidance algorithm with reference profiles must closely follow the refer- 
ence profile if the assumed reference profile is to guide the vehicle to the correct final state. 
A predictor-corrector effectively recomputes a new reference profile each time it is executed, 
so the guidance execution rate can be much lower than that for analytic algorithms. 
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2.2 DISPERSIONS 


Before the guidance algorithm can be designed, the possible dispersions that may affect 
the trajectory must be considered. The Shuttle entry guidance system is required to reach 
the Terminal Area Energy Management (TAEM) interface with less than a 2.5 nautical mile 
position error from the target aim point. The dispersions of significance to an entry vehicle 
trajectory include: 

1. Vehicle characteristics 

a. Mass 

b. Aerodynamics 

c. Maneuver rates 

2. Environment characteristics 

a. Atmospheric density 

b. Atmospheric winds 

c. Atmospheric properties influencing aerodynamic flow regimes (temperature, 
mean free path, etc.) 

3. Initial entry state vector 

a. Velocity 

b. Flight path angle 

c. Heading 

4. Propagation errors in navigation state vector 

Of these potential dispersion sources, only the vehicle mass, aerodynamics, and the 
atmospheric density and winds will be significant. By the early 1990's, almost perfect navi- 
gation can be expected through use of the Global Positioning System (GPS). If the deorbit 
burn guidance and control systems are assumed to correctly guide to the navigated state 
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and there are no navigation errors, then the dispersions in the initial entry state vector are 
negligible. 

The vehicle mass should be known accurately, so for this study, a 3cr error of ±5% is 
assumed. Experience from the Space Shuttle program shows that the vehicle aerodynamics 
should be known to within +5% for the force coefficients on the first flight. Only the stability 
derivatives and control effectiveness were missed significantly [7]. Even though the force 
coefficients may be known to excellent accuracy, reduced control effectiveness can reduce 
the possible trim angle of attack range reducing the maximum L/D achievable. Therefore, 
for this study, a ±10% dispersion in the lift and drag coefficients is considered. It should be 
noted that the first few flights of a new vehicle are usually targeted to the middle of the foot- 
print to maximize margin and allow for accurate determination of the vehicle characteristics 
before the full ranging capability of the vehicle is used. After the first few flights, the aero- 
dynamic characteristics should be known to within a few percent, so only about a ±3% dis- 
persion must be considered. 

The atmospheric dispersions were obtained from two sources. Reference [8] specifies 
the atmospheric dispersions to which aerospace vehicles must be designed. The average of 
the steady state winds at four geographic locations is shown in Figure 2 on page 79. This 
model was incorporated into the simulator environment with a magnitude scale factor to 
simulate less than worst-case winds. The wind direction was selected for each run made 
with winds and held constant throughout the trajectory. Reference [8] specifies Reference 
[9] as the source for atmospheric density dispersions. However, the recent Shuttle flights 
have provided estimated density data of a quality never before available. Atmospheric den- 
sity profiles derived from Shuttle accelerometer measurements of the normal force acceler- 
ation and the estimated normal force coefficient and relative velocity vector are presented in 
Reference [10]. Figure 3 on page 80, taken from that report, shows the envelope of the 


27 



derived density profiles for the first 12 Shuttle flights. Of particular interest is the range of 
dispersions seen: -47% to +12%. Figures 4 on page 81 and 5 on page 82 show the density 
profiles for the STS-1 and STS-9 Shuttle flights. High frequency density shear components 
and constant density biases from the standard atmosphere are seen. For this study, con- 
stant density biases of ±30% and the Shuttle derived density profiles from Reference [10] 
were used. 


2,3 REFERENCE TRAJECTORIES 


The size of the footprint for a particular vehicle is determined by the range in vehicle 
L/D and the constraints placed on the trajectory such as heat rate limits. The edges of the 
footprint correspond to the use of maximum or minimum L/D. Maximum downrange or 
crossrange, for example, requires maximum L/D, while minimum downrange requires mini- 
mum L/D. 

The determination of the optimal angle of attack and bank angle control histories for 
maximum crossrange and downrange has been the topic of many papers [11] [12] [13] . 
Wagner [12] used several optimization techniques to evaluate the maximum crossrange 
achievable for a multiphase bank angle history flown at maximum L/D. The multiphase bank 
profiles considered are shown in Figure 6 on page 83. It is seen that as the number of phas- 
es increases, the multiphase profile approaches the optimal continuous profile also shown in 
this figure. It was determined that a three-phase bank angle profile as illustrated in Figure 6 
achieved almost the same crossrange as a continuous bank profile. This is shown in 
Figure 7 on page 84 reproduced here from that paper. Further, as the number of phases 
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increases, the optimum bank angle profile approaches a continuous profile that is almost lin- 
ear with velocity as shown in Figure 8 on page 85. It was also Lhown that flying at the maxi- 
mum L/D maximizes the crossrange attained. 

This result is confirmed in Reference [13] which utilized a nonlinear programming tech- 
nique to optimize the Space Shuttle trajectory for the maximum downrange and maximum 
crossrange cases. The maximum downrange trajectory requires flying at zero bank angle 
and at the angle of attack corresponding to maximum L/D as shown in Figure 9 on page 86. 
The control histories for the maximum crossrange case are shown in Figures 10 on page 87 
and 11 on page 88. Again, the optimal control history is the angle of attack corresponding to 
maximum L/D and an almost linear bank angle profile with velocity. 

Optimized trajectories for the ERV were reported in Reference [14], These trajectories 
were determined using the Program to Optimize Simulated Trajectories (POST) [15] and 
imposed the following constraints on the trajectories: 

1. Maximum heat rate of 125 BTU/sq ft/sec 

2. Maximum heat load of 150K BTU/sq ft 

The achievable footprint with these constraints, reported in Reference [14], is shown here in 
Figure 12 on page 89 . Subsequently, the heat load limit was increased to 175K BTU/sq ft 
resulting in the larger footprint shown in Figure 12. As will be seen, these footprints omit a 
large area in the minimum downrange region that is achievable within the heating con- 
straints. Also shown is the footprint of the Space Shuttle which has a maximum hypersonic 
L/D of 1.2 as compared with 1.8 for the ERV. 

Figure 13 on page 90 shows the altitude history for the maximum downrange, maximum 
crossrange, and minimum downrange cases. Figures 14 on page 91 and 15 on page 92 
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show the bank angle and angle of attack histories for these trajectories. Figures 16 on page 
93 and 17 on page 94 show the heat rate and heat load histories for these cases. 

Figure 15 shows that the constant angle of attack corresponding to maximum L/D is 
flown for the edge of the footprint except for the minimum downrange case. For the mini- 
mum downrange case, the angle of attack corresponding to the minimum L/D on the back 
side of the L/D curve (high drag coefficient) is flown early, followed by a ramp in angle of 
attack starting at 1500 seconds after entry interface. This ramp corresponds to the vehicle 
actually turning around and flying slightly back uprange, so maximum L/D is desired later to 
maximize the distance flown uprange.. The angle of attack for the maximum downrange 
case is slightly greater than that for maximum L/D because this trajectory exceeds the heat 
load limit if flown at maximum L/D. The maximum downrange region of the footprint is 
therefore limited by the heat load limit set for the ERV. If the limit were relaxed, flight at 
maximum L/D would allow a longer downrange trajectory. 

Figure 14 shows that the bank angle profile for maximum crossrange is approximately 
linear with time which is almost linear with velocity, which suggests that a linear bank angle 
profile with velocity is sufficient. The maximum downrange case has a constant bank angle 
of zero which is again linear with velocity. The minimum downrange case does not have a 
linear bank profile. As was mentioned previously, for this case, the vehicle turns around and 
flys back uprange. 

The results of these studies suggest that use of a constant angle of attack profile and a 
linear bank with velocity profile will capture a large portion of the achievable footprint. As 
will be seen in the results, these profiles suffice to capture most of the footprint reported in 
Reference [14] and additionally reach a large area in the minimum downrange region out- 
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side the reported footprint. Only a small area of the reported footprint in the minimum 
downrange region is unachievable. 

Also of interest are the peaks in heat rate seen in Figure 16. Because the peaks in heat 
rate are very short, explicit control of the heat rate should be possible in the maximum heat 
rate regions without significantly impacting the guidance. 


2.4 GUIDANCE APPROACH 

The guidance design will attempt to maximize the size of the footprint while flying a con- 
stant angle of attack profile and a linear bank angle with velocity profile. The predictor algo- 
rithm integrates the equations of motion forward in time using the assumed control profile 
and the necessary environment and vehicle models. The corrector then determines (using 
multiple predicted trajectories with various control histories) the sensitivities of the final 
state constraints to the control variables. The sensitivities are then used to compute the 
required control variable values to reach the desired final state conditions. Heat rate control 
is provided locally during the regions of maximum heating without significantly affecting the 
assumed control histories. Also, in-flight measurements are utilized to increase the accura- 
cy of the predicted trajectories by compensating for off-nominal conditions. 

Such a simple profile for the maximum downrange and crossrange cases simplifies the 
modeling of the control histories in the predictor. The only remaining question is how much 
of the footprint this profile will capture. As will be seen in Subsection “4.2 Open-Loop 
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Footprint” on page 57, such a profile achieves almost complete coverage of the achievable 
footprint. 

Also of concern is the linearity and convergence properties of the final state constraints 
with the control variables. As will be seen, over almost all of the footprint except near the 
edges, the constraints are highly linear and convergent with the control variables. Opera- 
tionally, only about 75% of the achievable footprint is used to ensure guidance margin. 
Thus, the question of nonconvergence near the edges is avoided. 


32 



3.0 GUIDANCE DESIGN 


3.1 INTRODUCTION 


This section describes the impiementational details of the guidance scheme described 
in the previous section. The equations of motion and environment and vehicle character- 
istics modeled in the predictor algorithm are described. The corrector algorithm to control 
the final state constraints with the two available control variables is derived. Also derived 
are the heat rate control and in-flight measurement algorithms. The heat rate control algo- 
rithm provides control of the peaks in stagnation heat rate during the early portion of entry. 
The in-flight measurement algorithm utilizes accelerations measured by the navigation sys- 
tem to more accurately model the expected environment and vehicle characteristics in the 
predictor algorithm. Because the predictor-corrector algorithm is computationally intensive, 
areas where significant execution time savings have been or can be realized are indicated. 
Program listings of the algorithm coded in the HAL/S computer language are presented in 
“Appendix B. ALGORITHM PROGRAM LISTINGS” on page 135. 

As will be seen, the only inputs to the guidance system are the environment and vehicle 
models, the assumed control profiles, and the navigated state vector. The state vector is an 
input to any guidance system. The other inputs are developed for the analysis of* any new 
vehicle. Therefore, the guidance system is highly transportable between vehicles because 
only the vehicle characteristics and aerodynamics model must be changed for a new vehi- 
cle. 
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3.2 UNIT TARGET VECTOR 


The target aim point to which the vehicle is to be guided is specified by the longitude 
and geodetic latitude of the Terminal Area Energy Management (TAEM) interface point which 
occurs at 80K feet for the Shuttle. This point is selected based on the guidance algorithm 
employed during the TAEM guidance phase. TAEM guidance provides precise control of 
vehicle energy during the final stages of entry to guide to a specified runway with acceptable 
energy. For computational ease, the longitude and geodetic latitude are converted to a tar- 
get unit vector in Earth-fixed coordinates by first computing the geocentric latitude from, 



( 1 ) 


( 2 ) 


The unit target vector is then computed from, 




cos(<j> c ) COS(/l) 
cos (<f> c ) sin(A) 
sin(<£c) 


(3) 


Alternatively, 



/* r 



where, 

i z = sin(</> c ) 

'x = cos(A) 7 1 - /f 

i y = sign(A) J 1 - /? - i 2 z 


( 4 ) 


( 5 ) 

( 6 ) 
( 7 ) 


3.3 COMMANDED ATTITUDE COMPUTATION 


Because the predictor can not be executed as frequently as analytic guidance algo- 
rithms early in the entry, and because it in fact does not have to be executed as frequently, 
it is necessary to update the commands sent to the vehicle autopilot more frequently than 
the predictor-corrector execution rate. Typically, this would be done at the rate of current 
analytic guidance algorithms, e.g., the Space Shuttle rate of .52 hz. The commanded bank 
angle, 4> cmd , is computed for the linear bank with velocity profile as shown in Figure 18 on 
page 95 from the desjred bank angle, <j> d , and the current navigated inertial velocity magni- 
tude, V, , 

4>cn, d = <t> d ( 8 ) 

V EI ~ v f 
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to yield the near-optimal linear bank with velocity profile. The desired angle of attack con- 
trol history is a constant angle of attack, and therefore, 

<*c md = (9) 

As implemented in the current design, the guidance algorithm executive is executed at 1.0 
hz. The attitude commands are updated at this frequency using Eq. (8) and (9) . The pred- 
ictor-corrector algorithm is executed at .02 hz. during the entire entry phase, although it is 
practical to run it much more frequently late in the trajectory when the length of the trajecto- 
ry to be predicted is short. The possible execution rate of the predictor-corrector for a typi- 
cal flight computer is addressed in Subsection “4.8 Algorithm Execution Time” on page 66. 


3.4 CORRECTOR ALGORITHM 


The corrector algorithm is executed to update the commanded attitude control history to 
be flown. The guidance algorithm controls to two final state constraints, downrange error 
and crossrange error, using two control variables, a constant angle of attack and the inter- 
cept of the bank profile at the entry interface velocity as shown in Figure 18 on page 95 and 
expressed in Eq. (8) . 


Expanding the downrange and crossrange errors in a Taylor series expansion of the 
control variables and neglecting the second-order and higher terms yields, 


A DR e 


dDR e 

dcc d 


Aa rf 


+ 


dDR e 


A 4>d + 


( 10 ) 
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. dCR e 3CR e . 

AC/?„ = — r Aa rf + — — — A <j> d + ... 


da H 




( 11 ) 


To intercept the target, the change in the constraint errors must null the predicted errors, or, 

A DR e = —DR e . (12) 

A CR e = -CR e (13) 

Equations (10) through (13) provide a set of two simultaneous equations in two unknowns, 


r dDR„ dDR. 


dct d 

dCR„ 


8<f>d 

dCR„ 


l— da d d<f> 0 


Aa rf l 

~—DR~ 

_&4> d ] 

-CR e _ 


(14) 


which are solved for the control variable changes required, 

, 8DR e dCR. , , 

Aa rf = (-n- CR o — DR e ) I det 


S<f>d 


34> d 


. , , dCR e „„ dDR e „ x , , 

A 4> d = (-T— =- DR e — CR.) I det 


3a„ 


d a. 


(15) 


(16) 


where det is the determinant of the matrix in Eq. (14). The partial derivatives are approxi- 
mated by finite difference equations of the form, 

dDR e DR e (cf> d = <t> 3 ) - DRl<j> d ~4>,) 


d<f>a 


4> 3 


(17) 


There are four partial derivatives that must be evaluated. They can be evaluated from three 
predicted trajectories with control histories selected as: 
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1. a, 


= a'd, = <P'd 

2. a 2 = a'„ + <5a d = <j>’ d 

3. a 3 = a' rf <f> 3 = <f>’ d + 

where the primes denote the control variables from the previous guidance solution. The 
new guidance commands are then, 

= a' d + Aa d (18) 

<f>d = 4>'d + &4>d ( 19 ) 


Protection must be provided for the. case where the determinant in Eq. (15) and (16) is 
small or identically zero which corresponds to a loss of control authority of the control vari- 
ables over the control constraints. In this case, no change is made to the control variables, 
and the guidance command from the previous cycle is used. As the vehicle approaches the 
TAEM interface altitude, the control authority decreases. Large control variable changes 
become necessary to null the constraint errors in the short flight time remaining. This prob- 
lem can be avoided in one of two ways. First, the guidance commands can be frozen at a 
selected point before the termination altitude. For entry guidance, this approach is not pre- 
ferred because the vehicle still has not landed. Alternatively, the target aim point can be 
lowered below the TAEM interface altitude point at which TAEM guidance is activated. The 
decreasing control authority problem is therefore reduced. 

For the simulated trajectories in this report, the first approach is employed because it is 
desired to evaluate guidance performance by considering the dispersions in the final state at 
the TAEM interface altitude. Because the guidance algorithm controls only the final state 
and not the intermediate states, it is necessary to target for the point at which the guidance 
is terminated. 
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3.5 PREDICTOR ALGORITHM 


3.5.1 Introduction 

The predictor algorithm is a simplified three-degree-of-freedom (3-DOF) trajectory simu- 
lator complete with models for those environment and vehicle characteristics necessary to 
model the- translational equations of motion of the vehicle. Because the predictor is compu- 
tationally intensive, the algorithm must be carefully designed to minimize computation, and 
the coding of the algorithm in a particular computer language should make use of any lan- 
guage-specific features to reduce computational requirements. Also, because the corrector 
only utilizes the final state vector errors to correct the control variables, only the accuracy of 
the predicted final state vector need be considered in selecting those effects to be modeled. 

The environmental effects of concern for the long trajectories flown by entry vehicles 
over large altitude and velocity ranges are: 

1. Variation of atmospheric properties with altitude 

2. Earth oblateness effect on gravity vector 

3. Effect of atmospheric rotation with Earth on relative velocity vector 

4. Movement of runway due to Earth rotation 

The vehicle characteristics of importance are: 

1. Vehicle mass 

2. Aerodynamic coefficient variation with flight regime 

3. Aerodynamic coefficient variation with angle of attack 

4. Control history during trajectory 

Dispersions to be considered are: 
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1. Vehicle mass variation from nominal 

2. Winds 

3. Atmospheric density variation from nominal atmosphere 

4. Aerodynamic coefficient variation from nominal 

These dispersions can be measured in-flight because they affect the sensed acceleration 
measured by the vehicle's inertial navigation system. The estimation of these dispersions is 
discussed in Subsection “3.6 Estimators 5 ’ on page 48. 

The predictor performs the following computations upon being called by the corrector 
with a desired control variable history: 

1. Initialize the predictor state to the navigated state vector 

2. Compute any ancillary parameters from the state vector 

3. Compute the total acceleration vector from the predictor state vector and the envi- 
ronment and vehicle models using the control variable profiles specified by the cor- 
rector 

4. Integrate the equations of motion forward in time one time step 

5. Check the predictor termination conditions 

a. Repeat steps 3 and 4 if the conditions are not met 

b. Continue on to step 6 if the conditions are met 

6. .Compute and return to the corrector the final predicted state errors from the target 
state vector and the predicted final state vector 

3.5.2 Equations of Motion 

The corrector provides a time-homogeneous navigated state vector comprised of, 

1. The GMT time tag of the state vector, t GMT 

2. The inertial position vector, R, 

3. The inertial velocity vector, V ( 
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Also provided is the control variable history to be followed for the prediction. The equations 
of motion to be integrated are, 

dR, 


d t 


= V, 


( 20 ) 


d V, 
dt 


= A 


( 21 ) 


The acceleration is computed from the atmosphere and vehicle models as follows, 

A = -jfi r — 9i + 3 aero (22) 


The gravitational acceleration, g lt is computed including the J 2 term as, 


a = ig 

\ R ,\ 2 

(23) 

where, 


X = X + i j, R r“" «i s z‘)X + 2 z u> 

2 \R, I 2 

(24) 

and, 


Z — l R , * Ipole 

(25) 

The aerodynamic acceleration, a Wf01 is computed from, 


^ aero ^fift hi ft & drag ^drag 

(26) 

where, 


__ C L q S 
a„H - m 

(27) 
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(28) 


&drag 
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V* = 


Vs* Vi, 


V/? — V, — oj eanh x R, 




' drag 




',« = (VaoX'/J cos(</>) + i lat sin(</>) 


h, t = 


l„ — 


'r X 'drag 


I 'r X 'drag I 


IK/ 1 


The acceleration due to lift, a, (ft , is more easily computed from, 


q ^ drag 


(29) : 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 


since the nominal lift-to-drag ratio, L/D, is corrected using in-flight accelerometer measure- 
ments of the actual vehicle sensed aerodynamic accelerations. 

The atmospheric density, p, is computed by the atmosphere model using the position 
vector, R t . The 1962 U.S. Standard Atmosphere model is employed and is described in Ref- 
erence [16], If another atmosphere model is selected as being a more accurate estimate of 
the day-of-flight atmosphere, this model would replace the 1962 U.S. Standard Atmosphere 
model. An operational vehicle might employ monthly or seasonal atmospheres from such 
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sources as the GRAM Atmosphere [9] or even day-of-flight measurements to more accu- 
rately model the expected atmosphere in the predictions. The level of accuracy required in 
the atmosphere model will depend on the vehicle ranging capability and the amount of that 
capability to be used for a particular entry. Entries to the edges of the footprint will demand 
a very accurate atmosphere model. 

The aerodynamic coefficients are highly vehicle dependent. To minimize computational 
requirements, they should be updated during the prediction as infrequently as possible. Of 
course, the update frequency required depends on the trajectory flown and the rate of 
change of the aerodynamic coefficients with flight regime change. The aerodynamic coeffi- 
cient model for the ERV is presented in "Appendix A. ERV AERODYNAMICS MODEL” on 
page 133. 

The density, p, from the atmosphere model and the lift-to-drag ratio, L/D, from the aero- 
dynamic model are both corrected by in-flight measurements as covered in Subsection “3.6 
Estimators” on page 48. The estimated dispersions are compensated for using the following 
equations, 


P = K p P« d (37) 



where the density and lift-to-drag ratio scale factors, K a and Kj__, are provided by the estima- 
tor and are held constant throughout the prediction being made. 

The control history to be followed is the constant angle of attack, a rf , and the linear bank 
angle with velocity, <p. The latter is computed from, 
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<t> = 4>d 


(39) 


V,-Vf 
V B -V t 

where (j) d is the intercept of the linear bank angle profile at the entry interface velocity, V EI . 
Because the entry interface and final velocities are not known a priori, and because small 
variations in them have little effect on the predicted trajectory compared with the selected 
control variables' values, the velocities are selected as constant values that cover all 
expected dispersions in the entry and final velocities. These values are, 

V B = 26,000 ft/sec 

V f = 1,000 ft/sec 


3.5.3 Integration of the Equations of Motion 


The equations of motion are integrated using the 4th order Runge-Kutta algorithm with a 
variable time step to minimize the number of time steps required to integrate the trajectory 
to the final state. The 4th order Runge-Kutta algorithm requires four evaluations of the 
acceleration per time step, but permits a time step more than four times as large as an algo- 
rithm requiring only one acceleration evaluation per time step. The Runge-Kutta solution 
[17] for the differential equations of motion of the form, 


dR, 

~dT 


v, 


(40) 


dV, 

~di 


f(t , R,, V,) 


(41) 


is, 
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R,(t + At) = R,(t) + AL (k 0 + 2K, + 2K 2 + K 3 ) 
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Vtf + At) = V,(t) + f (K' 0 + 2K<, + -2KV+ K' 3 ) 
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(45) 

(46) 

(47) 
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(49) 

(50) 

(51) 


The time step is varied inversely with the total acceleration on the vehicle. This method 
of time step control was selected because of its simplicity. The time step control equation is 
of the form, 

At = (52) 

14 I 
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and the time step is limited between a minimum and maximum value, 


At = midval (At min , A t, A t m J 


( 53 ) 


The optimization of the integration algorithm is important in developing a flight quality 
algorithm, but is beyond the scope of this study. Higher-order integration algorithms with 
time step control methods [17] may yield significant reductions in the required computation 
time. 


3.5.4 Termination Conditions for the Predictor 

After each integration time step, the predicted state is compared with the termination 
condition. The termination condition is defined by the altitude of TAEM interface (80K feet). 
Because the predicted state at the TAEM interface altitude may have a relatively large alti- 
tude rate and range rate, the predictor must be terminated accurately to provide an altitude- 
homogeneous set of predicted state errors. Also, the variable time step control may allow 
large integration time steps if the acceleration is low near the final state, further complicat- 
ing the task of terminating accurately. Reasonable altitude homogeneity is ensured by forc- 
ing use of the minimum integration time step starting some safe altitude above the 
termination altitude. 


3.5.5 Final State Error Computation 

The final state errors are computed from the unit target vector and the predicted final 
state vector Because the target is fixed to the Earth and moves a significant distance dur- 


46 



ing the long entry trajectory, the rotation of the Earth must be considered. This is done by 
transforming the final state vector from inertial to Earth-fixed coordinates with the rotation 
matrix Mf ' which is computed from the predicted termination time, the known orientation of 
the Earth at some epoch time, and the known rotation rate of the Earth. This computation is 
performed in the Earth-Fixed-From-Reference subroutine of the predictor-corrector which 
may actually be a GN&C utility function also employed by the navigation principal function. 

The downrange and crossrange errors are defined as shown in Figure 19 on page 96. 
The errors are computed by first computing the downrange (in-plane) and crossrange (per- 


pendicular) directions as follows, 
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The downrange and crossrange errors are then, 

DR e = R equator cos -1 (/f x if F pl3ne ) sign((i FF x if F plane ) • i% rpen ) (59) 

CR e = R equator cos ~\if n F plane • /f) sign((if„ F plane x i FF ) • (i E p F rpen x if F p/ane )) (60) 
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These errors have the dimensions of R equator and are converted to nautical miles for ease of 
interpretation. 


3.5.6 Algorithm Coding 

A few comments regarding implementation of the predictor are appropriate. The com- 
putations required to update the aerodynamic coefficients are the major computational load 
for the predictor. It was found that it is not necessary to update the aerodynamics on each 
of the four acceleration evaluations of the 4th order Runge-Kutta algorithm. They are there- 
fore only evaluated once each integration time step. The computational load could be 
reduced further if they are only updated when the independent variables (altitude, viscous 
interaction parameter, and Mach Number) change by a significant amount from the previous 
update. Also, although not done in this implementation, the aerodynamic coefficients should 
be curve-fit if possible to avoid a table lookup and interpolation implementation. It is noted 
in Figures 20 on page 97 and 21 on page- 98 that the aerodynamic coefficients do not change 
very much below 300K feet until the Mach Number decreases below 2, so perhaps, two 
tables or curve-fits would suffice instead of the thirty tables currently used. 


3.6 ESTIMATORS 


The final state predicted by the predictor algorithm for a particular control history is a 
function of the assumed environment and vehicle characteristics. The accuracy of the pre- 
dicted final state can be increased, and hence, the guidance margin increased, if in-flight 
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measurements are utilized to make the assumed models more accurately reflect the condi- 
tions actually experienced by the vehicle. 

The accelerations modeled in the predictor are due to gravity and the aerodynamic forc- 
es. The gravity acceleration can be modeled to sufficient accuracy using standard gravity 
models. However, the aerodynamic accelerations are subject to significant variations due to 
uncertainties in the atmospheric density, atmospheric winds, vehicle aerodynamics, and 
vehicle mass. These uncertainties can be compensated for in the predictor by applying a 
multiplicative scale factor to the lift and drag accelerations modeled in the predictor that is 
equal to the ratio of the actual accelerations experienced to the predicted accelerations at 
any point in the trajectory. 


The measured lift and drag accelerations are derived from the inertial measurement 
system sensed acceleration assuming a zero sideslip angle as follows, 


A 

^drag 



( 61 ) 


A ^ A 2 

^ lift \, df ® 3/ d drag 


( 62 ) 


where the inertial acceleration, a /5 is computed by back-differencing the accumulated sensed 
velocity counts from the inertial measurement unit. 


3 / 


]/jmu ^imu past 

At~. 


( 63 ) 


In the predictor, the aerodynamic accelerations are, 


^drag 


m 2 ^ 


( 64 ) 
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Q a drag 


( 65 ) 


& lift — 


Data from the Shuttle program [10] shows that the primary dispersion affecting the 
aerodynamic acceleration is in the atmospheric density. Further, over large altitude ranges, 
this dispersion can be modeled to an accuracy sufficient for the prediction process as a con- 
stant multiplicative bias. Therefore, for implementational purposes, the dispersion in the 
aerodynamic accelerations due to the atmospheric uncertainties will be lumped into a densi- 
ty scale factor as follows, 

A 



( 66 ) 


where, 



(67) 


and the values for the nominal vehicle characteristics and the nominal atmospheric density 
are determined using the predictor models for the vehicle state at the time of the measure- 
ment. Because the nominal ballistic coefficient is assumed in deriving the measured densi- 
ty, and the measured acceleration is due to the actual ballistic coefficient, uncertainties in 
the ballistic coefficient will be reflected in the measured density. The equation for the drag 
acceleration in the predictor is then, 

ddrag = (~Trr] L, T K P Pstd (68) 


or substituting for K ^ from Eq. (66) yields, 


^drag 



( 69 ) 


Substituting forp from Eq. (67) then yields, 
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A 

^drag 


(70) 


^drag 


so the modeled drag is corrected for the dispersed drag coefficient, density, relative velocity, 
and vehicle mass. 


In general, the measured drag acceleration is a noisy signal and will exhibit short term 
variations due to short lived local atmospheric dispersions [10]. Filtering of the density 
scale factor is therefore necessary and is implemented using a first-order filter, 

A 

K p = (1-K,) K p + K, (71) 


which has a time constant , x p , of, 


Af 

ln(1 - K,) 


(72) 


where Af is the sample rate of the measured drag acceleration, and K 1 is the filter gain. A 
similar lift-to-drag ratio scale factor is derived and applied to the lift acceleration, 



where, 



(74) 
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and, 


a a 

d drag 

Again, filtering is necessary, 


( 75 ) 
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(76) 


yielding a time constant, of, 

o 

T = Af 

7 ln(1 - K z ) 


(77) 


A time constant of 25 seconds was selected for both the density and L/D filters. This value 
filtered out the high frequency density shear components seen in the Shuttle profiles while 
still providing adequate response to long term disturbances. 


3.7 HEAT RATE CONTROL 


The primary trajectory constraint on entry vehicles is the maximum heat rate the vehicle 
can withstand. In general, the thermal protection system material is selected to withstand 
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the maximum local heat rate on any particular portion of the vehicle, and the material thick- 
ness is selected to withstand the total integrated heat load c/er the trajectory. Accurate 
pre-flight predictions of the expected heat rate during entry can significantly reduce the ther- 
mal protection system weight yielding significant performance increases for an entire mis- 
sion. 

Inspecting the reference trajectories in Figure 16 on page 93 shows that sharp peaks in 
the heat rate occur. If these peaks are accurately controlled, and this control can be accom- 
plished using only short term departures from the predictor assumed control history, no sig- 
nificant departure will occur from the desired trajectory. 

Heat rate control can be accomplished using either angle of attack, bank angle, or a 
combination of both. Of these, bank angle alone is preferred because a constant angle of 
attack trajectory is assumed and because angle of attack changes the vehicle drag coeffi- 
cient resulting in a rapid change in energy rate and a rapid departure from the desired tra- 
jectory. Also, most entry vehicles restrict the angle of attack range during maximum heat 
rate regions to reduce the area on the vehicle that must be protected from the high heat 
rate. Although the ERV does not need to restrict the angle of attack range, and hence, the 
guidance does not provide for such a capabilty, the restriction can be handled by replacing 
the constant angle of attack control history by a reference angle of attack control history 
about which a constant angle of attack bias is applied for control. 

Heat rate control is accomplished by computing the incremental bank angle required to 
fly along the specified heat rate boundary (assumed to be a constant heat rate for any flight 
regime) and then modulating bank* angle according to the guidance value or the guidance 
value plus the incremental lift for heat rate control, whichever requires more lift up. Hence, 
no effort is made to pull the vehicle down into the atmosphere to follow the heat rate bound- 


53 



ary; instead, lift up is applied if the vehicle is flying "too low". The incremental lift for heat 
rate control is computed to provide a second-order control response as follows, 

K .. 

cos(A 4>) = -=T- (Q - Q des ) + -J- (Q - 0| im ) (78) 

q q 

To fly along a constant heat rate boundary, 

Qum = constant (79) 


and the desired rate of change of heat rate, Q des , is, 
Qdes ~ 0 


(80) 


so, 

K„ . 

cos(A <f>) = -±-Q + (Q - Q, im ) (81) 

q q 

The stagnation heat rate is determined using the Engineering Correlation Formula [18] for a 
one foot radius reference sphere as, 
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The time rate of change of heat rate, Q, is determined by back-differencing the heat rate 
between guidance cycles, 


Q = 


Q Opast 
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(83) 


The equations of motion assuming small flight path angle yield, 
C, a S 

h = — - — cos (</>) - g 


(84) 
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Considering only the perturbations due to the incremental lift, cos(A </>) , from Eq. (81) yields, 


h - 
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Qlim) = 0 


(85) 


Proper selection of the gains K ^ and K j is accomplished by linearizing- Eq. (85) in altitude 
and assuming that the time rate of change of V R is small compared to the change in V/T. 
With these assumptions. 


Q 


17700 


( V h V 05 

V 10000 ) 


d 'fp 
d h 


h 


and, 


Q 


17700 


( v* y« 

V 10000 ) 


d_h 

d h d t 


( 86 ) 


(87) 


Therefore, the homogeneous second-order differential equation in altitude is, 

h + KK-h + KK-h = 0 (88) 


where, 
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The natural frequency and damping ratio of the second-order differential equation are, 


JK K g 

(90) 
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or alternatively, for a desired natural frequency and damping ratio, K ^ and K - are selected 
as, 
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(92) 


Ki = 


Kv = 


K 

K 


(93): 


The derivative in Eq. (89) can be evaluated assuming an exponential atmosphere of the form, 
P = e - ^) (94) 


yielding, 

e ~ ( h / 2 "s) (95) 

d h 2 h s 

This logic is contained in the guidance algorithm in the Heat Rate Control subroutine. The 
incremental lift required for heat rate control is provided to the Attitude Command subrou- 
tine which adds it into the guidance command if it requires more lift up than the guidance 
command. This occurs when the incremental lift given by Eq. (81) is greater than zero, 

cos(A <j>) > 0 (96) 

Appropriate values of the natural frequency and damping ratio were determined parametri- 


caily as, 



C O n = 

0.10 rad 
sec 

(97) 

C = 

1.00 

(98) 
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4.0 PERFORMANCE 


4.1 SIMULATOR 


Open-loop and closed-loop entry trajectories were simulated for the Entry Research 
Vehicle (ERV) using a derivative of the 6-DOF Aeroassist Flight Experiment Simulator (AFES- 
IM) [19] developed at The Charles Stark Draper Laboratory which is coded in the HAL com- 
puter language. For this study, the aerodynamic model described in “Appendix A. ERV 
AERODYNAMICS MODEL.” on page 133 and the wind model shown in Figure 2 on page 79 
were incorporated into the AFESIM. The characteristics of the ERV [14] are listed in Table 1 
on page 73. The entry conditions with which all trajectories were initialized are also listed in 
Table 1 on page 73. Because only the performance characteristics of the guidance were 
being evaluated, the simulator was operated in the 3-DOF mode. 


4.2 OPEN-LOOP FOOTPRINT 


Open-loop trajectories were run using the constant angle of attack and linear bank with 
velocity profiles to determine the portion of the footprint achievable. All trajectories were 
terminated at the TAEM interface altitude of 80K feet, so the footprint can be increased about 
100 nautical miles in all directions due to the range flown below 80K feet. 
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Figure 22 on page 99 shows the lift-to-drag ratio, L/D, for the ERV at Mach 10 versus 
angle of attack, a. It is seen that maximum L/D is obtained at an angle of attack of 15 
degrees. It is desirable to fly on the back side of the L/D curve (angle of attack greater than 
15 degrees) so as to maximize the drag coefficient for a given L/D. This reduces heating by 
causing a quicker loss of velocity early in entry than flying at the same L/D on the front side 
of the L/D curve. The L/D versus angle of attack curve shows the same shape with the max- 
imum L/D at 15 degrees for all flight regimes with only a variation in the magnitude of L/D 
across the angle of attack range. Therefore, angle of attack is modulated between 15 and 50 
degrees for the footprint with 15 degrees corresponding to maximum L/D and 50 degrees 
corresponding to minimum L/D. 

The open-loop footprint is shown in Figure 23 on page 100. Also shown for comparison 
is the reported footprint for a heat load limit of 175K BTU/sq ft (shown earlier in Figure 12 on 
page 89). That footprint included the range flown below 80K feet, hence the slight differ- 
ences. It is seen that almost the entire reportedly achievable footprint is captured with the 
assumed control profile. Most importantly, all of the maximum crossrange region is reached 
when the range flown below 80K feet is included. Also, most of the minimum downrange 
region of the footprint was captured even though the control profiles used do not correspond 
to the optimal profiles determined using POST and shown in Figures 14 on page 91 and 15 
on page 92. Additionally, the footprint reported in Reference [14] does not include the large 
area in the minimum downrange region that the open-loop trajectories reached. The small 
area not reached in the minimum downrange region by the control profiles is relatively 
unimportant because the downrange ranging capability of the vehicle can be adjusted by 
changing the deorbit time. A vehicle in low earth orbit travels at about four nautical miles 
per second, so downrange is easily adjusted while on-orbit. 
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Because the predictor-corrector guidance algorithm will follow the same control histo- 
ries as used to generate the open-loop footprint for a nominal trajectory, the guidance algo- 
rithm can reach all of the open-loop footprint for nominal conditions. It is seen that the 
achievable footprint is bounded by the heat rate and heat load limits imposed on the ERV. 
At least an additional 2000 nautical miles of ranging capability in the downrange direction 
exists if the heat limits are relaxed. 


4.3 EFFECT OF DISPERSIONS ON FOOTPRINT 


The effect of dispersions on the achievable footprint was determined by repeating the 
open-loop trajectories with the dispersions discussed in Subsection “2.2 Dispersions” on 
page 26. The worst-case (3<r) dispersions are summarized in Table 2 on page 73. Table 3 
on page 74 shows the dispersions in downrange and crossrange for three of the control his- 
tories in the maximum downrange region of the footprint. It is seen that only variations in 
the lift and drag coefficients cause significant dispersions in the final state. Also, it is seen 
that the effect of a +10% C L dispersion is the same as that of a -10% C D . dispersion. This is 
expected because both dispersions cause the same increase in the vehicle L/D. The same 
occurs for a -10% C L dispersion and a +10% C 0 dispersion, both of which decrease the 
vehicle L/D. 

The effects of the dispersions on trajectories to the maximum crossrange region of the 
footprint are seen in Table 4 on page 75. Again, it is seen that aerodynamic dispersions 
have the greatest effect. A dispersion that increases L/D increases the range, while a dis- 
persion that decreases L/D decreases the range. 
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Table 5 on page 76 shows the effects of the dispersions on the minimum downrange 
region of the footprint. The worst-case range dispersions again occur for the aerodynamic 
dispersions. 


4,4 ESTIMATOR PERFORMANCE 


Figure 24 on page 101 shows the time response of the density filter with a 25 second 
time constant for the STS-9 atmosphere. This trajectory also has dispersions of +1.9% in 
C D , -3.2% in mass, and a 63.8% crosswind. Therefore, the filter output does not follow the 
actual density dispersion also shown in the figure. When the acceleration level is below 0.07 
g's, the measurements are not incorporated, so the filter is inactive before 300 seconds and 
from 600 to 850 seconds. As the velocity drops, the wind becomes a greater contributor to 
the measured density error, hence the divergence in the measured density ratio starting at 
1000 seconds. Figure 25 on page 102 shows the response of the L/D filter with a 25 second 
time constant for a -1.9% C L and a +1.9% C D dispersion. Again, the winds affect the meas- 
urement by creating errors in the navigated angle of attack, so the estimated L/D ratio is 
slightly in error. 

The use of an air data system like the Shuttle Entry Air Data System (SEADS) could sig- 
nificantly improve the estimation process by providing accurate estimates of the angle of 
attack, atmospheric density, and wind magnitude and direction. More accurate estimates 
will increase the guidance margin, thereby increasing the achievable footprint for dispersed 
trajectories. 
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4.5 CLOSED-LOOP PERFORMANCE 


Based on the results of the open-loop trajectories with dispersions, worst-case disper- 
sions were selected for each of three regions of the footprint: maximum downrange, maxi- 
mum crossrange, and minimum downrange. Closed-loop trajectories with the 
predictor-corrector guidance algorithm were then run to the three regions of the footprint. 
The three- target points selected for the closed-loop performance evaluation are shown in 
Figure 23 on page 100. The 3<r errors defined in Table 2 on page 73 were scaled such that 
the total error due to multiple error sources would still represent a 3cr dispersion so as to 
test the guidance system for reasonably probable dispersion cases [20]. To run all disper- 
sions at their 3<x levels would be unrealistic. 

The nominal and dispersed results for trajectories to each of the three regions are listed 
in Tables 6 on page 77 through 8 on page 77 . Plots of selected parameters from these 
cases are included. Figures 26 on page 103 through 31 on page 108 present the altitude, 
velocity, heat rate, heat load, downrange, and crossrange time histories for the nominal 
maximum downrange trajectory. Figures 32 on page 109 through 37 on page 114 present the 
altitude, velocity, heat rate, heat load, downrange, and crossrange time histories for the 
nominal maximum crossrange trajectory. Figures 38 on page 115 through 43 on page 120 
present the altitude, velocity, heat rate, heat load, downrange, and crossrange time histories 
for the nominal minimum downrange trajectory. In each of these cases, it is seen that the 
heat rate does not approach the heat rate limit, so no incremental bank angle is needed for 
heat rate control. 

The control histories for the nominal maximum downrange trajectory and the dispersed 
case listed second in Table 6 on page 77 are presented in Figure 44 on page 121 and 
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Figure 45 on page 122 . Figure 44 shows the angle of attack histories for the nominal and 
dispersed maximum downrange cases. It is seen that a two degree change in angle of 
attack is required early in the trajectory increasing to four degrees by the end of the trajec- 
tory. Figure 45 shows the bank angle histories for the nominal and dispersed maximum 
downrange cases. The bank angle required shows no change from zero degrees for this 
case. 

The control histories for the nominal maximum crossrange trajectory and the dispersed 
case listed second in Table 7 on page 77 are presented in Figure 46 on page 123 and 

Figure 47 on page 124 . Again, it is seen that a four degree change in angle of attack is 

required for the dispersed case. The bank angle history shows no change for the dispersed 
case from that of the nominal case. 

The control histories for the nominal minimum downrange trajectory and the dispersed 
case listed second in Table 8 on page 77 are presented in Figure 48 on page 125 and 

Figure 49 on page 126 . For this case, approximately a one degree change in angle of 

attack is required. No change is required in the bank angle profile. 

The required change in angle of attack for each of the dispersed cases shown was pri- 
marily due to the change in the vehicle L/D as this was shown to be the primary dispersion 
source in Subsection “4.3 Effect of Dispersions on Footprint” on page 59. The breaking point 
of the guidance occurs when the vehicle does not have enough L/D range to overcome the 
loss in L/D due to aerodynamic dispersions. As mentioned previously, the Shuttle entry gui- 
dance algorithm was required to guide to the TAEM interface aim point to within 2.5 nautical 
miles of position. The results presented for the nominal and dispersed cases show that this 
requirement is met with the predictor-corrector guidance algorithm. Also, the trajectory 
plots show that the algorithm achieves this performance with very infrequent guidance 
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updates (.02 hz.) and with very small control variable changes from the nominal constant 
angle of attack and linear bank with velocity profiles. Most importantly, almost all of the 
achievable footprint is captured using the predictor-corrector algorithm. 


4.6 HEAT RATE CONTROL PERFORMANCE 


The closed-loop trajectories shown previously did not require heat rate control because 
the maximum heat rate experienced was significantly lower than the limit imposed on the 
ERV. The time responses for bank angle, angle of attack, and heat rate for the beginning of 
a typical trajectory with and without heat rate control are shown in Figures 50 on page 127 
through 52 on page 129. The resulting bank angle versus velocity profile is shown in 
Figure 53 on page 130. These trajectories are for the middle of the footprint where the peak 
heat rate does not exceed the limit for the ERV. Therefore, for illustrative purposes, the heat 
rate limit was reduced to 100 BTU/sq ft/sec. Comparing the trajectories with and without 
heat rate control, it is seen that the heat rate control takes place over a fairly long time 
range, but requires a significant departure from the linear bank profile over only a very short 
velocity range. The impact on the trajectory is therefore small, and the predictor-corrector 
stays converged on almost the same control history even though the vehicle does not follow 
the assumed control profile during the heat rate control area. 
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4.7 OVERCONTROL 


For those trajectories not at the edge of the footprint, excess vehicle capability exists 
that can be utilized to increase guidance margin for dispersions that may occur later in the 
trajectory. For example, a 13,800 nautical mile downrange trajectory for the ERV only 
requires flying at 20 degrees angle of attack instead of 15 degrees for the nominal trajectory. 
The ERV can modulate angle of attack between 15 degrees (maximum L/D) and 50 degrees 
(minimum L/D) on the back side of the L/D curve, so the modulation capability is not equally 
centered about the commanded angle of attack if flying at 20 degrees. By flying at 15 
degrees (maximum L/D) early in the trajectory, guidance can center the remaining guidance 
capability equally about the aim point to cover dispersions in all directions, not just those 
that require less L/D to reach the target point. This approach is referred to here as overcon- 
trol or command biasing. 

Overcontrol can be implemented in several ways. First, the command can be biased 
from the desired command when that command is not in the center of the modulation range. 
As the vehicle flies a biased angle of attack, for example, the predicted final state will differ 
from that for the unbiased command in such a direction that the next guidance command will 
be moved in the direction opposite to the bias. By biasing in the proper direction, the com- 
mand can be driven toward the center of the modulation range. If the guidance requires an 
L/D higher than that in the middle of the L/D range, flying at an even higher L/D will drive 
the required L/D toward the middle. Secondly, the target aim point can be moved from the 
nominal aim point early in the entry. For example, for a trajectory to the maximum down- 
range region of the footprint, the target aim point can be moved even farther downrange. Of 
course, at some point in the trajectory, the aim point must be moved back to the desired 
point. 
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The first approach was implemented in the predictor-corrector algorithm by biasing the 
angle of attack by five degrees when it was more than two degrees away from 30 degrees. 
The biasing was terminated at an inertial velocity of 13,500 feet per second so as to allow 
the guidance to fly the proper control history near the end of the trajectory to reach the tar- 
get aim point. 

Figure 54 on page 131 compares the angle of attack control history for a 13,760 nautical 
mile downrange trajectory with the dispersions used for the closed-loop trajectories shown 
earlier. Without overcontrol, the vehicle misses the target aim point by 19.20 nautical miles. 
This occurs because the wind contribution to the dispersion increases as the vehicle velocity 
drops, so the multiplicative scale factor on density does not properly model this dispersion. 
As the wind contribution increases, a higher L/D is required, and the angle of attack is driv- 
en to 15 degrees or maximum L/D. Because maximum L/D was not utilized earlier in the 
trajectory, the vehicle did not reach the target. Late in the trajectory, the predictor-corrector 
goes unconverged as control authority is exhausted, causing the angle of attack to jump 
between 15 and 30 degrees. By this point, the target aim point was unreachable anyway 
due to the dispersions. 

With command biasing, the commanded angle of attack early in the trajectory is that 
corresponding to maximum L/D or 15 degrees. It is seen that biasing drives the commanded 
angle of attack to 25 degrees once the biasing is terminated at a velocity of 13,500 feet per 
second or a time of 3,700 seconds. Later, when the wind dispersion drives the angle of 
attack toward 15 degrees, there is significant margin remaining, and the angle of attack is 
only driven to 24.5 degrees by the dispersion. With command biasing, the miss distance at 
TAEM interface is only 0.27 nautical miles. Therefore, guidance margin is increased by 
using overcontrol. More of the theoretically achievable footprint is attainable for dispersed 
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cases. An even larger magnitude dispersion could have been handled late in the trajectory 
since the angle of attack was not driven to that for maximum L/D. 

Further work is needed in this area to determine the proper way to utilize overcontrol to 
maximize guidance margin for the expected dispersions. The probability of the various dis- 
persions occurring and the histories of those dispersions along a trajectory must be consid- 
ered. For example, if a "thick" atmosphere is encountered early in the trajectory equal to 
the worst-case expected dispersion, it is highly unlikely that the atmosphere will get "thick- 
er" later in the trajectory. Therefore, it is unnecessary to preserve guidance margin in the 
direction needed to cover a "thicker" atmosphere beyond that already required for the 
expected worst-case atmosphere. Such considerations should be taken into account in the 
design of the overcontrol algorithm. 


4.8 ALGORITHM EXECUTION TIME 


An estimate of the execution time required for the predictor-corrector algorithm was 
made using the execution time estimate feature of the HAL compiler. The estimate is for the 
AP101 Shuttle flight computer. Figure 55 on page 132 shows the execution time required in 
seconds as a function of the time to the TAEM interface point for a maximum downrange tra- 
jectory. It is seen that early in the entry when the trajectory to be predicted is long, the pre- 
dictor requires 43.7 seconds of CPU time. When only 500 seconds to the TAEM interface 
point remains, the required time drops to 4.5 seconds. This figure can also be interpreted as 
the minimum update interval for the predictor-corrector. Also, the guidance command will 
be computed and sent to the vehicle autopilot a period of time after the start of the guidance 
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cycle equal to the required execution time. It is seen that early in the entry, a significant 
delay occurs between the start of the guidance cycle and the computation of the guidance 
command. This delay was not simulated in the closed-loop trajectories, but will have a mini- 
mal effect on the guidance margin because the guidance is not trying to fly a reference tra- 
jectory like the analytic guidance algorithms. The predictor-corrector is numerically 
computing a trajectory that will fly directly to the target aim point. Any error that builds up 
between the start of the guidance cycle and the issuing of the guidance command can be 
nulled easily since the entry is long, and the error will shrink as the delay decreases with 
decreasing time to the TAEM interface point. 

The Shuttle AP101 CPU is the product of early 1970's technology and is significantly 
slower than flight computers that might be employed in future entry vehicles. The 80C86 
CPU for example is two to five times faster than the AP101 CPU, so the execution time 
required shown - in Figure 55 on page 132 can be scaled down by a factor of two to five. 
Computers utilizing parallel processing architecture could predict the three required trajecto- 
ries simultaneously in three CPUs, cutting the required execution time by a factor of three. If 
scaled by a factor of four due to the faster CPU and a factor of three due to parallel process- 
ing architecture, the maximum time required drops to 3.6 seconds, and the time with 500 
seconds remaining to the TAEM interface point drops to 0.4 seconds. The predictor-correc- 
tor is therefore a viable guidance scheme for future entry vehicles. 
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5.0 FUTURE RESEARCH TOPICS AND CONCLUSIONS 


5.1 FUTURE RESEARCH TOPICS 


Several topics for further algorithm development and optimization are discussed. These 

are: 

1. Further reductions in CPU execution time 

2. Use of an air data system for in-flight measurements 

3. Use of overcontrol to increase guidance margin 

4. Control of more than two state constraints 

Optimization of the predictor algorithm and the integration scheme can yield significant 
reductions in execution time beyond that already attained. Simplifying the aerodynamic 
model can yield a great reduction in execution time and an equally important reduction in 
the computer core required. The current model has 30 tables, each with 51 breakpoints over 
the angle of attack range. A curve fit of the aerodynamic coefficients over the angle of 
attack range and the flow regimes would reduce the core required to store the model data 
and the computations required for each lookup. 

The estimator algorithm was shown to be effective in determining the dispersions from 
in-flight measurements. However, the estimator is unable to differentiate between density 
dispersions and atmospheric winds. Figure 24 on page 101 showed that the multiplicative 
density scale factor did not accurately model the wind contribution to the drag acceleration 
because the relative contribution of the wind to the dispersion increases as the vehicle 
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velocity drops late in the trajectory. An air data system could provide an independent meas- 
urement of the atmospheric winds, improving the estimation process and increasing guid- 
ance margin by increasing the accuracy of the predicted trajectories. 

The concept of overcontrol was introduced and shown to be effective for at least one 
dispersed case. Further investigations should be made to determine how much overcontrol 
is optimal for the expected dispersions. It may be possible to use the sensitivities of the 
constraints to the control variables to determine a proper amount of command biasing for 
any particular dispersion at any point in the trajectory. 

Only the downrange error and crossrange error at TAEM interface are controlled in the 
current design. The vehicle energy is not controlled which can allow significant dispersions 
in the ranging capability during the TAEM phase of entry. Approaches include redefining the 
TAEM aim point in terms of a desired energy level or utilizing a third control variable to pro- 
vide control over an energy level constraint. The Space Shuttle makes use of a split rudder 
as a speedbrake to provide a large energy control capability. Such an approach could be 
utilized with the predictor-corrector by computing the sensitivity of the three constraints to 
the three control variables. This would require four predictions instead of the three currently 
needed, but the fourth prediction could be made only during the latter part of entry to clean 
up any dispersions in energy level that occur during the entry due to dispersions. The CPU 
execution time would then increase by one-third over that currently projected when the third 
constraint is controlled. 
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5.2 CONCLUSIONS 


A predictor-corrector entry guidance algorithm has been demonstrated that exhibits 
excellent performance and almost complete coverage of the achievable footprint. This algo- 
rithm employs a simple control variable history to achieve near-optimal guidance for the 
maximum downrange and maximum crossrange trajectories. Explicit heat rate control is 
employed without significantly impacting the achievable footprint. This is achieved because 
unlike previous guidance algorithms that included a long heat rate control phase with no 
active targeting, the proposed algorithm always actively targets to the aim point and only 
controls heat rate in the short high heat rate regions as required. 

The algorithm has been demonstrated to handle atmospheric and aerodynamic disper- 
sions within the capability of the vehicle. The required computer execution time is shown to 
be within the capability of new flight computers. 

Algorithm adaptability is provided through the utilization of in-flight measurements to 
improve the accuracy of the predicted trajectory. Algorithm maintenance is simplified 
because there are no reference trajectories used, and there are a minimum of 
vehicle/mission-specific input parameters (l-loads). Transportability of the algorithm 
between different entry vehicles is provided by eliminating vehicle-specific entry phases 
other than the heat rate control phase which only requires the input of a heat rate limit. The 
guidance algorithm does require a vehicle aerodynamic model, but this is developed in the 
normal vehicle definition phase anyway. 

In summary, an entry guidance algorithm has been developed that achieves near-opti- 
mal performance while maximizing flexibility, adaptability, and transportability. Although 
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more computationally intensive than analytic algorithms, execution of the predictor-corrector 
is within the capability of current flight computers. It is hoped that this guidance approach 
will significantly reduce the development and maintenance costs for new entry guidance sys- 
tems. 
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Table T. Characteristics of the ERV and Trajectory Entry Conditions 


Mass 

186.0 

slugs 

Reference area 

177.6 

sq ft 

Mean Aerodynamic Chord 

25.0 

ft 

Altitude 

400,000.0 

ft 

loertial Velocity 

25,778.843 

ft/sec 

Flight Path Angle 

-0.996 

deg 

Inclination 

28.50 

deg 

Latitude 

-28.071 

deg 

Longitude 

-69.313 

deg 

Vacuum Apogee 

150.0 

n.m. 

Vacuum Perigee 

20.0 

n.m. 


Table 2. Dispersions Used in Performance Study 


Dispersion 

Symbol 

Magnitude 

(%) 

Direction 

(deg) 

Aerodynamics 

Lift Coefficient 

CL+ 
C L~ 

+ 10% 
-10% 


Drag Coefficient 

CD+ 

CD~ 

+ 10% 
-10% 


Vehicle Properties 

Mass 

M+ 

+ 5% 
-5% 


Atmospheric Properties 

Density 

P + 

P 

+ 30% 
-30% 


Tailwind 

TW 

99% 

61.5 deg 

Positive Crosswind 

CW+ 

99% 

151.5 deg 

Headwind 

HW 

99% 

241.5 deg 

Negative Crosswind 

cw~ 

99% 

331.5 deg 
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Table 3. Dispersed Cases for Maximum Downrange Region 


(deg) 

4>d 

(deg) 

Dispersion 

Downrange 

(n.m.) 

Crossrange 

(n.m.) 

15 

0 

NOMINAL 

14817 

497 

15 

0 

CL+ 

16445 

456 

15 

0 

CL~ 

13229 

420 

15 

0 

CD+ 

13242 

422 

15 

0 

CD~ 

16810 

427 

15 

0 

M+ . 

14890 

498 

15 

0 

M~ 

14739 

496 

15 

0 

P + 

14615 

494 

15 

0 

P 

15089 

498 

15 

0 

TW 

14882 

498 

15 

0 

CW+ 

14829 

497 

15 

0 

HW 

14751 

496 

15 

0 

cw~~ 

14802 

497 

20 

0 

NOMINAL 

13849 

463 

20 

0 

CL+ 

15388 

499 

20 

0 

CL~ 

12355 

344 

20 

0 

CD+ 

12384 

347 

20 

0 

CD- 

15711 

491 

20 

0 

M+ 

13909 

466 

20 

0 

M~ 

13785 

460 

20 

0 

P + 

13659 

453 

20 

0 

P 

14105 

475 

20 

0 

TW 

13901 

466 

20 

0 

CW+ 

13857 

464 

20 

0 

HW 

13796 

461 

20 

0 

CW~ 

13838 

463 

25 

0 

NOMINAL 

12116 

321 

25 

0 

CL+ 

13415 

439 

25 

0 

CL~ 

10822 

178 

25 

0 

CD+ 

10858 

183 

25 

0 

CD~ 

13715 

456 

25 

0 

M+ 

13161 

325 

25 

0 

M~ 

12068 

316 

25 

0 

P + 

11950 

305 

25 

0 

P 

12339 

342 

25 

0 

TW 

12155 

324 

25 

0 

CW+ 

12121 

321 

25 

0 

HW 

12076 

317 

25 

0 

CW~ 

12108 

320 
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Table 4. Dispersed Cases for Maximum Crossrange Region 


(deg) 

4> d 

(deg) 

Dispersion 

Downrange 

(n.m.) 

Crossrange 

(n.m.) 

15 

60 

NOMINAL 

8308 

1822 

15 

60 

CL+ 

9035 

2149 

15 

60 

CL~ 

7602 

1506 

15 

60 

CD+ 

7605 

1531 

15 

60 

CD~ 

9200 

2193 

15 

60 

M+ 

8342 

1825 

15 

60 

M~ 

8271 

1819 

15 

60 

P + 

8197 

1823 

15 

60 

P 

8465 

1820 

15 

60 

TW 

8260 

1737 

15 

60 

CW+ 

8473 

1890 

15 

60 

HW 

8313 

1866 

15 

60 

cw~ 

7984 

1692 ! 

20 

60 

NOMINAL 

7791 

1661 

20 

60 

CL+ 

8476 

1966 

20 

60 

CL~ 

7125 

1370 

20 

60 

CD+ 

7132 

1392 

20 

60 

CD- 

8622 

2007 

20' 

60 

M+ 

7820 

1664 

20 

60 

M~ 

7761 

1659 

20 

60 

P + 

7689 

1662 

20 

60 

P 

7934 

1660 

20 

60 

TW 

7766 

1607 

20 

60 

CW + 

7846 

1686 

20 

60 

HW 

7775 

1673 

20 

60 

CW~ 

7580 

1572 

25 

60 

NOMINAL 

6940 

1344 

25 

60 

CL+ 

7539 

1602 

25 

60 

CL~ 

6355 

1100 

25 

60 

CD+ 

6363 

1119 

25 

60 

CD~ 

7661 

1635 

25 

60 

M+ 

6963 

1345 

25 

60 

M~ 

6915 

1342 

25 

60 

P + 

6848 

1344 

25 

60 

P 

7070 

1342 

25 

60 

TW 

6933 

1317 

25 

60 

CW+ 

6958 

1355 

25 

60 

HW 

6921 

1347 

25 

60 

cw~ 1 

6819 

1290 










Table 5. Dispersed Cases for Minimum Downrange Region 


(deg) 

<Pd 

(deg) 

Dispersion 

Downrange 

(n.m.) 

Crossrange 

(n.m.) 

30 

90 

NOMINAL 

2982 

938 

30 

90 

CL+ 

3014 

1102 

30 

90 

CL~ 

2934 

778 

30 

90 

CD+ 

2914 

793 

30 

90 

CD- 

3045 

1119 

30 

90 

M+ 

2996 

937 

30 

90 

M~ 

2969 

938 

30 

90 

P + 

2914 

943 

30 

90 

P 

3078 

930 

30 

90 

TW 

2993 

914 

30 

90 

CW+ 

2992 

941 

30 

90 

HW 

2996 

939 

30 

90 

cw~~ 

2982 

907 

30 

80 

NOMINAL 

3851 

1041 

30 

80 

CL+ 

4011 

1233 

30 

80 

CL- 

3680 

858 

30 

80 

CD+ 

3668 

875 

30 

80 

CD~ 

4060 

1254 

30 

80 

M+ 

3865 

1041 

30 

80 

M~ 

3865 

1041 

30 

80 

P + 

3778 

1046 


80 

P 

3954 

1034 

30 

'80 

TW 

3856 

1020 

30 

80 

CW+ 

3841 

1044 

30 

80 

HW 

3834 

1042 

30 

80 

CW~ 

3828 

1007 

30 

70 

NOMINAL 

4932 

1075 

30 

70 

CL+ 

5262 

1279 

30 

70 

CL~ 

4602 

881 

30 

70 

CD+ 

4600 

898 

30 

70 

CD~ 

5337 

1304 

30 

70 

M+ 

4949 

1075 

30 

70 

M~ 

4915 

1075 

30 

70 

P + 

4854 

1078 

30 

70 

P 

5043 

1069 

30 

70 

TW 

4934 

1057 

30 

70 

CW+ 

4925 

1077 

30 

70 

HW 

4914 

1074 

30 

70 

CW~ 

4884 

1040 
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Table 6. Maximum Downrange Region Closed-Loop Results 


Dispersions 

Final State 

CL 

CD 

Mass 

P 

Wind 

<t> 9 

X 

Error 

<%) 

(%) 

(%) 

(%) 

(%) 

(deg) 

(deg) 

(n.m.) 

TARGET POINT 

+ 9.800 

+ 144.700 

- 



NOMINAL 


+ 9.799 

+ 144.702 

0.13 

- 1.9 

+ 1.9 

-3.2 

+ 19.1 

63.8 HW 

+ 9.799 

+ 144.701 

0.08 

- 1.9 

+ 1.9 

-3.2 

STS 1 

63.8 HW 

+ 9.803 

+ 144.694 

0.30 

- 1.9 

+ 1.9 

-3.2 

STS 9 

63.8 HW 

+ 9.802 

+ 144.696 

0.27 

- 1.9 

+ 19 

-3.2 

STS 11 

63.8 HW 

+ 9.801 

+ 144.698 

0.13 


Table 7. Maximum Crossrange Region Closed-Loop Results 


Dispersions 

Final State 

CL 

CD 

Mass 

P 

Wind 

<f> 9 

X 

Error 

<%) 

(%) 

<%) 

(%) 

(%) 

(deg) 

(deg) 

(n.m.) 

TARGET POINT 

- 2.300 

+ 55.000 

- 



NOMINAL 


- 2.301 

+ 55.000 

0.06 

- 1.9 

+ 1.9 

-3.2 

+ 19.1 

63.8 CW~ 

- 2.295 

+ 54.998 - 

0.32 

- 1.9 

+ 1.9 

- 3.2 

STS 1 

63.8 ON~ 

- 2.304 

+ 55.001 

0.25 

- 1.9 

+ 1.9 

- 3.2 

STS 9 

63.8 CW~ 

- 2.301 

+ 54.990 

0.60 

- 1.9 

+ 1.9 

- 3.2 

STS 11 

63.8 CW~ 

- 2.299 

+ 55.000 

0.06 


Table 8. Minimum Downrange Region Closed-Loop Results 


Dispersions 

Final State 





<f> g 


X 

Error 



1 — 


(deg) 


(deg) 

(n.m.) 

TARGET POINT 

-14.900 

+ 

9.800 

- 



NOMINAL 


-14.897 

+ 

9.798 

0.22 

1 1.9 

- 1.9 

+ 3.2 -19.1 

63.8 TW 

-14.902 


9.801 

0.13 | 

+ 1.9 

- 1.9 

+ 3.2 STS 1 

63.8 TW 

-14.901 

+ 

9.800 

0.06 

+ 1.9 

- 1.9 

+ 3.2 STS 9 

63.8 TW 

-14.900 

■ + 

9.800 

0.00 

+ 1.9 

- 1.9 

+ 3.2 STS11 

63.8 TW 

-14.899 

4* 

9.799 

0.08 
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Figure 1. Three-View Drawing of the ERV 
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Figure 2. Atmospheric Wind Profile 
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Figure 3. Envelope of Density Profiles Derived from Shuttle Flights 
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Figure 4. STS-1 Density Profile Comparison 
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Figure 5. STS-9 Density Profile Comparison 
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■ METHOD 1 

METHOD 2 

POINT WHERE A* • 90° 



NONDIMENS IONAI VELOCITY. V 

Reproduced from Reference (12) 


Figure 6. Multiphase Bank Angle Program for L/D = 1.5 
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* MULTIPHASE ROLL PROGRAM (JACKSON'S EQUATIONS) 


o TOV 
° TOh 
4 TO ^ 
• TOV 
■TOh 


aaz 

80. 000 FT 
90° 

0.02 

80.000 


*\ 


TOh ■ 80. 000 FT 


MULTIPHASE ROLL PROGRAMS SOLVED ON 
INTEGRATED TRAJECTORY PROGRAM (METHOD 1) 

EULER IAN ROLL ANGLE FROM SLYE'S 
EQUATIONS SOLVED ON INTEGRATED 
TRAJECTORY PROGRAM (METHOD 2) 

STEEPEST DESCENT SOLUTION (METHOD 3) 



CONTINUOUS ROLL 
PROGRAMS 


Reproduced from Reference (12) 


Figure 7. Crossrange Versus Number of Bank Steps 
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Figure 8. Comparison of Optimum Bank Angle Programs 
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Figure 9. Optimum Shuttle Angle of Attack Profile for Maximum Downrange 
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Figure 10. Optimum Shuttle Bank Angle Profile for Maximum Crossrange 
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Figure 11. Optimum Shuttle Angle of Attack Profile for Maximum Crossrange 
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Figure 12. Landing Footprint for the ERV 
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Reproduced from Reference (14) 


Figure 13. Altitude Histories for the Entry Missions of the ERV 


90 





100 


90 


80 


70 

Bank 

60 

angle, 

50 

deg 

40 


30 


20 


10 


/ 

I - I 


0 


/ 


-0- 


\ 


A 


\!\ 

M 

I Vi 


— 0- Maximum downrange 
— >— Maximum crossrange 
Minimum downrange 


! K 

I 

I - 


\ 


\ 


\ 


2 3 

Time, sec 


e- 


■e- 


5 x 10 3 


Reproduced from Reference (14) 


Figure 14. Bank Angle Histories for the Entry Missions of the ERV 
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Figure 15. Angle of Attack Histories for the Entry Missions of the ERV 
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Figure 16. Heat Rate Histories for the Entry Missions of the ERV 
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Figure 17. Heat Load Histories for the Entry Missions of the ERV 
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Figure 18. Bank Angle Versus Velocity Profile 
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Figure 19. Definitions of Downrange and Crossrange Errors 





Figure 20. Predicted Lift Coefficient Profile for the ERV 
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Figure 21. Predicted L/D Profile for the ERV 
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Figure 22. Predicted L/D versus Angle of Attack Profile for the ERV 
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Figure 23. ERV Open-Loop Footprint with the Control Profile 
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Figure 24. Time Response of the Density Filter 
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Figure 25. Time Response of the L/D Filter 
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Figure 26. Closed-Loop Altitude History for the Maximum Downrange Case 
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Figure 29. Closed-Loop Heat Load History for the Maximum Downrange Case 
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Figure 30. Closed-Loop Downrarrge History for the Maximum Downrange Case 



IOOST 
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Figure 31. Closed-Loop Crossrange History for the Maximum Downrange Case 
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Figure 32. Closed-Loop Altitude History for the Maximum Crossrange Case 
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Figure 33. Closed-Loop Velocity History for the Maximum Crossrange Case 
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Figure 34. Closed-Loop Heat Rate History for the Maximum Crossrange Case 


111 





Figure 35. Closed-Loop Heat Load History for the Maximum Crossrange Case 
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Figure 36. Closed-Loop Downrange History for the Maximum Crossrange Case 
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Figure 37. Closed-Loop Crossrange History for the Maximum Crossrange Case 
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Figure 38. Closed-Loop Altitude History for the Minimum Downrange Case 
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Figure 40. Closed-Loop Heat Rate History for the Minimum Downrange Case 
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Figure 41. Closed-Loop Heat Load History for the Minimum Downrange Case 
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Figure 42. Closed-Loop Downrange History for the Minimum Downrange Case 
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Figure 43. Closed-Loop Crossrange History for the Minimum Downrange Case 
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Figure 44. Angle of Attack Comparison for the Maximum Downrange Case 
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Figure 45. Bank Angle Comparison for the Maximum Downrange Case 


122 




Figure 46. Angle of Attack Comparison for the Maximum Crossrange Case 
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Figure 47. Bank Angle Comparison for the Maximum Crossrange Case 
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Figure 48. Angle of Attack Comparison for the Minimum Downrange Case 


125 





126 




Figure 50. Bank Angle Versus Time Comparison for Heat Rate Control 
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Figure 52. Heat Rate Versus Time Comparison for Heat Rate Control 
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Figure 53. Bank Angle Versus Velocity Comparison for Heat Rate Control 
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Figure 54. Angle of Attack Versus Time Comparison with Overcontrol 
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Figure 55. Required Execution Time for the Predictor-Corrector 
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APPENDIX A. ERV AERODYNAMICS MODEL 


The aerodynamics of the ERV were reported in Reference [14], and the. longitudinal per- 
formance coefficients, C L and L/D, are shown in Figures 20 on page 97 and 21 on page 98. 
Figure 22 on page 99 shows a typical L/D versus angle of attack profile. This profile is for a 
Mach Number of 10, but across the flow regimes, the maximum L/D always occurs at an 
angle of attack of approximately 15 degrees. This data was incorporated into the aerodyna- 
mic model of the simulator and into the aerodynamic model of the predictor. It is seen that 
the aerodynamic flow regimes are a function of: 

1. Mach Number, 7W 

2. Viscous Interaction Parameter, V 

3. Altitude, h 

The Mach Number, M, is computed from, 

M = ~ (99) 


where the speed of sound, C s , is computed from, 



The viscous interaction parameter, V, is computed from, 



( 100 ) 


( 101 ) 


where, 
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r xo.5 


c' = (y— ) 

' static 


T statlc + 122.1 X 10-( 5 / r sfaf/c> 
r + 122.1 X 10~< 5 / r '> 


— 11.0 


( 102 ) 


and, 


— = 0.468 + 0.532 -T^+ 0.195 


^“static 


M 2 


(103) 


1 static 


The Reynolds Number, fie, is calculated from, 

Re . £ y » c ~ 


where the coefficient of viscosity for air, ju, is given by, 

P T& 


= 


/ 2 . 
;iatic 


s + r sfat ,> 


(104) 


(105) 
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APPENDIX B. ALGORITHM PROGRAM LISTINGS 


Compiled listings of the flight software principal functions for the predictor-corrector gui- 
dance algorithm as coded for use in the 6-DOF Aeroassist Flight Experiment Simulator 
(AFESIM) follow. The algorithms are coded in the HAL/S computer language. The principal 
functions are: 

1. IL_LOAD - Values for all constants and l-loads 

2. FSW SEQ - Flight Software Sequencer 

3. ORB_NAV - Orbit Navigation Algorithm 

4. AEROGUID - Predictor-Corrector Guidance Algorithm 

At the beginning of each principal function is a description of the function and the 
input/output parameters. At the end of each principal function is a cross reference table list- 
ing the program line at which each variable is referenced or computed. 
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PBBCffDENG F AO* BLANK NOT FILMED 


HAL/S STD 360-24.20 


INTERMETRICS 


INC. 



STMT . SOURCE 

211 Ml UNLOAD: 

211 Ml PROCEDURE | 

C| 

Cl FUNCTION: VALUES OF I-LOADS AND CONSTANTS 

C| INPUTS: NONE 

C| OUTPUTS: ALL I-LOADS AND CONSTANTS LISTED 

C| COMMENTS: NONE 

C| 


C| 

Cl MATH CONSTANTS AND CONVERSION FACTORS 

C 

212 Ml DECLARE PI SCALAR DOUBLE CONSTANT! 3.1415926535897932385 

213 Ml DECLARE OLG_TO_SEC SCALAR DOUBLE CONSTANT! 3600 ) | 


214 Ml DECLARE SEC_TO_DEG SCALAR DOUBLE CONSTANT! 1 / DEG_TO_SEC)| 

215 Ml DECLARE DEG_TO_RAD SCALAR DOUBLE CONSTANT! PI / 180 )» 

216 Ml DECLARE RAD_TO_DEG SCALAR DOUBLE CONSTANT! 1 / DEG_TO_RAD)> 

217 Ml DECLARE StC_TOJlAD SCALAR DOUBLE CONSTANT! DEG_JO_RAD / 3600 )\ 

218 Ml DECLARE RAD_TO_SEC SCALAR DOUBLE CONSTANT! 1 / SEC_TO_RAD If 

219 Ml DECLARE F1_TO_M SCALAR DOUBLE CONSTANT! 0 .3048 i | 

220 Ml DECLARE M_TO_FT SCALAR DOUBLE CONSTANT! 1 / FT_TO_M>* 

221 Ml DECLARE FT_TO_NM SCALAR DOUBLE CONSTANT! FT_TO_M / 1852)* 

222 Ml DECLARE NM_TO_FT SCALAR DOUBLE CONSTANT! 1 / FT_TO_W)l 

223 Ml DECLARE G_TO_FPS2 SCALAR DOUBLE CONSTANT! 9.80665 M_TO_FTU 

224 Ml DECLARE FPS2_TO_G SCALAR DOUBLE CONSTANT! 1 / GJTO_FPS2 >> 

225 Ml DECLARE LBM_TO„KG SCALAR DOUBLE CONSTANT! .45359237 J \ 

226 Ml DECLARE KG_TO_LBM SCALAR DOUBLE CONSTANT! 1 / LBM_TO_KG)» 

227 Ml DECLARE SLUG_TO_KG SCALAR DOUBLE CONSTANT! LBM_TO_KG G_TO_FPS2)j 

228 Ml DECLARE KG_TO_SLUG SCALAR DOUBLE CONSTANT! KG_TO_LBM FPS2_TO_G)> 

229 Ml DECLARE LBF_TO_N SCALAR DOUBLE CONSTANT! 4.4482216152605 

230 Ml DECLARE N_TO_LBF SCALAR DOUBLE CONSTANT!! / LBF_TO_N ) \ 


APRIL 27, 1987 


15:0:48.60 


CURRENT SCOPE 

XL_LOAD 

IL_LOAD 

IL_LOAD 
IL LOAD 
IL LOAD 
IL_L0AD 
IL LOAD 
IL_L0AD 

IL_L0AD 
IL LOAD 
IL_L0AD 

IL_LOAD 

I L_L0AD 

IL_L0AD 

I L^. LOAD 

IL_LOAD 

XL_L0AD 

IL.LOAD 

IL.LOAD 

IL_LOAD 

IL_LOAD 

IL.LOAD 

XL_LOAD 

IL„LOAD 

IL_LOAD 

IL.LOAD 

IL.LOAD 

IL_LOAD 

IL_LOAD 


IL_ LOAD 



CO 

Oo 


HAL/S 

STD 360-24.20 INTERMETRICS, INC. 

APRIL 27, 1987 

15:0:48.60 

STMT 

SOURCE 


CURREOT SCOPE 

ci 



1 IL LOAD 

cl 

FSW SEQ VARIABLES 


I IL LOAD 

cl 



| IL_LOAD 

251 Ml 

AERO_OAP_CNT = 1) 


1 IL_L0AD 

232 Ml 

AERO_DAP_PHS = 0) 


I IL_LOAD 

233 Ml 

AERO_GUID_CNT = 5) 


| IL_L0AD 

234 Ml 

AERO_GUID_PHS = 0) 


I IL_LOAD 

235 Ml 

ORB_NAV_CNT = 5) 


| I L_LOAD 

236 Ml 

ORB_NAV_PHS = 0) ( 


| IL__LOAD 

C| 



| IL LOAD 

cl 

EPOCH DATA 


| IL LOAD 

Cl 



I IL_LOAD 

el 

* 


1 

237 Ml 

EF TO REF AT EPOCH = MATRIX 1+1.0, +0.0, +0.0, +0.0, +1.0, +0.0, 

♦0.0, +0.0, +1.0)) 

| IL LOAD 

sl 

3D0UBLE ,3,3 


1 

238 Ml 

TEPOCH = 0) 


1 IL_LOAD 

cl 



1 IL LOAD 

Cl 

EARTH PHYSICAL PARAMETERS 


| ILJ.OAD 

Cl 



| IL_LOAD 

239 Ml 

EARTH_FLAT = 1 / 298.3) 


| IL_LOAD 

240 Ml 

EARTH_J2 = 1082. 7E-6) 


| IL_L0AD 

El 

3 


1 

241 Ml 

EARTH_MU = ( 3 . 986012E14 ) / !.3048) ) 


| IL_LOAD 

E| 

_ 


1 

242 Ml 

EARTH POLE = VECTOR 1 2.8991969340471790E-3 , -5. 1580036678323420E-5 


| IL LOAD 

sl 

3D0UBLE ,3 


1 

242 Ml 

9 . 9999579598 948170E-1 ) ) 


1 IL_LOAD 

243 M| 

EARTH_R = (6378166.0) / .3048) 


I ILJ.OAD 

244 Ml 

EARTH_RATE = 7. 292114883223324E-5) 


| IL_LOAD 

El 

_ 


1 

245 Ml 

WE_NAV = EARTH_RATE EARTH_POLE ) 


| IL_LOAD 



HAL/S STD 360-24.20 INTERMETRICS, INC 


APRIL 27 » 1987 15:0:48.60 


STMT 


SOURCE 


CURRENT SCOPE 


CO 

CD 


Cl I XL LOAD 

Cj PREDICTOR-CORRECTOR I-LOADS EMISSION-SPECIFIC ) I I L_LOAD 

C| — | IL_LOAD 


C| | IL_LOAD 

C| TARGET AIM POINT FOR MAXIMUM OOWNRANGE CASE I IL.LOAD 
C| OOHNRANGE = 13849 N.M. I IL_LOAD 
Cl CROSSRANGE = 497 N.M. I ILJ.OAD 
C| J XL_LOAD 


C | I IL_LOAD 

C| INITIAL CONDITIONS: I IL_LOAD 

C | t | UNLOAD 

C| LATITUDE -28.071 DEG I IL_LOAD 

C| LONGITUDE -69.313 DEG I IL_LOAD 

C| INCLINATION 28.50 DEG I XL_LOAD 

C| FLIGHT PATH ANGLE -0.996 DEG I IL_LOAD 

C| INERTIAL VELOCITY 25778.843 FT/SEC I IL_LOAD 

Cl ALTITUDE 400000.0 FT I IL_LOAD 

C | | IL_LOAO 


Cl I IL_LOAD 

C| GEODETIC LATITUDE OF TAEM INTERFACE AIM POINT I IL_LOAD 

C| — | IL_LOAD 

246 Ml LAT.TARGET = 3.084 f t IL_LOAD 


C| | IL_LOAD 

C| LONGITUDE OF TAEM INTERFACE AIM POINT I IL_L0AD 

C| | IL_LOAD 

247 Ml LONG_TARGET = 157. 659 f I IL_LOAD 


C| | ILJ-OAD 

C| INITIAL CONTROL VALUES FOR MAXIMUM OOHNRANGE CASE I IL_LOAD 

Cj | IL_LOAD 

248 Ml PHI_EI * 0.0» I IL_LOAD 

249 Ml ALPHA_EI = 20. Of I IL_LOAD 


C| | IL LOAD 

C| ESTIMATOR FILTER GAINS JTAU =25.0 SECONDS) l IL_LOAD 

C| | IL_LOAD 

250 Ml K_RHO_FILTER_GAIN s .03 921 > I IL_LOAD 

251 Ml L_OVER_D_FILTER_GAIN = .03921 > I IL_LOAD 


c | I IL LOAD 

Cl VEHICLE MASS (SLUGS) 1 IL_LOAD 

C| I IL_LOAD 


252 Ml MASS_NAV = 186. Of 


| IL_LOA0 



avoi ii 1 ! 3 

QV01~TI | H0i33MH03-a013ia3Hd SNISn 3amillV 03QNVWW03 30 S3IV0dn 13 

OV01“1I j N33M139 3AI1H33X3 a0133aa03-aOI3I<13Hd 30 SNOI1TI33X3 30 a39mN |3 

QV01 II | S3 

0V01“1I | <0'T = ain30a3V~lQ |H 092 

QV01“1I | I 3 

OVOIJII | 3AI1H03X3 a0133aH03-a013I03ad 30 SN0I1/133X3 N33M139 !N3W3a0NS 3WI1 1 3 

avoi ii | 1 3 

QV01“1I I <VZZT = 33a~S |W 6S2 

QVOl“lI | (3 

avOl“TI | Aa3 30 V39V 33N3a333a S3 

avo~Tn | 13 

avoi“n | <t = 33aow~Aj.iAva9 ;u es z 

avoi _ n l S3 

OV01“1I | W931 2r HUM 1300W A±IAVa9 13 

avo-Tn 1 1 3 

avo-Tn | ‘ZO’O = 33Nvaino~Nna - 9 su lsz 

avo-Tn j 1 3 

0V01 _ 1I | H3NVain9 31VAI13V 01 H3IHM IV 13A31-9 |D 

avoi _ n | 1 3 

avo-Tn l <o-ooooi ♦ wavi~nv = svia - W3VX - nv jw 9sz 

avo-Tn i <0-00009 = W3Vi“nv |w ssz 

avoi"n | <o-oooooi = aircT3Z33aj~nv jw wz 

avo-pii l <o-ooooo<7 = nx3“nv |w ss z 

avoi“u | S3 

QV01 - II | N0IlVa931NI aoj d31S 3NI1 WOWINIW S3SR 13 

OVO-TlI | ONV «S139aVi * 0NVWWO3 S3Z33a3 <Slint) 33NV0in9 H3IHM IV SSOnilllV j3 

QV01*~1I I S3 

avoi - i i 1 S3 

avo-Tn 1 ( Q39NVH3 AmvwaoN low saaiawvavd N9IS30J sovoi-i aoxoaaaoo-aoioioaad 6 3 

OVOI II J 53 

3d03s iNaaano aoanos iwis 

09‘9t» : 0 : SI Z96I ‘Z2 tlddV 3NI ‘ S 3 I U 3 H II H N I 0Z'«-09£ OiS S/1VH 



HAL/S STD 360-24.20 INTERMETRICS » INC. APRIL 27, 1987 15:0:48.60 

STMT. SOURCE CURRENT SCOPE 

C| | IL LOAD 

Cl LINEAR BANK WITH VELOCITY PROFILE CONSTANTS I IL LOAD 

C| | ILJlOAD 

262 Ml PHIJ)ES_MAX = 180.0 ) | ILJ.OAD 

263 Ml PHI_MAX = 90.0* | IL_LOAD 

264 Ml V_FINAL_MAG = 1000.0 5 | IL_LOAD 

265 Ml V_INITIALJ1AG = 26000, Oj | IL_LOAD 

266 Ml V_MAG_C MANGE = V_INITIAL_MAG - V_FINAL_MAG) I ILJLOAD 

Cj t | il_LOAD 

C| CONSTANT ANGLE OF ATTACK PROFILE CONSTANTS I IL LOAD 

C| l IL_L0AD 

267 Ml ALPHA_MAX = 45.0) I ILJ.OAD 

268 Ml ALPHA_MIN =15.0) I IL_LOAD 

•CJ | IL_LOAD 

C| VARIABLE TIME STEP CONTROL CONSTANTS I IL_LOAD 

C| | ILJ.OAD 

269 Ml DELTA_T_PREDJ3AIN = 200.0) | IL_LOAD 

270 Ml DE LT A_T_P RE D_MAX = 20.0) | IL_LOAD 

271 Ml DELTAJT_PRED_MIN =2.0) | IL_L0AD 

C| | ILJLOAD 

C| HEAT RATE CONTROL CONSTANTS j IL LOAD 

C| J IL_L0AD 

272 Ml HS = 23500.0) I ILJLOAD 

273 Ml OMEGA __QD0T = 0.10) I IL_LOAD 

274 Ml QDOT_LIMIT =125.0) I IL_LOAD 

275 Ml RH0_SL =0.002378) I IL_LOAD 

276 Ml 2ETA_QDOT = 1.00) I IL_L0AD 

277 Ml CLOSE IL_LOAD) l IL_LOAD 

**** BLOCK SUMMARY **** 


COMPOOL VARIABLES USED 


ORIGINAL PAGE IS 
OF POOR QUALITY 



HAL/S STD 360-24. 20 


INTERMETRICS 


INC. 


APRIL 27, 1987 


15:0:48.60 


STMT SOURCE * CURRENT SCOPE 

AERO_DAP_CNT*, AERO_DAP_PHS», AERO_GUID_CNT*, AERO_GUID PHS*, ORB_NAV_CNT* , ORB_NAV_PHS* , EF_TO_REF AT EPOCH*, T_EPOCH* 

EARTH_FLAT*, EARTH_J2*, EARTH_MU*, EARTH^POLE*, EARTH_r1«, EARTH_RATE*» HE_NAV», EARTH_RATE , EARTH_P0LE7 LAT TARGET* 

LONGJTARGET*, PHI_EI*, ALPHABET*, K_RHO_FILT£R_GAIN*, L_OVER_dJfILTER_GAIN*, MASS_NAV*, ALT__EXIT*, ALT_FREEZE_GUID* , ALT TAEM* 
ALT_TAEM_BIAS* , ALT_TAEM, G_RUN_GUI DANCE * , GRAVITY_MODE L* , S_REF*> DT_AEROGUID* » GUID_PASS_LIM*, PHI oIs_MAX*, PHI MAX* 

V_FINAL_MAG*» V_INITIAL_MAG* , V_MAG_CHANGE* , V_INITIAL_MAG, V_FINAL MAG, ALPHA MAX*, ALPHA_MIN*, DELTA_T_PRED GAIN* 

DE LT A_J"_PRE D_MAX* , DELTA_T_PRED_MIN* , HS*, OMEGA_QDOT*, QDOT_LIMIT*, RHO_SL*, ZETA_QD0T* 

"■5 


lx 

ro 
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if 
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*«** COMPILATION LAYOUT **** 

ILJ’OOL: EXTERNAL COMPOOLj 
IL_LOAD: PROCEDURE i 


CO 




HAL/S STP 360-2*. 20 


INTERMETRICS , INC 


APRIL 27 , 1987 15:0:*8.60 


SYMBOL « CROSS REFERENCE TABLE LISTING : 

(CROSS REFERENCE FLAG KEY: * = ASSIGNMENT, 2 = REFERENCE » I = SUBSCRIPT USE, 0 = DEFINITION 5 


-U 


DCL 

NAME 

TYPE 

ATTRIBUTES & CROSS REFERENCE 





21 

AERO OAP CNT 

INTEGER 

SINGLE,, 

.ALIGNED, INITIAL 

XREF: 0 

0021 

* 

0231 

NOT REFERENCED 

20 

AERO_DAP PHS 

INTEGER 

SINGLE , 

ALIGNED, INITIAL 

XREF: 0 

0020 

* 

0232 

NOT REFERENCED 

23 

AERO GUID CNT 

INTEGER 

SINGLE , 

ALIGNED, INITIAL 

XREF: 0 

0023 

* 

0233 

NOT REFERENCED 

22 

AERO GUID PHS 

INTEGER 

SINGLE , 

ALIGNED, INITIAL 

XREF: 0 

0022 

* 

023* 

NOT REFERENCED 

97 

ALPHA El 

SCALAR 

SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0097 

* 

02*9 

NOT REFERENCED 

209 

ALPHA MAX 

SCALAR 

SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0209 

* 

0267 

NOT REFERENCED 

209 

ALPHA MIN 

SCALAR 

SINGLE , 

ALIGNED, INITIAL 

XREF: 0 

0209 

* 

0268 

NOT REFERENCED 

79 

ALT EXIT 

SCALAR 

SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0079 

* 

0253 

NOT REFERENCED 

209 

ALT FREEZE GUID 

SCALAR 

SINGLE , 

ALIGNED, INITIAL 

XREF: 0 

0209 

* 

025* 

NOT REFERENCED 

209 

ALT TAEM 

SCALAR 

SINGLE , 

ALIGNED, INITIAL 

XREF: 0 

0209 

* 

0255 

2 0256 

209 

ALT TAEM BIAS 

SCALAR 

SINGLE , 

ALIGNED, INITIAL 

XREF: 0 

0209 

* 

0256 

NOT REFERENCED 

215 

DEG_TO_RAD 

SCALAR 

DOUBLE , 
2 0217 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 0215 2 0216 

213 

DEG TO SEC 

SCALAR 

DOUBLE , 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 0213 2 021* 

209 

DELTA T_PRED_GAIN 

SCALAR 

SINGLE , 

ALIGNED, INITIAL 

XREF; 0 

0209 

* 

0269 

NOT REFERENCED 

209 

D E LT A_T_P RE D_MAX 

SCALAR 

SINGLE , 

ALIGNED, INITIAL 

XREF: 0 

0209 

* 

0270 

NOT REFERENCED 

209 

DELTA_T_PRED_MIN 

SCALAR 

SINGLE, 

ALIGNED, INITIAL 

XREF; 0 

0209 

* 

0271 

NOT REFERENCED 

209 

DT AEROGUID 

SCALAR 

SINGLE , 

ALIGNED, INITIAL 

XREF; 0 

0209 

* 

0260 

NOT REFERENCED 

159 

EARTH FLAT 

SCALAR 

DOUBLE , 

ALIGNED, INITIAL 

XREF: 0 

0159 

* 

0239 

NOT REFERENCED 

160 

EARTH J2 

SCALAR 

DOUBLE, 

ALIGNED, INITIAL 

XREF: 0 

0160 

* 

02*0 

NOT REFERENCED 

161 

EARTH MU 

SCALAR 

DOUBLE , 

ALIGNED, INITIAL 

XREF: 0 

0161 

* 

02*1 

NOT REFERENCED 

162 

EARTH POLE 

3 - VECTOR 

DOUBLE, 

ALIGNED, INITIAL 

XREF: 0 

0162 

* 

02*2 

2 02*5 

163 

EARTH_R 

SCALAR 

DOUBLE, 

ALIGNED, INITIAL 

XREF; 0 

0163 

* 

02*3 

NOT REFERENCED 

16* 

EARTH_RATE 

SCALAR 

DOUBLE, 

ALIGNED, INITIAL 

XREF: 0 

016* 

* 

02** 

2 02*5 

* 

EF_TO_REF_AT_EPOCH 

3X3 MATRIX 

DOUBLE, 

ALIGNED, INITIAL 

XREF: 0 

000* 

* 

0237 

NOT REFERENCED 

22* 

FPS2 TO G 

SCALAR 

DOUBLE, 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 022* 2 0228 

219 

FT_TO_M 

SCALAR 

DOUBLE, 
2 0221 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 0219 2 0220 

221 

FT_TQ_MM 

SCALAR 

DOUBLE, 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 0221 2 0222 

209 

G_RUN_GUI DANCE 

SCALAR 

SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0209 

* 

0257 

NOT REFERENCED 

223 

G_TO_FPS2 

SCALAR 

DOUBLE, 
2 0227 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 0223 2 022* 

165 

GRAVITY_MODEL 

INTEGER 

SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0165 

* 

0258 

NOT REFERENCED 

209 

GUID PASS LIM 

SCALAR 

SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0209 

* 

0261 

NOT REFERENCED 

209 

HS 

SCALAR 

SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0209 

* 

0272 

NOT REFERENCED 

211 

UNLOAD 

PROCEDURE 

XREF: 

0 0211 NOT REFERENCED 





163 

K RHO FILTER GAIN 

SCALAR 

SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0168 

* 

0250 

NOT REFERENCED 

226 

KG TO LBM 

SCALAR 

DOUBLE , 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 0226 2 0228 

228 

KGTOSLUG 

SCALAR 

DOUBLE, 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 0228 




NOT REFERENCED 






171 

L_OVER_D_FI LTE REGAIN 

SCALAR 

SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0171 

* 

0251 

NOT REFERENCED 

209 

LAT_TARGET 

SCALAR 

SINGLE , 

ALIGNED, INITIAL 

XREF; 0 

0209 

* 

02*6 

NOT REFERENCED 

229 

LBF_TO_N 

SCALAR 

DOUBLE , 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 0229 2 0230 

225 

LBM_TOJ<G 

SCALAR 

DOUBLE , 
2 0227 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 0225 2 0226 

209 

LONG_TARGET 

SCALAR 

SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0209 

* 

02*7 

NOT REFERENCED 

220 

M_TO_FT 

SCALAR 

DOUBLE , 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 0220 2 0223 

10 

MASS_NAV 

SCALAR 

SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0010 

* 

0252 

NOT REFERENCED 

230, 

N TO LBF 

SCALAR 

DOUBLE, 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 0230 




NOT REFERENCED 








lVNKMUO 
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DCL NAME TYPE ATTRIBUTES ft CROSS REFERENCE 


222 

NM_TO_FT 

SCALAR 

DOUBLE , 

ALIGNED, 

STATIC, 

CONSTANT 



XREF: 

0 0222 




NOT REFERENCED 







209 

OMEGA qOOT 

SCALAR 

SINGLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0209 

4 

0273 

NOT REFERENCED 

35 

ORB NAV CNT 

INTEGER 

SINGLE , 

ALIGNED, 

INITIAL 

XREFs 0 

0035 

4 

0235 

NOT REFERENCED 

34 

ORB NAV PHS 

INTEGER 

SINGLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0034 

4 

0236 

NOT REFERENCED 

209 

PHI DES MAX 

SCALAR 

SINGLE , 

ALIGNED, 

INITIAL 

XREF: 0 

0209 

4 

0262 

NOT REFERENCED 

98 

PHI El 

SCALAR 

SINGLE , 

ALIGNED, 

INITIAL 

XREF: 0 

0098 

4 

0248 

NOT REFERENCED 

209 

PHI_MAX 

SCALAR 

SINGLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0209 

4 

0263 

NOT REFERENCED 

212 

PI 

SCALAR 

DOUBLE, 

ALIGNEO, 

STATIC, 

CONSTANT 



XREF? 

0 0212 2 0215 

209 

QDOT LIMIT 

SCALAR 

SINGLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0209 

4 

0274 

NOT REFERENCED 

216 

RADJTO_DEG 

SCALAR 

DOUBLE , 

ALIGNED, 

STATIC, 

CONSTANT 



XREF: 

0 0216 




NOT REFERENCED 







218 

RAD_TO_SEC 

SCALAR 

DOUBLE, 

ALIGNED, 

STATIC, 

CONSTANT 



XREF: 

0 0218 




NOT REFERENCED 







209 

RHO_SL 

SCALAR 

SINGLE , 

ALIGNED, 

INITIAL 

XREF: 0 

0209 

4 

0275 

NOT REFERENCED 

174 

S REF 

SCALAR 

SINGLE , 

ALIGNED, 

INITIAL 

XREFs 0 

0174 

4 

0259 

NOT REFERENCED 

214 

SEC_TO_DEG 

SCALAR 

DOUBLE, 

ALIGNED, 

STATIC, 

CONSTANT 



XREF: 

0 0214 




NOT REFERENCED 







217 

SEC TO RAD 

SCALAR 

DOLBLE, 

ALIGNED, 

STATIC, 

CONSTANT 



XREF: 

0 0217 2 0218 

227 

SLUG_TQ_KG 

SCALAR 

DOUBLE, 

ALIGNED, 

STATIC, 

CONSTANT 



XREF: 

0 0227 




NOT REFERENCED 







6 

T_EPOCH 

SCALAR 

DOUBLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0006 

4 

0238 

NOT REFERENCED 

209 

V FINAL MAG 

SCALAR 

SINGLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0209 

4 

0264 

2 0266 

209 

V_INITIAL MAG 

SCALAR 

SINGLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0209 

4 

0265 

2 0266 

209 

V MAG CHANGE 

SCALAR 

SINGLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0209 

4 

0266 

NOT REFERENCED 

176 

HE_NAV 

3 - VECTOR 

SINGLE, 

ALIGNED, 

INITIAL 

XREFs 0 

0176 

4 

0245 

NOT REFERENCED 

209 

ZETA_QDOT 

SCALAR 

SINGLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0209 

4 

0276 

NOT REFERENCED 



ORIGINAL, PAGE 



HAL/S STD 360-24.20 


INTERMETRICS 


INC 


APRIL 27 s 1987 


15:1:9.08 


STMT SOURCE CURRENT SCOPE 

518 Ml FSW_SEQ: I FSH_SEQ 

518 Ml PROCEDURE l \ FSW_SEQ 

Cj j fsh_seq 

C| FUNCTION: EXECUTE FLIGHT SOFTWARE PRINCIPAL FUNCTIONS AT I FSW SEQ 

Cf PROPER RATE AND IN PROPER ORDER WHEN FUNCTIONS ACTIVE I FSW_SEQ 

C| INPUTS: ORB_NAV_ACT - ORBIT NAVIGATION ACTIVE FLAG I FSW_SEQ 

C| AERO_GUID_ACT - PREDICTOR-CORRECTOR ACTIVE FLAG I FSH_SEQ 

C| AERO_DAP_ACT - DIGITAL AUTOPILOT ACTIVE FLAG | FSW_SEQ 

C| OUTPUTS: NONE | FSW_SEQ 

C| COMMENTS: NONE | FSW_SEQ 

Cj | FSW_SEQ 


C| | FSW SEQ 

C| LOCAL VARIABLES | FSH_SEQ 

c| I fshIseq 

519 Ml DECLARE FSW^PASS INTEGER DOUBLE INITIAUO)> I FSW_SEQ 




APRIL 27, 1987 15sl?9.08 

CURRENT SCOPE 

i 

I FSW_SEQ 
| FSW_SEQ 


El . j 

522 Ml IF A E RO_GUI D_ ACT = ON AND MODI FSW_PASS» AERO_GUID_CNT ) = AEROJ3UID_PHS THEN I FSN_SEQ 

523 Ml CALL AEROJ3UID* I FSH„SEQ 

« 

El . j 

524 Ml IF AE RO_D AP_ ACT - ON AND MODI FSW_PASS» AERO_DAP_CNT } = AERO_OAP_PHS THEN I FSH_SEQ 

525 Ml CALL AERO_DAP» I FSW_SEQ 

526 Ml FSM_PASS = FSH_PASS + H I FSH_SEQ 

527 Ml CLOSE FSM_SEQi i FSH_SEQ 

***# BLOCK SUMMARY #*** 

EXTERNAL PROCEDURES CALLED 

ORB NAV > AERO_GUID» AERO_DAP 

COMPOOL VARIABLES USED 


ORB_NAV_ACT » ORB_NAV_CNT» ORB_NAV_PHS„ AERO_GUID_ACT» AERO_GUID_CNT, AERO_GUID_PHS s AERO_DAP_ACT, AERO_DAP_CNT, AERO_DAP_PHS 


HAL/S STD 360-24.20 INTERMETRICS p INC. 

STMT SOURCE 

E| 

520 Ml IF ORB_NAV_ACT = ON AND M0D( FSW_PASS, ORB_NAV_CNT ) = ORB_NAV_PHS THEN 

521 Ml CALL ORB NAVj 


ORIGINAL PAGE IS 

OF POOR QUALITY 



HAL/S STD 360-24. 20 INTERMETRICS , 

**** COMPILATION LAY O.U T **** 

IL_POOL : EXTERNAL COMPOOL ) 

FSH_POOL: EXTERNAL COMPOOL > 

ORB_NAV: EXTERNAL PROCEDURE* 

£ AERO GUID: EXTERNAL PROCEDURE i 

oo 

AERO_DAP: EXTERNAL PROCEDURE » 

FSH_SEQ: PROCEDURE j 


INC., APRIL 27, 1987 !5:1;9.08 
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APRIL 27, 1987 15:1:9.08 


SYMBOL £ CROSS REFERENCE TABLE LISTING: 


(CROSS REFERENCE FLAG KEY: 4 = ASSIGNMENT , 2 = REFERENCE, 1 = SUBSCRIPT USE, 0 = DEFINITION) 


DCL 

NAME 

TYPE 

517 

AERO DAP 

PROCEDURE 

255 

AERQ_DAP ACT 

BIT(l) 

21 

AERO_DAP_CNT 

INTEGER 

20 

AERO_DAP_PHS 

INTEGER 

516 

AERO_GUID 

PROCEDURE 

256 

AERO_GUID_ACT 

BIT( 1) 

23 

AERO_GUID_CNT 

INTEGER 

22 

AERO GUID PHS 

INTEGER 

519 

FSH_PASS 

INTEGER 

518 

FSH_SEQ 

PROCEDURE 

515 

ORB_NAV 

PROCEDURE 

262 

ORB_NAV_ACT 

BIT(l) 

35 

ORB~NAV~"CNT 

INTEGER 

34 

ORB_NAV PHS 

INTEGER 


ATTRIBUTES & CROSS REFERENCE 

EXTERNAL, VERSIONS XREFs 0 

ALIGNED, INITIAL XREF : 0 0255 2 

SINGLE, ALIGNED, INITIAL XREFs 0 
SINGLE, ALIGNED, INITIAL XREFs 0 
EXTERNAL, VERSIONS XREF: 0 

ALIGNED, INITIAL XREFs 0 0256 2 

SINGLE, ALIGNED, INITIAL XREFs 0 
SINGLE, ALIGNED, INITIAL XREF; 0 
DOUBLE, ALIGNED, STATIC, INITIAL 
6 0526 

XREF: 0 0518 NOT REFERENCED 
EXTERNAL, VERSIONS XREF; 0 

ALIGNED, INITIAL XREF: 0 0262 2 

SINGLE, ALIGNED, INITIAL XREFs 0 
SINGLE, ALIGNED, INITIAL XREFs 0 


0517 2 0525 

0524 

0021 2 0524 

0020 2 0524 

0516 2 0523 

0522 

0023 2 0522 

0022 2 0522 

XREFs 0 0519 2 0520 2 0522 


0515 2 0521 

0520 

0035 2 0520 

0034 2 0520 


2 0524 


OF POOR QUALfTY 



HAL/S STD 360-24.20 


INTERMETRICS 


INC 


APRIL 27, 1987 


11:19:28. 56 


C.TJ 

o 


STMT 


SOURCE 


CURRENT SCOPE 


523 Ml ORBJJAV; 


0RB_NAV 


523 Ml PROCEDURE ) 


0RB_NAV 


r 1 




cl 

FUNCTION: MAINTAIN ESTIMATE OF VEHICLE STATE VECTOR AND COMPUTE 

c| 

STATE VECTOR 

DERIVED PARAMETERS ' 

cl 

INPUTS: T ATTITUDE 

- 

TIME TAG OF STATE VECTOR 

cl 

RI 

- 

POSITION VECTOR 

cl 

VI 

- 

VELOCITY VECTOR 

cl 

AI 

- 

ACCELERATION VECTOR 

cl 

QIB 

- 

ATTITUDE QUATERNION 

cl 

PHI 

- 

BANK ANGLE 

Cl 

OUTPUTS: T_NAV 

- 

TIME TAG OF STATE VECTOR 

cl 

R_NAV 

- 

POSITION VECTOR 

Cl 

V NAV 

- 

VELOCITY VECTOR 

cl 

A_NAV 

- 

ACCELERATION VECTOR 

Cl 

Q B TO I 

- 

ATTITUDE QUATERNION 

Cl 

R NAV MAG 

- 

MAGNITUDE OF POSITION VECTOR 

c! 

UNITER 

- 

UNIT VECTOR IN DIRECTION OF POSITION VECTOR 

cl 

ALT_NAV 

- 

ALTITUDE ABOVE FISHER ELLIPSOID 

Cl 

V_NAV_MAG 

- 

MAGNITUDE OF VELOCITY VECTOR 

cl 

V_REL_NAV 

- 

RELATIVE VELOCITY VECTOR 

cl 

V_REL_MAG 

- 

MAGNITUDE OF RELATIVE VELOCITY VECTOR 

Cl 

RDOT NAV 

- 

RADIAL VELOCITY MAGNITUDE 

cl 

G_LOAD 

- 

SENSED ACCELERATION MAGNITUDE IN G'S 

cl 

ALPHA_NAV 

- 

ANGLE OF ATTACK 

cl 

BETA_NAV 

- 

SIDESLIP ANGLE 

Cl 

PHI_NAV 

- 

BANK ANGLE 

cl 

ENTRY COMPLETE 

- TAEM INTERFACE FLAG 

Cl 

COMMENTS; PERFECT NAVIGATION IS ASSUMED, SO THE STATE VECTOR FROM 

Cl 

cl 

FROM THE ENVIRONMENT MODEL IS COPIED. 


0RB_NAV 
ORB NAV 
ORB NAV 
ORB_NAV 
0RB_NAV 
0RB_NAV 
0RB_MAV 
ORB NAV 
ORB_NAV 
ORB_NAV 
ORB_NAV 
ORB__NAV 
0RB_NAV 
ORB_NAV 
ORB_NAV 
ORB_NAV 
0RB_NAV 
0RB_NAV 
ORB NAV 
0RB_NAV 
ORB_NAV 
0RB_NAV 
ORB_NAV 
ORB_NAV 
ORB_NAV 
ORB_NAV 
ORB_NAV 
ORB_NAV 
ORB_NAV 


Cl 

Cl LOCAL VARIABLES 
Cl 


| ORB NAV 
| ORB_NAV 
I ORB_NAV 


DECLARE VREL_BODY VECTOR(3) SINGLE INITIAL! OH 


524 Ml 


I ORB„NAV 



AVN flyo I *IHd = AVN IHd |H Of 5 

AVN~aao t j 3 

AVNOdO | 339NV MNV8 0V3d |3 

AVN 8«0 | 1 3 

AVN~«aO I = l“0l“9"b IN 625 

I 13 

AVN~flHO 1 1 3 

AVN_8do I NoiNuaxvnb aorumv avsa |o 

AVN 8H0 l 1 3 

AVN oao I <IV = AVN - V |W 825 

I 13 

AVN flao | *IA = AVN - A JW Z2S 

I 13 

AVN“«ao | <ia = AVN~a |h 92s 

I • 13 

AVN”flaO | 1 3 

AVN flaO | aO_L33A 31V1S 0V3d |3 

AVN 8d0 | 1 3 

avn - oao 1 b «3aniiu.v“i = AVN - nwi"i in 525 

AVN~8aO | - 1 3 

AVN - aaO | 9V1 3HI1 OV3a 1 3 

AVN sao | 1 3 

avn - oao | 1 3 

AVN'aaO I , aoi33A 3J.V1S lN3NNOaiAN3 AdOO 13 

avn aao | 1 3 

3d03S iNaaano aoanos mis 

95*82 *61 =11 Z86I « LZ UddV “3NI ‘ S 3 I U 3 H H H N I 02*t>2-09£ (US S/1VH 
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11:19:28.56 

STMT 

SOURCE 


CURRENT SCOPE 


c 

c 

c 

531 H 
E 

532 M 
E 

533 M 
E 

534 M 

E 

534 M 
E 

535 M 
E 

536 M 
E 

537 M 
E 

538 h 
E 

539 M 

C 

C 

C 

E 

540 M 

541 M 
S 

542 M 
S 

C 

C 

c 

543 M 
E 

544 M 


COMPUTE STATE VECTOR DERIVED PARAMETERS 
T_NAV = T_IMU_NAVj 
R_NAV_MAG = ABVALt R_NAV 
UNIT_R = R_NAV / R_NAV_MAG i 

Z 

ALT_NAV = R_NAV„MAG - 11 - EARTH_FLAT 5 EARTH.R / SQRTU ♦ 111 - EARTH_FLAT ) -1)11 

2 

EARTH_POLEJ ))> 

VNAV_MAG a ABVALl V_NAV)j 
V_REL_NAV = V_NAV - ( WE_NAV * RNAV ) \ 

V_REL_MAG = ABVAU V_REL_NAVU 
RDOT_NAV = V_NAV . UNIT_Rj 
G_L0AD = ABVALt A_N A V ) FPS2_TO_G» 

ANGLE OF ATTACK AND SIDESLIP ANGLE 

VREL_BODY = SQFORMt SQPOSEI Q_B_TO_I ), V_REL_NAV)j 

ALPHA_NAV = SARCT AN2 ( VRE L_BOD Y , VREL_BODY I RAD_T0 DEG> 

5 1 

BETA_NAV = ARCSIN(VREL_BODY / V_REL_MAG) RAD_TO_DEG \ 

2 

FLAG SIGNALING TAEH INTERFACE OR SKIP OUT 

IF UALT_NAV > ALT_EXIT ) AND ( RDOT_NAV > 0)) OR ( ALT_NAV < ALT_TAEM ) THEN 
AERQ_BRAKE_COMPLETE = TRUE \ 


i UNIT_R 


ORB_NAV 

ORB_NAV 

0RB_NAV 

ORB_NAV 

ORBJ'IAV 

ORB_NAV 

ORB_NAV 

ORB^NAV 

ORB_NAV 

ORBNAV 

ORB_NAV 

ORB_NAV 

ORB_NAV 

ORB_NAV 

ORB_NAV 

ORB_NAV 

ORB_NAV 

ORB_NAV 

ORB_NAV 

ORB_NAV 

ORB__NAV 

ORBNAV 

ORB_NAV 

ORB_NAV 


HAL/S 

STD 360-24.20 

INTERMETRICS, 

INC, 

APRIL 27, 1987 

11:19:28.56 

STMT 


SOURCE 



CURRENT SCOPE 

545 Ml 

CLOSE ORB_NAV \ 




| ORB_NAV 


**** BLOCK SUMMARY **** 


CTS 

GO 


EXTERNAL FUNCTIONS INVOKED 
SQFORM, SQPOSE, SARCTAN2 

COMPOOL VARIABLES USED 

T_IMU_NAV», ATTITUDE » R_NAV«, RI , V_NAV« , VI, A_NAV*» AI , Q_8_T0_I*> QIB, PHI_NAV*, PHI, T_NAV*» T_IMU_NAV, R_NAV_MAG # » R_NA V 
UNIT_R* , R_NAV_MAG, ALT_NAV», EARTH„FLAT, EARTH„R, UNIT_R, EARTH_POLE , V_NAV_MAG* , V_NAV, V_REL_NAV», HE_NAV, V_REL_MAG* 
V_REL_NAV, RDOT_NAV#, G_LOAD* , A_NAV , FPS2_TO_G, Q_B_TO_I , ALPHA_NAV* , RAD_TO_DEG, BETA_NAV# , V_REL_MAG, ALT_NAV > ALT_EXIT 
RDOT_NAV > ALT JTAEM, AERO_BRAKE_COt1PLETE# 


Sg 


Q 


SJ r- 
O *gr 

c 2“ 

> o 
C rr ‘ 

3 Uj 



cn 


HAL/S STD 260-24.20 INTERMETRICS, INC. APRIL 27, 1987 11 : 19 : 28.56 

**** COMPILATION LAYOUT **** > 

ENV_POOL: EXTERNAL COMPOOLj 
FSH_POOL: EXTERNAL COMPOOLj 
IL.POOL: EXTERNAL COMPOOLj 
DQFORM: EXTERNAL FUNCTION J 
Q_ERR_ANG: EXTERNAL FUNCTION J 
SARCTAN2: EXTERNAL FUNCTION j 
SQFQRM: EXTERNAL FUNCTION j 
SQMULT : EXTERNAL FUNCTION J 
SQPOSE: EXTERNAL FUNCTION J 
SRV_TO_QIL: EXTERNAL FUNCTION J 
ORB_NAV : PROCEDURE j 



HAL/S STD 360-2*. 20 


INTERMETRICS 


INC 


APRIL 27, 1987 


11;I9:28.56 


SYMBOL 4 CROSS REFERENCE TABLE LISTING; 

I CROSS REFERENCE FLAG KEY: * = ASSIGNMENT, 2 = REFERENCE, 1 = SUBSCRIPT USE, 0 = DEFINITION) 

DCL NAME TYPE ATTRIBUTES S CROSS REFERENCE 


88 

A NAV 

3 - VECTOR 

DOUBLE , 

ALIGNED, 

INITIAL XREF: 

0 0088 

* 

0528 

2 

0539 



78 

AERO BRAKE COMPLETE 

BIT! 1 ) 

ALIGNED; 

, INITIAL 

XREF; 0 0078 

* 05** 

NOT REFERENCED 



1 

AI 

3 - VECTOR 

DOUBLE , 

ALIGNED, 

INITIAL XREF: 

0 0001 

2 

0528 





75 

ALPHA NAV 

SCALAR 

SINGLE, 

ALIGNED, 

INITIAL XREF: 

0 0075 

* 

05*1 

NOT REFERENCED 

38* 

ALT EXIT 

SCALAR 

SINGLE , 

ALIGNED, 

INITIAL XREF: 

0 038* 

2 

05*3 





80 

ALT NAV 

SCALAR 

DOUBLE, 

‘ALIGNED, 

INITIAL XREF: 

0 0080 

* 

053* 

2 

05*3 



51* 

ALT TAEM 

SCALAR 

SINGLE » 

ALIGNED, 

INITIAL XREF: 

0 051* 

2 

05*3 





76 

BETA NAV 

SCALAR 

SINGLE, 

ALIGNED, 

INITIAL XREF: 

0 0076 

* 

05*2 

NOT REFERENCED 

518 

CARG 

SCALAR 

SINGLE , 

ALIGNED, 

INPUT-PARM 

XREF: 

0 

0518 





516 

OqFORM 

3 - VECTOR FUNCTION 

DOUBLE » 

INITIAL, 

EXTERNAL, VERSION =2 


XREF: 

0 

0516 






NOT REFERENCED 









*6* 

EARTH FLAT 

SCALAR 

DOUBLE, 

ALIGNED, 

INITIAL XREF: 

0 0*6* 

2 

053* 





*67 

EARTH_POLE 

3 - VECTOR 

DOUBLE, 

ALIGNED, 

INITIAL XREF: 

0 0*67 

2 

053* 





*68 

EARTH R 

SCALAR 

DOUBLE, 

ALIGNED, 

INITIAL XREF: 

0 0*68 

2 

053* 





15 

FPS2 TO G 

SCALAR 

DOUBLE, 

ALIGNED, 

CONSTANT 

XREF: 

0 

0015 

2 

0019 

2 

0539 

83 

G LOAD 

SCALAR 

DOUBLE, 

ALIGNED, 

INITIAL XREF: 

0 0083 

* 

0539 

NOT REFERENCED 

523 

ORB NAV 

PROCEDURE 

XREF : 

0 0523 

NOT REFERENCED 








520 

P 

* - VECTOR 

SINGLE, 

ALIGNED, 

INPUT-PARM 

XREF: 

0 

0520 





1 

PHI 

SCALAR 

SINGLE , 

ALIGNED » 

INITIAL XREF: 

0 0001 

2 

0530 





77 

PHI NAV 

SCALAR 

SINGLE, 

ALIGNED, 

INITIAL XREF: 

0 0077 

* 

0530 

NOT REFERENCED 

521 

Q 

* - VECTOR 

SINGLE , 

ALIGNED, 

INPUT-PARM 

XREF: 

0 

0521 





519 

Q 

* - VECTOR 

SINGLE , 

ALIGNED, 

INPUT-PARM 

XREF: 

0 

0519 





517 

q 

* - VECTOR 

SINGLE , 

ALIGNED, 

INPUT-PARM 

XREF: 

0 

0517 





520 

q 

* - VECTOR 

SINGLE, 

ALIGNED, 

INPUT-PARM 

XREF: 

0 

0520 





516 

q 

* - VECTOR 

DOUBLE , 

ALIGNED, 

INPUT-PARM 

XREF: 

0 

0516 





6* 

Q B TO I 

* - VECTOR 

SINGLE , 

ALIGNED, 

INITIAL XREF; 

0 006* 

* 

0529 

2 

05*0 



517 

q_ERR_ANG 

3 - VECTOR FUNCTION 

SINGLE, 

INITIAL, 

EXTERNAL, VERSI0N=2 


XREF: 

0 

0517 






NOT REFERENCED 









1 

qiB 

* - VECTOR 

DOUBLE , 

ALIGNED, 

INITIAL XREF: 

0 0001 

2 

0529 





522 

R 

3 - VECTOR 

SINGLE, 

ALIGNED, 

INPUT-PARM 

XREF; 

0 

0522 





87 

R_NAV 

3 - VECTOR 

DOUBLE , 

ALIGNED, 

INITIAL XREF; 

0 0087 

* 

0526 

2 

0532 

2 

0533 




2 0536 










89 

R_NAV_MAG 

SCALAR 

DOUBLE , 

ALIGNED, 

INITIAL XREF: 

0 0089 

* 

0532 

2 

0533 

2 

053* 

7 

RAD_TO_DEG 

SCALAR 

DOUBLE, 

ALIGNED, 

CONSTANT 

XREF: 

0 

0007 

2 

05*1 

2 

05*2 

90 

RDOT_NAV 

SCALAR 

DOUBLE, 

ALIGNED, 

INITIAL XREF: 

0 0090 

* 

0538 

2 

05*3 



1 

Rl 

3 - VECTOR 

DOUBLE , 

ALIGNED, 

INITIAL XREF: 

0 0001 

2 

0526 





518 

SARCTAN2 

SCALAR FUNCTION 

SINGLE , 

INITIAL, 

EXTERNAL, VERSI0N=2 


XREF: 

0 

0518 

2 

05*1 

518 

SARG 

SCALAR 

SINGLE , 

ALIGNED, 

INPUT-PARM 

XREF; 

0 

0518 





519 

SQFORM 

3 - VECTOR FUNCTION 

SINGLE, 

IF&TIAL, 

EXTERNAL, VERSIONS 


XREF: 

0 

0519 

2 

05*0 

520 

SqilULT 

* - VECTOR FUNCTION 

SINGLE , 

INITIAL, 

EXTERNAL, VERSI0N=3 


XREF: 

0 

0520 






NOT REFERENCED 









521 

SQPOSE 

* - VECTOR FUNCTION 

SINGLE, 

INITIAL, 

EXTERNAL, VERSIONS 


XREF: 

0 

0521 

2 

o 

o 

522 

SRV_TO_qiL 

* - VECTOR FUNCTION 

SINGLE, 

INITIAL, 

EXTERNAL, VERSION =2 


XREF: 

0 

0522 






NOT REFERENCED 









1 

T_ATTITUDE 

SCALAR 

DOUBLE , 

ALIGNED, 

INITIAL XREF: 

0 0001 

2 

0525 





67 

T IMU_NAV 

SCALAR 

DOUBLE , 

ALIGNED, 

INITIAL XREF: 

0 0067 

* 

0525 

2 

0531 



92 

t_nav 

SCALAR 

DOUBLE , 

ALIGNED, 

INITIAL XREF: 

0 0092 

* 

0531 

NOT REFERENCED 

516 

U 

3 - VECTOR 

DOUBLE, 

ALIGNED s 

INPUT-PARM 

XREF: 

0 

0516 





519 

U 

3 - VECTOR 

SINGLE , 

ALIGNED, 

INPUT-PARM 

XREF: 

0 

0519 





93 

UNIT_R 

3 - VECTOR 

DOUBLE , 

ALIGNED, 

INITIAL XREF: 

0 0093 

* 

0533 

2 

053* 

2 

0538 


ORIGINAL PAGE IS 

OF POOR QUALITY 
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DCL NAME TYPE ATTRIBUTES & CROSS REFERENCE 


522 

V 

3 - VECTOR 

SINGLE, ALIGNED, 

INPUT-PARM 


XREF: 

0 

0522 



94 

VJJAV 

3 - VECTOR 

DOUBLE, ALIGNED, 

INITIAL 

XREF: 

0 

0094 

4 

0527 

2 0535 

2 0536 




2 0538 









95 

V_NAV_MAG 

SCALAR 

DOUBLE , ALIGNED, 

INITIAL 

XREF: 

0 

0095 

4 

0535 

NOT REFERENCED 

96 

V_REL_MAG 

SCALAR 

DOUBLE, ALIGNED, 

INITIAL 

XREF: 

0 

0096 

4 

0537 

2 0542 


97 

V„REL_NAV 

3 - VECTOR 

DOUBLE, ALIGNED, 

INITIAL 

XREF: 

0 

0097 

4 

0536 

2 0537 

2 0540 

1 

VI 

3 - VECTOR 

DOUBLE, ALIGNED, 

INITIAL 

XREF: 

0 

0001 

2 

0527 



524 

VREL_BQDY 

3 - VECTOR 

SINGLE, ALIGNED, 

STATIC, 

INITIAL 

XREF: 

0 

0524 

4 0540 

2 0541 




2 0542 









481 

NE_NAV 

3 - VECTOR 

SINGLE, ALIGNED, 

INITIAL 

XREF: 

0 

0481 

2 

0536 





HAL/S STD 360-24.20 

INTERMETRICS , 

INC. 

APRIL 27, 1987 

14:13:26.85 

STMT 

SOURCE 

* 


CURRENT SCOPE 


515 Ml 
515 Ml 

AERO_GUID; 
PROCEDURE * 


| AEROJ3UID 
I AERO_GUID 

Cl 

DATE: 

4/1/87 

1 AERO GUID 

cl 

FUNCTION: 

NUMERIC PREDICTOR/CORRECTOR ENTRY GUIDANCE ALGORITHM 

1 AERO GUID 

Cl 


FOR THE ENTRY RESEARCH VEHICLE IERV). 

l AERO GUID 

c! 

INPUTS: 

A NAV - SENSED INERTIAL ACCELERATION VECTOR 

| AERO GUID 

cl 


ALPHA NAV - ANGLE OF ATTACK 

1 AERO GUID 

c| 


ALT_NAV - ALTITUDE ABOVE FISHER ELLIPSOID 

I AERO GUID 

cl 


G_LOAD - SENSED ACCELERATION MAGNITUDE IN 6'S 

I AERO GUID 

cl 


R NAV - INERTIAL POSITION VECTOR 

| AERO GUID 

cl 


T_GMT - GREENWICH MEAN TIME 

| AEROJ5UID 

cl 


V_NAV - INERTIAL VELOCITY VECTOR 

j AERO GUID 

cl 


V NAV MAG - INERTIAL VELOCITY MAGNITUDE 

| AEftO_GUID 

cl 


V_REL_MAG - RELATIVE VELOCITY MAGNITUDE 

I AER0J5UID 

cl 


V_RE L_NAV - RELATIVE VELOCITY VECTOR 

j AERO GUID 

cl 

OUTPUTS: 

ALPHA CMD - COMMANDED ANGLE OF ATTACK 

| AERO GUID 

cl 


PHI_CMD - COMMANDED BANK ANGLE 

I AERO_GUID 

Cl 

DESIGNED BY; 

K. SPRATLIN 

I AERO_GUID 

cl 


C.S. DRAPER LABORATORY, INC. 

j AERO GUID 

cl 


MAIL STOP 2B 

| AERO GUID 

Cl 


555 TECHNOLOGY SQUARE 

I AER0_GUID 

Cl 


CAMBRIDGE , MA 02139 

I AERO GUID 

cl 


(617) 25B-2441 

j AERO GUID 

Cl 



1 AERO GUID 


ORIGINAL PAGE IS 
Of POOR QUALITY 



HAL/S 

STD 360-24.20 INTERMETRICS, INC. 

APRIL 27, 1987 

14:13:28.85 

STMT 

SOURCE 


CURRENT SCOPE 

cl 

cl 

Cl 

LOCAL VARIABLES 


I AEROJ3UID 
I AERO GUID 



I AERO_GUID 

516 Ml 

DECLARE ALPHA_DES SCALAR SINGLE 1 


I AERO_GUID 

517 Ml 

DECLARE CL_EST SCALAR SINGLE) 1 


I AERO_GUID 

518 Ml 

DECLARE CQ5PHI_QD0T SCALAR SINGLE INITIAL) -1 . 0 ) ) 


I AERO_GUID 

519 Ml 

DECLARE DELTA_T_PRED SCALAR SINGLE) 


1 AERO_GUID 

520 Ml 

DECLARE E F_F ROM_RE F_ AT_E POC H MATRIX) 3, 3) OOUBLE) 


1 AERO_GUID 

521 Ml 

DECLARE GUID_PASS INTEGER SINGLE) 


1 AERO_GUID 

522 Ml 

DECLARE I_TARGET_EF VECTOR) 3 ) DOUBLE) 


1 AERO_GUID 

523 Ml 

DECLARE INITIALIZE_GUIDANCE BOOLEAN INITIAL) TRUE ) ) 


I AERO_GUID 

524 Ml 

DECLARE PHI_DES SCALAR SINGLE) 


I AERO_GUID 

525 Mj 

cl 

cl 

cl 

DECLARE RHQ_NAV SCALAR SINGLE) 


I AERO_GUID 

ATMOSPHERIC PROPERTIES STRUCTURE 


I AERO GUID 
I AERO_GUID 



I AEROJ3UID 

526 Ml 

STRUCTURE ATMOSPROP: 


1 AERO_GUID 

526 Ml 

1 H SCALAR SINGLE, 


I AEROJ3UID 

526 Ml 

1 RHO SCALAR SINGLE, 


1 AERO_GUID 

526 Ml 

1 TS SCALAR SINGLE, 


I AEROJ5UID 

526 Ml 

l TM SCALAR SINGLE) 


i AERO_GUID 
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HAL/S STD 360-24.20 INTERMETRICS, INC. 

STMT SOURCE 

C| 

Cl COMPUTE BANK MAGNITUDE FOR HEAT RATE CONTROL 
Cl 

540 Ml 1 CALL HEAT_RATE_CONTROL ( 

C| 

C| UPDATE COMMANDED ATTITUDE 

C| 

541 Ml 1 CALL ATTITUDE_COMMANDj 

542 Ml END i 


APRIL 27, X987 14:13s 28.85 

CURRENT SCOPE 

| AEROJ3UID 
1 AERO^GUID 
I AEROJ3UID 

| AERO_GUID 

I AEROJ3UID 
I AERO_GUID 
I AER0_GUID 

I AERO_GUID 

I AERO_GUID 
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HAL/S STD 360-24.20 

INTERMETRICS, 

INC. 

APRIL 27, 1987 

14:13:28.85 

STMT 

SOURCE 



CURRENT SCOPE 


COMPOOL VARIABLES USED 

EF_TO_RE F_AT_E POCH , ALPHA_EI, PHI_EI, ALPHA_CMD*» ALPHA_MIN» ALPHA_MAX, PHI^CMD*, PHI_MAX» LONG_TARGET, DEG_TO_RAD, LAT_TARGET 
EARTH_FLAT, K_RHO_NAV», K_LOD_NAV* 

OUTER VARIABLES USED 

E F_FROM_REF_AT_EPOCH* » GUID_PASS», ALPHA_DES* , PHI_OES*» I_TARGET„EF« 


cr> 

co 



HAL/S STD 360-24.20 INTERMETRICS t INC. 

STMT SOURCE 

571 Ml FILTERS: 

571 Ml PROCEDURE ) 

ci : 

C| FUNCTION: IN-FLIGHT ACCELERATION MEASUREMENT FILTERING 
C I 

Cl 

c| LOCAL VARIABLES 

C| 

572 Ml DECLARE A_DRAG_MAG SCALAR SINGLE ) 

573 Ml DECLARE A_LIFT_MAG SCALAR SINGLE i 

574 Ml DECLARE ATMOS A TMOSPROP -STRUCTURE > 

575 Mi DECLARE CD_NOM SCALAR SINGLE j 

576 Ml DECLARE CL_NOM SCALAR SINGLE J 

577 Ml DECLARE LOD_MEAS SCALAR SINGLE » 

578 Ml DECLARE LOD_NOM SCALAR SINGLE \ 

579 Ml DECLARE MACH SCALAR SINGLE | 

580 Ml DECLARE RHO_MEAS SCALAR SINGLE ) 

5S1 Ml DECLARE V_J3AR SCALAR SINGLE \ 


APRIL 27, 1987 14:13:28.85 

CURRENT SCOPE 

S FILTERS 

I FILTERS 

I FILTERS 
I FILTERS 
I FILTERS 

I FILTERS 
I FILTERS 
I FILTERS 

I FILTERS 

I FILTERS 

I FILTERS 

I FILTERS 

I FILTERS 

I FILTERS 

I FILTERS 

I FILTERS 

I FILTERS 

I FILTERS 


C| 

C| LOOK UP OF NOMINAL DENSITY AND L/D 
C| 

El - 

582 Ml CALL USATMOS62( R_NAV, EARTH_POLE ) ASSIGN! ATMOS ) » 

El + 

583 Ml CALL AERO_PARAMETERS(V_REL_MAG, ATMOS) ASSIGN! V_BAR, MACH H 

584 Ml CALL LOOKUP! ALP HA_NAV, ATMOS. H„ V_BAR, MACH) ASSIGN! CL_NOM, CD_NOM)j 

585 Ml LOD_NOM = CL_NOM / CD_NOM) 



E| _ 

586 Ml A_DRAG_MAG = -!A_NAV . V_REL_NAV) / V_RE L_MAG J 


I FILTERS 
I FILTERS 
I FILTERS 

I 

I FILTERS 

I 

l FILTERS 

I FILTERS 

I FILTERS 

I FILTERS 
I FILTERS 
I FILTERS 

1 

I FILTERS 



CD 

cn 


HAL/S STD 360-24.20 
STMT 


INTERMETRICS, INC. 

SOURCE 


APRIL 27, 1987 14:13:28.85 

CURRENT SCOPE 


El 

587 Ml 

C| 

C| 

cl 

588 H| 
El 

589 Ml 

Cl 

C| 

Cl 

590 Ml 

591 Ml 

C| 

Cl 

C| 

592 Ml 

C| 

C| 

Cl 

593 Ml 

594 Ml 


A_LIFT_MAG = SQRTi A_NAV . AJWV - A_DRAG_MAG A_DRAG_MAG ) | 

COMPUTE MEASURED L/D AND DENSITY 

LOD_MEAS = A_LIFT_MAG / A_DRAG_MAG \ • 

2 

RH0_MEAS = 2 A_DRAG_MAG MASS_NAV / <CD_NOM S_REF V_REL_MAG ]) 

FILTER MEASURED L/D AND DENSITY RATIOS 

K_LOD_NAV = Cl - L_OVER_D_FILTER_GAIN I K_LOD_NAV ♦ L_OVER_D„FILTER_GAIN C LOD_MEAS / LOD_NOM)» 
K_RHO_NAV = (1 - K_RHO_FILTER_GAIN ) K_RHO_NAV ♦ K_RHO_FI LTER_GAIN C RHO_ME AS / ATMOS. RHOM 

COMPUTE FILTERED MEASURED DENSITY 

RHOJ4AV = K_RHO_NAV ATMOS. R HO* 

COMPUTE FILTERED ESTIMATED CL 

CL_EST = K_LOD_NAV CL_NOM* 

CLOSE FILTERS* 


I 

I FILTERS 

I FILTERS 
I FILTERS 
I FILTERS 

I FILTERS 


I FILTERS 

I FILTERS 
I FILTERS 
I FILTERS 

I FILTERS 

I FILTERS 

I FILTERS 
I FILTERS 
I FILTERS 

I FILTERS 

I FILTERS 
I FILTERS 
I FILTERS 

I FILTERS 

I FILTERS 


**## BLOCK SUMMARY **** 


OUTER PROCEDURES CALLED 

LOOKUP, AERO_PARAMETERS, USATMOS62 

COMPOOL VARIABLES USED 

R_NAV, EARTH_POLE> V_REL_MAG, ALPHA_NAV , A_NAV , V_RE L_NAV , MASS_NAV» S_REF, K LOD_NAV*, L_OVER_D FILTER GAIN, K LOD NAV 
K_RHO_NAV*, K_RHO_F I LTE REGAIN, K_RH0_NAV “ “ “ 

OUTER VARIABLES USED 
RHO_NAV#, CL_EST* 

OUTER STRUCTURE TEMPLATES USED 
ATMOSPROP 


o 
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tl 
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HAL/S STD 360-24.20 


INTERMETRICS 


INC. 


APRIL 27, 1987 


14:13:28.85 


STMT 
595 M| 

595 Ml 

C| 

Cl 

C| 

C| 

Cl 

Cl 

596 M| 

597 Ml 

598 Ml 

599 Ml 

600 Ml 

601 Ml 

602 Ml 

603 Ml 

604 M| 

605 Ml 

606 Ml 

607 Ml 

608 Ml 

609 Ml 

E| 

610 Ml 
El 

611 Ml 

612 Ml 
El 

613 Ml 

C| 

Cl 


SOURCE 

HEAT_RATE_CONTROL : 

PROCEDURE* 

FUNCTION: CONTROL PEAK HEAT RATE 


LOCAL VARIABLES 

DECLARE COSPHI_l SCALAR SINGLE* 

DECLARE COSPHI_2 SCALAR SINGLE* 

DECLARE FIRST_PASS BOOLEAN INITIAL I TRUE )* 
DECLARE HS_2 SCALAR SINGLE) 

DECLARE K_QDQT SCALAR SINGLE* 

DECLARE K_QDOT_RATE SCALAR SINGLE* 

DECLARE K1_GAIN SCALAR SINGLE* 

DECLARE Kl_QDOT SCALAR SINGLE* 

DECLARE OME G A_Q D OT_SQU A R E D SCALAR SINGLE* 
DECLARE QBAR SCALAR SINGLE* 

DECLARE 3DOT SCALAR SINGLE* 

DECLARE QOOT_PAST SCALAR SINGLE* 

DECLARE QDOT_RATE SCALAR SINGLE* 

DECLARE TWO_ZETA_OMEGA SCALAR SINGLE* 


3.05 

QDOT = 17700.0 SQRT! RHO_NAV ) !V_REL_MAG / 100001 * 

IF ( FIRST_PASS = TRUE) THEN 
DO* 


1 FIRST_PASS = FALSE* 


HEAT RATE CONTROL CONSTANTS 


CURRENT SCOPE 

I HEAT_RATE_CONTROL 

I HEATER A TE_CONTROL 

I HE AT_RATE_COWTROL 
I HEAT_RATE_CONTROL 
I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 
I HEAT_RATE_CONTROL 
I HEATJ?ATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEATJiATE_JDONTROL 

j HEAT_RATE_COWTROL 

I H E AT_R ATE _CONT ROL 

I H E AT_R ATE_C ONT ROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 


I 

I HEAT_RATE_CONTROL 

I 

I HEAT_RATE_CONTROL 
I HEAT_RATE_CONTROL 

I 

I HEAT_RATE_CONTROL 

I HE AT_R ATE_CONTROL 
I HEAT_RATE_CONTROL 



167 


HAL/S 

STMT 

CJ 

614 HI 

615 M| 

616 Ml 

617 Ml 

Cl 

C| 

cl 

618 Ml 

619 Ml 

620 Ml 

620 Ml 

621 Ml 

622 Mj 

623 Ml 


STD 360-24.20 


INTERMETRICS, INC 
SOURCE 


APRIL 27, 1987 


HS_2 = 2 HS) 

K1_GAIN = S_REF SQRT<RHO_SL) 17700.0 / I HS_2 MASS_NAV)» 
OMEGA_QDDT_SQUARED = OMEGAJJDOT OMEGA_QDOT) 
TW0_2ETA_0MEGA = 2 2ETA_QD0T OMEGA_QDOT) 


FIRST PASS INITIALI2ATI0N 

L QDOT_RATE = 0) 

END | 

ELSE 

QDOT_RATE = l QDOT - QDOT_PAST) / DT_AEROGUID> 
QDOT_PAST = QDOT ) 

IF ( QDOT > .5 QDOT_LIMIT ) THEN 
DO) 


C| 

C| 

Cl 

624 Ml 1 


ESTIMATED DYNAMIC PRESSURE 

QBAR = .5 RHO_NAV VRELMAG V_REL_MAG) 


C| 

C| 

Cl 

El 

625 M| 

626 Ml 

627 Ml 


HEAT RATE CONTROL GAINS 


3.05 

K1_QD0T = K1J3AIN CL_EST (V_REL_MAG / 10000.0) EXP( -( ALT_NAV / HS 2))) 

K_QDOT - OMEGA_QDOT_SQUARED / K1_QD0T> 

K_QDOT_RATE = TW0_2E T A_OME G A / K1_QD0T) 


Cl 

Cl 

C| 

628 Ml 

629 Ml 

630 Ml 


HEAT RATE CONTROL EQUATION 

L C0SPHI_1 = K_QDOT_RATE QDOT_RATE / QBAR \ 

L C0SPHI_2 = K_QDOT I QDOT - QDOT_LIMIT ) / QBAR) 

L COSPHI_QDOT = C0SPHI_1 + C0SPHI_2) 


14:13:28.85 

CURRENT SCOPE 

I HEAT_RATE_CONTROL 

I HEAT_RATE__CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 
I HEAT_RATE CONTROL 
I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I H£AT_RATE_CONTROL 

I HEAT_RATE_CONTROL 
I HEAT_RATE_CONTROL 
I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

j HEAT_RATE_CONTROL 
I HE AT_R ATE_CONTRQL 
I HEAT_RATE_CONTROL 

I 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 
I HlfAT_RATE_CONTROL 
I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 

I HEAT_RATE_CONTROL 
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HAL/S 

STD 360-24 , 20 INTERMETRICS, INC. 


APRIL 27, 1987 

14:13:28.85 


STMT 

SOURCE 



CURRENT SCOPE 


631 M| 

END} 



1 HE AT_R ATE_CONTROL 


632 Ml 

ELSE 

* 


8 HEAT_RAT£_CONTROL 


632 Ml 

COSPHI_QDOT = -1.0} 



8 HEAT_RATE_CONTROL 

05 

OO 

633 Ml 

CLOSE HEAT_RATE_CONTROL } 



I H E AT_R ATE_C ONTRO L 


*«»* BLOCK SUMMARY «##* 

COMPOOL VARIABLES USED 

V_RE L_MAG > HS» S_REF> RHO_SL, MASS^NAV, OMEGA_QDOT » ZETA_QDOT, DT_AEROGUID, QDOTJ.IMIT, ALT_NAV 

OUTER VARIABLES USED 

RHO_NAV , CL_EST, COSPHIJ3DOT* 



CD 

CO 


HAL/S STO 360-24.20 INTERMETRICS » INC. 

STMT SOURCE 

634 Ml ATTXTUDE.COMMAND? 

634 Ml PROCEDURE ( 

C| 

C| LOAD ALPHA COMMAND 
C| 

635 Ml ALPHA.CMD = ALPHA.DESi 

Cl - 

C| COMPUTE AND LIMIT PHI COMMAND 

C| 

636 Ml PHI_CMD = PHI.0ES ABSJ V.NAV.MAG - V.FINAL.MAG) / V.MAG.CHANGE \ 

637 Ml PHI.CMD = MIDVALM -PHI.MAX ) , PHI.CMD, PHI.MAX)* 

Cl 

C| ADJUST COMMANDED BANK ANGLE FOR HEAT RATE CONTROL 

C| 

638 M| IF tCOSPHI.QDOT > 0) THEN 

639 Ml DOj 

640 Ml X TEMPORARY COSPHX_CMD SCALAR SINGLE » 

64X Ml 1 COSPHI_CMD = C0S1PHI.CMD DEG.TO.RAD) + COSPHX.QOOTj 

642 Ml 1 COSPHI_CMD a MIDVALJ-1.0, COSPHI.CMD, +1.0 J> 

643 Ml X PHI_CMD = SIGN( PHI.CMD ) ARCCOS* COSPHI_CMD 1 RAD.TO.DEG> 

644 Ml ENDj 

C| 

C| LIMIT ALPHA AND PHI COMMANDS 

C| 

645 Ml ALPHA_CMD = MIDVAL! ALPHA.MIN » ALPHA.CMD, ALPHA.MAX ) i 

646 Ml PHI_CMD = MIDVALU -PHI.MAX ), PHI.CMD, PHI.MAXU 

647 Ml CLOSE ATTITUDE_COMMAND \ 


APRIL 27, 1987 14:13:28.85 

CURRENT SCOPE 

I ATTITUDE.COMMAND 

I ATTI TUDE_COMMAN D 

I ATTITUDE.COMMAND 
I ATTITUDE COMMAND 
I ATTITUDE _COMI1AND 

I ATTI TUD E.COMM AND 

I ATTITUDE.COMMAND 
I ATTI TUD E_C OMMAND 
I ATTITUDE.COMMAND 

I ATTI TUD E_C OMMAND 

I ATTITUDE .COMMAND 

I ATTITUDE COMMAND 
I ATTITUDE.COMMAND 
I ATTITUDE .COMMAND 

I ATTITUDE.COmAND 

I ATTITUDE.COmAND 

I ATTITUDE.COmAND 

I ATTI TUD E_COMMAND 

I ATTI TUD E.C OMMAND 

I ATTI TUD E.C OMMAND 

I ATTITUDE .COMMAND 

I ATTITUDE.COmAND 
I ATTITUDE.COmAND 
I ATTITUDE.COmAND 

I ATTITUDE.COmAND 

I ATTI TUD E.C OMMAND 

I ATTITUDE.COmAND 


**«* BLOCK SUMMARY **** 

COMPOOL VARIABLES USED 

ALPHA.CMD*, PHI.CMD*, V.NAV.MAG, V.FINAL.MAG, V.MAG.CHANGE , PHI.MAX, PHI.CMD, DEG.TO.RAD, RAD.TO.DEG, ALPHA.MIN, ALPHA.CMD 


OUTER VARIABLES USED 

ALPHA.DES, PHI.DES, COSPHI.QDOT 


ORIGINAL PAGE IS 
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HAL/S STD 360-24 . 20 


INTERMETRICS 


INC 


APRIL 27, 1987 


14:13:28.65 


STMT 



SOURCE 

CURRENT SCOPE 

648 

Ml 

CORRECTOR: 


I CORRECTOR 

648 

Ml 

cl 

cl 

cl 

Cl 

cl 
c 1 

PROCEDURE i 


1 CORRECTOR 


FUNCTION; 

PREDICTOR/CORRECTOR SEQUENCER 

I CORRECTOR 
I CORRECTOR 




j CORRECTOR 


LOCAL VARIABLES 

1 CORRECTOR 
1 CORRECTOR 





I CORRECTOR 

649 

Ml 

DECLARE 

ALPHA_TRY SCALAR SINGLE » 

\ CORRECTOR 

650 

Ml 

DECLARE 

DELTA_ALPHA SCALAR SINGLE CONSTANT 3 )) 

| CORRECTOR 

651 

Ml 

DECLARE 

DELTA_PHI SCALAR SINGLE CONSTANT! 3 Is 

I CORRECTOR 

652 

Ml 

DECLARE 

DETERM SCALAR SINGLE INITIAL! 0)| 

1 CORRECTOR 

653 

Ml 

DECLARE 

DDRE_DA SCALAR SINGLE * 

1 CORRECTOR 

654 

Ml 

DECLARE 

DDRE_DP SCALAR SINGLE s 

I CORRECTOR 

655 

Ml 

DECLARE 

DCREJJA SCALAR SINGLE $ 

| CORRECTOR 

656 

Ml 

DECLARE 

DCREJJP SCALAR SINGLE s 

I CORRECTOR 

657 

Ml 

DECLARE 

CRE SCALAR SINGLE 1 

| CORRECTOR 

658 

Ml 

DECLARE 

CRJERROR ARRAY! 3) SCALAR SINGLE » 

| CORRECTOR 

659 

Ml 

DECLARE 

CR_ERR SCALAR SINGLE \ 

I CORRECTOR 

660 

Ml 

DECLARE 

PHI_TRY SCALAR SINGLE s 

I CORRECTOR 

661 

Ml 

DECLARE 

PRED_EXIT BOOLEAN* 

8 CORRECTOR 

662 

Ml 

DECLARE 

DR_ERROR ARRAY! 3 I SCALAR SINGLE $ 

I CORRECTOR 

663 

Ml 

DECLARE 

DR_ERR SCALAR SINGLE * 

| CORRECTOR 

664 

Ml 

DECLARE 

DRE SCALAR SINGLE J 

I CORRECTOR 


665 

Ml 


DO FOR TEMPORARY I 

= 1 TO 3) 

1 CORRECTOR 

666 

Ml 

1 

DO CASE Is 


I CORRECTOR 

667 

Ml 

2 

DOS 


S CORRECTOR CASE 1 

668 

Ml 

3 

ALPHA_TRY 

= ALPHA_DESs 

8 CORRECTOR 

669 

Ml 

3 

PHX_TRY = 

PHI_DESJ 

j CORRECTOR 


CORRECTOR 



HAL/S 

STMT 

670 hi 

671 Ml 

672 Ml 

673 Ml 

674 Ml 

675 Ml 

676 Ml 

677 Ml 

678 Ml 

679 Ml 

C| 

C| 

Cl 

680 Ml 

C| 

ci 

cj 

681 Ml 
S| 

682 Ml 

si 

683 Ml 

C| 

Cj 

Cl 

684 Ml 

si 

685 M| 
S| 

686 Ml 
S| 

687 Ml 

si 


STD 360-24.20 INTERMETRICS » INC. 

SOURCE 

2 END) 

2 DO) 

3 ALPHA_TRY = ALPHA_DES * DELTA_ALPHA) 

3 PHI_TRY = PHIJJES) 

2 END) 

2 DO) 

3 ALPHA_TRY = ALPHA_DES) 

3 PHI_TRY = PHI_DES * DELTA_PHI ) 

2 END) 

1 END) 

CALL PREDICTOR HITH DESIRED CONTROL HISTORY 

1 CALL PREDICTOR) 

STORE FINAL STATE ERRORS 

1 DR_ERROR = DR_ERR) 

I 


APRIL 27, 1987 14:13:28.85 

CURRENT SCOPE 
I CORRECTOR 
| CORRECTOR CASE 2 
I CORRECTOR 
| CORRECTOR 
I CORRECTOR 
I CORRECTOR CASE 3 
I CORRECTOR 
l CORRECTOR 
I CORRECTOR 

I CORRECTOR DO CASE END 

I CORRECTOR 
I CORRECTOR 
I CORRECTOR 

I CORRECTOR 

| CORRECTOR 
| CORRECTOR 
j CORRECTOR 

I CORRECTOR 


1 CR_ERROR = CR_ERR) 
I 

END) 


COMPUTE PARTIALS 


Ui 

QL 

a 

a 

.DA 

= 

(DR_ERROR 

- DR. 

ERROR ) 

/ 




2 


1 


ODRE. 

.DP 

= 

(OR ERROR 

- DR. 

ERROR ) 

/ 




3 


1 


DCREJ5A 


(CR ERROR 

- CR_ERROR i 

/ 




2 


1 


DCRE_ 

.DP 

= 

1 CR_ERROR 

- CR_ 

.ERROR ) 

/ 


3 1 


DELTA_ALPHA) 

DELTA_PHI) 

DELTA_ALPHA) 

DELTA_PHI) 


CORRECTOR 

CORRECTOR 

CORRECTOR 

CORRECTOR 

CORRECTOR 

CORRECTOR 

CORRECTOR 

CORRECTOR 

CORRECTOR 



so 


HAL/S 

STD 360-24.20 

INTERMETRICS , 

STMT 


SOURCE 

688 Ml 

s| 

DRE = DR ERROR } 
1 


689 M| 
S| 

CRE = CR_ERROR \ 
1 



cj 

C| SOLVE SET OF 2 SIMULTANEOUS EQUATIONS 

C | 

690 Ml DETERM = 0CRE_DP DDRE_DA - DDRE_DP DCRE_DA} 

691 Ml IF ! DETERM -= 0) THEN 

692 M| DO} 

693 Ml 1 

694 Ml 1 

695 Ml I 

696 Ml 1 


INC. APRIL 27, 1987 14:13:28.85 

CURRENT SCOPE 

I CORRECTOR 

I 

| CORRECTOR 

I 

I CORRECTOR 
I CORRECTOR 
I CORRECTOR 

| CORRECTOR 

I CORRECTOR 

I CORRECTOR 

I CORRECTOR 

I CORRECTOR 

I CORRECTOR 

I CORRECTOR 


ALPHA_DES = ALPHA_DES + ( DDRE_DP CRE - DCRE_DP DRE) / DETERM} 
ALPHA_DES = MIDVAH ALPHA_MIN» ALPHA_DES» ALPHA_MAX J \ 

PHI_DES = PHI_DES - ( DDRE_DA CRE - DCRE_DA DRE) / DETERM} 
PHI_DES = MIDVAH C-PHI_DES_MAX), PHI_DES, PHI_DES_MAX ) } 


697 Ml 


END} 


I CORRECTOR 



HAL/S STO 360-24. 20 


INTERMETRICS 


INC 


APRIL 27, 1987 


14:13:28.85 


STMT 
698 M 

698 M 

C 

C 

c 

c 

c 

c 

699 M 

C 

C 

C 

700 h 

701 M 

702 M 

* 703 M 

^ 704 M 

CO 

705 M 

706 M 

707 M 

708 M 

709 M 

710 M 

711 M 

712 M 

713 M 

714 M 

715 M 

716 M 

717 M 

718 M 


SOURCE 

PREDICTOR: 

PROCEDURE ) 

FUNCTION: NUMERIC PREDICTOR ALGORITHM 


LOCAL FUNCTION 

DECLARE EARTH_FIXED_FROM_REFERENCE FUNCTION MATRIX! 3, 3) DOUBLE ) 
LOCAL VARIABLES 

DECLARE TOTA L_TI ME_STE PS INTEGER SINGLE ) 

DECLARE A_PRED VECTOR! 3) DOUBLE} 

DECLARE A_PRED_MAG SCALAR SINGLE} 

DECLARE ALPHA_PRED SCALAR SINGLE) 

DECLARE ATMOS ATMOSPROP -STRUCTURE ) 

DECLARE CD_PRED SCALAR SINGLE) 

DECLARE CL_PRED SCALAR SINGLE) 

DECLARE DOT SCALAR DOUBLE) 

DECLARE EF_FROM_REF MATRIX! 3, 3) DOUBLE) 

, DECLARE I_INPLANE VECTOR! 3 ) DOUBLE) 

DECLARE I_NORMAL VECTOR! 3 ) DOUBLE) 

DECLARE INTEG_LOOP SCALAR SINGLE) 

DECLARE IR_E VECTOR! 3 ) DOUBLE) 

DECLARE LOD_PRED SCALAR SINGLE) 

DECLARE MACH_PRED SCALAR SINGLE) 

DECLARE PHI_PRED SCALAR SINGLE) 

DECLARE R_PRED VECTOR! 3 I DOUBLE) 

DECLARE R_MAG_PRED SCALAR DOUBLE) 

DECLARE RDOT_PRED SCALAR SINGLE) 


CURRENT SCOPE 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 


PREDICTOR 



HAL/S STD 360-24 , 20 


INTERMETRICS 


I N C 


APRIL 27, 1987 


14: 13 ;28.85 


■^1 


i>tni SOURCE 

719 Ml DECLARE T_PRED SCALAR DOUBLE » 

720 Ml DECLARE V_BAR_PRED SCALAR SINGLE i 

721 Ml DECLARE V_MAG_PRED SCALAR DOUBLE } 

.722 Ml DECLARE V_PRED VECT0R(3) DOUBLE j 

723 Ml DECLARE V_REL_MAG_PRED SCALAR SINGLE \ 

724 Ml DECLARE V_REL_PRED VECT0R(3) SINGLE j 

725 Ml DECLARE VR_E VECT0R(3) DOUBLE \ 


C| 

C| INITIALIZE PREDICTOR STATE VECTOR 
C| 

726 M| T_PRED = T_GMT > 

El 

727 Ml R_PRED = R_NAVj 

E| 

728 Ml R„MAG_PRED = ABVAL( R_PRED ) \ 

El 

729 Ml V_PRED = V_NAV> 

El 

730 Ml V_MAG_PRED = ABVALf V_PRED ) j 

El - 

731 Ml V_REL_PRED = V_PRED - IWE_NAV * R_PRED)j 
El 

732 M| V_REL_MAG_PREO = ABVAL( V_REL_PRED ) \ 

C| 

cl ANGLE OF ATTACK FOR PREDICTION 
C| 

733 Ml ALPHA„PRED = ALPMA_TRYj 

C| 

cl INITIALIZE 1962 U.S. STANDARD ATMOSPHERE 
C| 

E ! ♦ 

734 Ml CALL USATMOS62! R_PRED f EARTH_POLE ) ASSIGN! ATMOS ) l 

C| 

c| FORCE 1ST TIME STEP TO BE MINIMUM TIME STEP 
C| 


CURRENT SCOPE 
I PREDICTOR 
I PREDICTOR 
I PREDICTOR 
I PREDICTOR 
I PREDICTOR 
I PREDICTOR 
I PREDICTOR 


I PREDICTOR 
I PREDICTOR 
| PREDICTOR 

I PREDICTOR 

I 

I PREDICTOR 

I 

I PREDICTOR 

I 

I PREDICTOR 

I 

I PREDICTOR 

I 

I PREDICTOR 

I 

1 PREDICTOR 

I PREDICTOR 
I PREDICTOR 
I PREDICTOR 

I PREDICTOR 

I PREDICTOR 
I PREDICTOR 
I PREDICTOR 

I 

I PREDICTOR 

I PREDICTOR 
I PREDICTOR 
I PREDICTOR 


HAL/S STD 360-2*+. 20 


INTERNE TRICS 


INC. 


APRIL 27 » 1987 


14:13:28.85 


-vl 

cn 


STMT SOURCE 

El 

735 M| A_PRED = VECTOR « 999999. , 999999., 999999. )j 

S| SSINGLE ,3 

C| 

C| PREDICTOR LOOP 
C| 


CURRENT SCOPE 

I 

I PREDICTOR 

I 

I PREDICTOR 
I PREDICTOR 
I PREDICTOR 


736 Ml DO FOR TEMPORARY TIMESING RE ME NT = 1 TO 5000) , 

Cl 

Cl COMPUTE TIME STEP FOR 4TH ORDER RUNGA-KUTTA INTEGRATION 
Cl 

737 Ml 1 TQTAL_TIME_STEPS = TIME..INCREMENT) 

738 Ml 1 IF < ATMOS. H <= ALT_TAEM_BIAS ) THEN 

739 Ml 1 DELTA_T_PRED = DE LT A_T_PRED_MIN j 

740 Ml 1 ELSE 

740 Ml 1 DO) 

El 

741 M| 2 A_PRED_MAG = ABVAL! A„PRED )> 

742 Ml 2 IF (A_PRED_MAG -s OJ THEN 

743 Ml 2 DELTA_T_PRED = DE LT A_T_PRE D_GAIN / A_PRED_MAG| 

744 Ml 2 ELSE 

744 Ml 2 DELTA_T_PR£D = DELTA_TJ>RED_MAX» 

745 M l 2 DELTA_T_PRED = MIDVAL( DELTA_T_PRED_MIN» DELTA_T_PRED, DELTA_T_PRED_MAX ) ) 

746 Ml 1 END ( 

C| 

C| AERODYNAMIC PROPERTIES LOOK-UP 
Cl - 

El ♦ 

747 Ml 1 CALL AERO_PARAMETERS( V_REL_MAG_PRED , ATMOS i ASSIGN* V_BAR_PRE0 , MACH_PRED )> 

748 Ml 1 CALL LOOKUPS ALPHA_PRED, ATMOS. H, V_BAR_PRED, MACH_PRED ) ASSIGN J CL_PRED , CD_PRED ) ) 

C I 

C| ESTIMATED L/D 

C| r~ 

LOD_PRED = K_LOD_NAV CL_PRED / CD_PRED) 


I PREDICTOR 

I PREDICTOR 
I PREDICTOR 
I PREDICTOR 

| PREDICTOR 

I PREDICTOR 

| PREDICTOR 

i PREDICTOR 

I PREDICTOR 

I 

I PREDICTOR 

S PREDICTOR 

I PREDICTOR 

I PREDICTOR 

I PREDICTOR 

I PREDICTOR 

I PREDICTOR 

I PREDICTOR 
I PREDICTOR 
I PREDICTOR 

I 

I PREDICTOR 

I PREDICTOR 

I PREDICTOR 
I PREDICTOR 
I PREDICTOR 


749 Ml 1 


I PREDICTOR 


ORIGINAL PAGE' 
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CO 


HAL/S STB 360-24.20 INTERMETRICS, 

STMT SOURCE 

C| 

C| EQUATIONS OF MOTION FOR ERV ENTRY 

Cl 

750 Ml 1 DO FOR INTEG_LOOP = 1 TO 4j 

751 Ml 2 TEMPORARY AERO_ACCEL VECTORI 3 I SINGLE* 

752 Ml 2 TEMPORARY CPHI SCALAR SINGLE * 

753 Ml 2 TEMPORARY DRAG_ ACCEL SCALAR SINGLE* 

754 Ml 2 TEMPORARY GRAV_ACCEL VECTORI 3 ) OOUBLE* 

755 Ml 2 TEMPORARY HS_NORM_PRED SCALAR SINGLE I 

756 Ml 2 TEMPORARY I_LAT VECTOR! 3 ) SINGLE* 

757 Ml 2 TEMPORARY I_LIFT VECTORI 3 J SINGLE* 

758 Ml 2 TEMPORARY I_VEL VECTORI 3) SINGLE* 

759 Ml 2 TEMPORARY LIFT_ACCEL SCALAR SINGLE* 

760 Ml 2 TEMPORARY RHO_PRED SCALAR SINGLE* 

761 Ml 2 TEMPORARY SPHI SCALAR SINGLE* 

762 Ml 2 TEMPORARY U_PRED VECTORI 3) DOUBLE* 

763 Ml 2 TEMPORARY Z_PRED SCALAR DOUBLE* 


INC. APRIL 27, 1987 14:13:28.85 

CURRENT SCOPE 

I PREDICTOR 
I PREDICTOR 
I PREDICTOR 

I PREDICTOR 

I PREDICTOR 

I PREDICTOR 

, I PREDICTOR 

I PREDICTOR 
I PREDICTOR 
I PREDICTOR 
I PREDICTOR 
I PREDICTOR 
I PREDICTOR 
I PREDICTOR 
I PREDICTOR 
I PREDICTOR 
I PREDICTOR 


C| 

C| ESTIMATED DENSITY 

C| 

764 Ml 2 RHO_PRED = K_RHO_NAV ATMOS . RHO * 

C| 

C| BANK ANGLE MODEL 
Cl - 

765 Ml 2 PHI_PRED = PHI_TRY ABSI V_MAG_PRED - V_FINAL_MAG ) / V_MAG_CHANGE * 

766 Ml 2 PHI_PRED = MIDVALI ( -PHI_MAX ) , PHI„PRED, PHIJ1AXJ DEG_TO_RAD* 

767 Ml 2 CPHI = COSI PHI_PRE0 J * 

768 Ml 2 SPHI = SIN! PHI_PRED ) * 

C| 

C| AERODYNAMIC ACCELERATION 
C| 


I PREDICTOR 
I PREDICTOR 
1 PREDICTOR 

I PREDICTOR 

I PREDICTOR 
I PREDICTOR 
I PREDICTOR 

I PREDICTOR 

I PREDICTOR 

I PREDICTOR 

I PREDICTOR 

I PREDICTOR 
I PREDICTOR 
I PREDICTOR 


9 



HAL/S STD 360-24.20 


INTERMETRICS 


INC 


APRIL 27, 1987 


14:13:28.85 


*v| 


STMT 


SOURCE 


E 

769 M 

770 Ml 
El 

771 M 
El 

772 Ml 
El 

773 Ml 
El 

774 Ml 

C| 

Cl 

Cl 

El 

775 Ml 
El 

776 Ml 
El 

777 Ml 
El 

777 M| 
El 

778 Ml 

C| 

Cl 

Cl 

El 

779 Ml 

Cl 

Cl 

C| 

780 Ml 

Cl 

C| 

Cl 

El 

781 Ml 


2 

2 DRAG_ ACCEL = .5 RHO_PRED V_REL_MAG_PRED CD_PRED S_REF / MASS_NAVj 

2 LIFT_ACCEL = LOD_PRED DRAG_ACCELj 

2 I_VEL = V_REL_PRED / V_RE L_MAG_PRED ( 

2 IJ.AT = UNIT(I_VEl * R_PRED)» 

2 I_LIFT = UNIT! I_LAT * I_VEL) CPHI ♦ J_LAT SPHI> 

2 AERO_ACCEL = LIFT_ACCEL I_LIFT - DRAG_ACCEL I_VEL> 

GRAVITY ACCELERATION HITH J2 TERM 

2 U_PRED = R_PRED / R_MAG_PREDj 

2 Z_PRE0 = U_PRED . EARTHJ*OLE ) 

2 2 

2 U_PRED = U_PRED + 13 EARTH_J2 / 2) { EARTH_R / R_MAG_PRED } 111 - 5 Z_PRE0 ) U_PRED ♦ 2 

2 Z_PRED EARTH„POLE ) i 

Z 

2 GRAV_ACCEL = -I EARTH_MU / R_MAG_PRED I U_PREDj 

TOTAL ACCELERATION 

2 A_PRED = AERO_ACCEL + GRAV_ACCEL* 

CALL RUNGA-KUTTA INTEGRATOR 
2 CALL INTEGRATOR \ 

STATE PARAMETERS 

2 R„MAG_PRED = ABVAH R_PRED ) i 


CURRENT SCOPE 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 


PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 


PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 

PREDICTOR 


PREDICTOR 


ORIGINAL PAGE IS 
OF POOR QUALITY 



HAL/'S STD 360-g4.20 


INTERHETRICS 


INC 


APRIL 27, 1987 


14:13:28.65 


STMT SOURCE CURRENT SCOPE 

Ej j 

782 Ml 2 V_MAG_PRED = ABVAU V_PRED ) | | PREDICTOR 

C| j PREDICTOR 

C| RELATIVE VELOCITY | PREDICTOR 

j C | | PREDICTOR 

El I 

783 Ml 2 V_REL_PRED = V_PRED - (HE_NAV * R_PRED)J I PREDICTOR 

El - I 

784 M | 2 V_RE L__MAG_PRED = ABVAU V_REL_PRED ) \ I PREDICTOR 

C| j PREDICTOR 

C| 1962 U.S. STANDARD ATMOSPHERE , I PREDICTOR 

C| j PREDICTOR 

El ♦ | 

785 Ml 2 CALL USATMOS62C R_PRED , EARTH_POLE ) ASSIGW ATMOS ) l I PREDICTOR 

786 Ml 1 ' END) | PREDICTOR 

C I . | PREDICTOR 

C| STATE PARAMETER^ | PREDICTOR 

C| | PREDICTOR 

787 Ml 1 T.PRED = T__PRED + DELTA_T_PRED> I PREDICTOR 

E| | 

788 Ml 1 RDOT_PRED = V_PRED . R_PRED / R_MAG_PRED» I PREDICTOR 

C| | PREDICTOR 

C| CHECK FOR ATMOSPHERIC EXIT I PREDICTOR 

C| | PREDICTOR 

789 Ml 1 IF (( ATMOS. H > 400000) AND IRDOTJ>RED > OH THEN I PREDICTOR 

E| | 

790 Ml 1 PRED_EXIT = TRUE \ \ PREDICTOR 

E| | 

791 Ml 1 IF |PRED_EXIT = TRUE) THEN I PREDICTOR 

792 Ml 1 EXIT ) | PREDICTOR 

Cl j PREDICTOR 

C| CHECK FOR TAEM INTERFACE I PREDICTOR 

C| | PREDICTOR 

793 Ml 1 IF ( ATMOS. H <= ALT_TAEM ) THEN I PREDICTOR 

794 Ml 1 EXIT & I PREDICTOR 

795 Ml 


EN0» 


PREDICTOR 



HAL/S STD 360-24.20 


INTERMETRICS 


INC 


APRIL 27, 1987 


14:13:28.85 


STMT SOURCE 

C| 

C| COMPUTE DOWNRANGE AND CROSSRANGE ERRORS 

cl 

El * 

796 Ml EF_FROM_REF = EARTH_FIXED_FRQM_REFERENCEI T_PRED ) l 

El * 

797 Ml IR„E = UNIT ( E F_FR0M_RE F R_PRED)J 

E| * - 

798 Ml VR_E « EFFROMREF !V_PRE0 - HE_NAV * R_PRED)» 

El - - 

799 Ml I_NORMAL = UNIT(IR_E * VR_E)» 


800 Ml I_INPLANE = UNIT! I_TARGET_EF - ( I_TARGET_EF . I_NORMAL ) I_NORMAL)» 

E| - - 

801 Ml DOT = I_INPLANE . I_TARGET_EF» 

802 Ml IF (ABS(DOT) > 1) THEN 

803 Ml DOT = SIGN! DOT)* 


804 Ml CR_ERR = EARTH_R FT_TO_NM ARCCOS(DOT) SIGN! ( I_INPLANE * I_TARGET_EF) . ! I_NORMAL * I_INPLANE ) ) * 
E| 

805 Ml DOT = IR_E . I_INPLANE l 

806 Ml IF ( ABS! DOT 1 > 1J THEN 

807 Ml DOT = SIGN! DOTH 


DRJERR = EARTH_R FT_TO_NM ARCCOS(DOT) SIGN!! IR_E * I_INPLANE ) . I_NORMAL)» 


CURRENT SCOPE 


PREDICTOR 


PREDICTOR 


PREDICTOR 


PREDICTOR 


PREDICTOR 


PREDICTOR 


PREDICTOR 


PREDICTOR 


PREDICTOR 

o © 

PREDICTOR 

Tl ^ 

PREDICTOR 

T3 § 

O <2 


O “2 

PREDICTOR 

» r* 


o 5 

PREDICTOR 

C 2 

PREDICTOR 

> C3 

r- m 

PREDICTOR 




PREDICTOR 



El 

808 Ml 



HAL/S 


STD 360-24.20 


INTERMETRICS 


INC 


A*>RIL 27, 1987 * 14:13:28.85 


STHT SOURCE CURRENT SCOPE 

809 Ml INTEGRATORS j INTEGRATOR 

809 Ml PROCEDURE } j INTEGRATOR 

c| I INTEGRAT0R 

C| FUNCTION: 4TH OROER RUNGA-KUTTA INTEGRATOR ALGORITHM | INTEGRATOR 

d I integrator 

c| I integrator 

C| LOCAL VARIABLES | INTEGRATOR 

C| j INTEGRATOR 

810 Ml DECLARE ACCUM_ ACCEL VECT0R13) DOUBLE } | INTEGRATOR 

811 Ml DECLARE ACCUM_VEL VECTORS > DOUBLE} | INTEGRATOR 

812 Ml DECLARE ORIGJ>OS VECTOR( 3 ) DOUBLE} 1 | INTEGRATOR 

813 Ml DECLARE ORIG_VEL VECTOR(3J DOUBLE} , | INTEGRATOR 

814 Ml DO CASE INTEG.LOOP} J INTEGRATOR 

8 IB Ml 1 DO} | INTEGRATOR CASE 1 

E| | 

83.6 Ml 2 ORIG_POS = R_PRED} | INTEGRATOR 

E| | 

817 M| 2 ORXG_VEL = V_PRED} | INTEGRATOR 

E| | 

818 M| 2 ACCUM_VEL = V_PRED} | INTEGRATOR 

El | 

819 M| 2 ACCUM_ACCEL - A_PRED} | INTEGRATOR 

E| | 

820 Ml 2 R_PRED = ORIGJ>OS + .5 DELTA_T_PRED V_PRED} | INTEGRATOR 

E| - i 

821 M| 2 V_PRED = ORIGJ/EL * .5 DELTA_T_PRED A_PRED} | INTEGRATOR 

822 Ml 1 END} | INTEGRATOR 

823 Ml 1 DO} | INTEGRATOR CASE 2 

E| | 

824 Ml 2 ACCUM_VEL * ACCUM_VEL + 2 V_PRED} | INTEGRATOR 

E| j 

825 Ml 2 ACCUM_ ACCEL = ACCUM_ ACCEL ♦ 2 A_PRED} | INTEGRATOR 



HAL/S STD 360-24.20 INTERMETRICS, INC. APRIL 27, 1987 14:13; 28.85 


STMT 



SOURCE' 

CURRENT SCOPE 

826 

E| 

Ml 

2 

R_PRED = ORIG_POS ♦ .5 DELTA_T_PRED V_PRED) 

1 

I INTEGRATOR 

827 

El 

Ml 

2 

V_PRED = 0RIG_V£L * .5 DELTA_T_PRED A_PRED| 

1 

1 INTEGRATOR 

828 

Ml 

1 

END) 

I INTEGRATOR 

829 

Ml 

1 

DO) 

1 INTEGRATOR CASE 3 

830 

El 

M| 

2 

ACCUM_VE L = ACCUM_VE L ♦ 2 V_PRED) 

1 INTEGRATOR 

831 

El 

Ml 

2 

ACCUM_ ACCEL = ACCUM_ACCEL ♦ 2 A_PRED J 

1 

1 INTEGRATOR 

832 

El 

Ml 

2 

R_PRED = ORIG_POS * DEITA_T_PRED V_PRED) 

1 

1 INTEGRATOR 

833 

El 

Ml 

2 

V_PRED = ORIG_VEL 4- DELTA_T_PRED A_PRED) 

1 

I INTEGRATOR 

834 

Ml 

1 

END) 

1 INTEGRATOR 

835 

Ml 

1 

DO) 

I INTEGRATOR CASE 4 

836 

El 

Ml 

2 

R_PRED = ORXG_POS * CACCUM_VEL + V_PRED i DELTA_T_PRED / 6) 

1 

1 INTEGRATOR 

837 

El 

Ml 

2 

V_PRED = 0RIG_VEL * IACCUM_ACCEL 4 A_PRED) DELTA_T_PRED / 6) 

1 

1 INTEGRATOR 

838 

Ml 

1 

END) 

1 INTEGRATOR 

839 

Ml 


END) 

1 INTEGRATOR DO CASE END 

840 

Ml 

CLOSE INTEGRATOR ) 

1 INTEGRATOR 

**** 

B 1 

0 

CK SUMMARY «*** 



OUTER VARIABLES USED 

INTEG_LOOP > R_PRED, V_PRED» A_PRED, RJ>RED», DELTAJT_PRED, V_PRED* 



~< U) 



oo 

ro 


HAL/S STD 360-24.20 INTERMETRICS , INC. APRIL 27, 1987 14:13:28.85 


STMT 



SOURCE 

CURRENT SCOPE 

841 

Ml 

EARTH_FIXED 

_FROM_REFERENCE : 

1 EARTH_FIXED_FROM_REFER 

841 

MS 

FUNCTION! T 1 

MATRIX! 3, 3 ) DOUBLE ) 

S EARTH_FIXEO_FROM_REFER 

| EARTH FIXED FROM REFER 
1 EARTH FIXED FROM_REFER 


C J 
Cl 

DATE: 

01/01/80 


ci 

FUNCTION: 

COMPUTE (LEFT) REFERENCE INERTIAL TO EARTH-FIXED 

| EARTH FIXED FROM_REFER 


Ci 


TRANSFORMATION MATRIX OR VARIOUS APPROXIMATIONS THEREOF „ 

I EARTH FIXED FROM REFER 


Cl 

INPUTS: 

T - THE NUMBER OF SECONDS ELAPSED SINCE TIME ORIGIN ' 

| EARTH FIXED FROM_REFER 


cl 

OUTPUTS: 

EARTH-FIXED FROM REFERENCE ROTATION MATRIX 

| EARTH FIXED FROM REFER 


Cl 

COMMENTS: 

NONE 

| EARTH FIXED FROM REFER 


Cl 

r 1 

REFERENCE: 

UNDOCUMENTED 

| EARTH FIXED FROM REFER 
I EARTH_FIXED_FROM_REFER 

I EARTH_FIXED_FROM_REFER 

842 

c I 
Ml 

DECLARE 

T SCALAR DOUBLE) 

843 

Ml 

DECLARE 

CLAMBDA SCALAR DOUBLE) 

1 EARTH_FIXE0_FROM_REFER 

844 

Ml 

DECLARE 

SLAMBDA SCALAR DOUBLE) 

I EARTH_FIXED_FROM_REFER 


845 Ml CLAMBDA = COSUT - T_EPOCH ) EARTH_RATE ) * 

846 Ml SLAMBDA = SINMT - T_EPOCH ) EARTH„RATE)) 

847 Ml RETURN (MATRIX (CLAMBDA, SLAMBDA, 0, -SLAMBDA, CLAMBDA , 0, 0, 0, 1) 

S| 3>DOUBLE,3,3 

E| * 

847 Ml EF_FROH_REF_AT_EPOCH)» 


j EARTH_FIXED_FROM_REFER 
t EARTH_FIXED_FROM_REFER 
I EARTH_FIXED FROM_REFER 

I 

I 

l EARTH_FIXED_FROM_REFER 


848 Ml CLOSE EARTH_FIXED_FROM_REFERENCE ) I EARTH_FIXED_FROM_REFER 


«**# BLOCK SUMMARY **#* 

COMPOOL VARIABLES USED 
T^EPOCH, £ARTH_RATE 

OUTER VARIABLES USED 

EF_FROM_REF_AT_EPOCH 



HAL/S STD 360-24.20 INTERMETRICS, INC. APRIL 27, 1987 14:13:28.65 

STMT SOURCE CURRENT SCOPE 

849 M| CLOSE PREOICTOR* | PREDICTOR 


***# BLOCK SUMMARY **** 

OUTER PROCEDURES CALLED 

USATMOS62 » AERO_PARAMETERS, LOOKUP 

COMPOOL VARIABLES USED 

T_GMT , R_NAV, V_NAV t WE_NAV, EARTH_POLE , ALT„TAEM_BIAS, DELTA_T_PREDJ1IN, DELTA_T PRED GAIN, DELTA_T_PRED MAX, K_LOD NAV 
K_RHO_NAV» V_FINAL_MAG , V_MAG_C MANGE , PHI_MAX, OEG_TO_RAD» S_REF, MASSNAV , EARTH_J2, EARTH_R, EARTH_MU, ALT.TAEM, FT_TO_F#1 

OUTER VARIABLES USED 

ALPHAJTRY, DELTA_T_PRED#, DELTA_T_PRED » PHI_TRY , PRED,EXIT*, PRED_EXIT , I_TARGET_EF , CR_ERR#» DR_ERR* 


OUTER STRUCTURE TEMPLATES USED 
ATMOSPROP 



HAL/S STD 360-24, 2G INTERMETRICS, INC. 

STMT SOURCE 

850 Ml CLOSE CORRECTOR f, 

**** BLOCK SUMMARY **** 

COMPOOL VARIABLES USED 

ALPHA_MIN, ALPHA_MAX, PHI_DES_MAX 

OUTER VARIABLES USED 

ALPHA_DES» PHIJJES, ALPHA_DES*» PHI_DES* 


APRIL 27, 1987 14s 13; 28.85 

CURRENT SCOPE 
8 CORRECTOR 



HAL/S STD 360-24.20 


INTERMETRICS 


INC 


APRIL 27 „ 1987 14? 13:28.85 


STMT SOURCE 

851 Ml AERO_PARAMETERS: 

851 Ml PROCEDURE (V_REL_MAG, ATMOS l ASSIGN! V_BAR , MACH )) 

Cl 

C| FUNCTION: COMPUTE AERODYNAMIC FLOW REGIME PARAMETERS 
Cl 

C| 

C| LOCAL VARIABLES 

Cj 

852 Ml DECLARE V_REL_MAG SCALAR SINGLE ) 

853 Ml DECLARE ATMOS ATMOSPROP -STRUCTURE j 

854 Ml DECLARE MACH SCALAR SINGLE l 

855 Ml DECLARE VBAR SCALAR SINGLE \ 

856 Ml DECLARE C„PRIME SCALAR SINGLE > 

857 Ml DECLARE GAMMA_VBAR SCALAR SINGLE » 

858 Ml DECLARE REYNOLDS_NUMBER SINGLE \ 

859 Ml DECLARE $PEED_OF_SOUND SCALAR SINGLE » 

860 Ml DECLARE T_PRIME SCALAR SINGLE \ 

861 Ml DECLARE T_WALL SCALAR SINGLE » 

862 Ml DECLARE VISCOSITY SCALAR SINGLE $ 

C| 

C| LOCAL CONSTANTS 

C| 

863 Ml DECLARE C_BAR SCALAR SINGLE CONSTANT! 25.0 )> 

864 Ml DECLARE DEG_R_TO_DEG_K SCALAR SINGLE CONSTANT! 9 / 5)j 

865 Ml DECLARE GAMMA SCALAR SINGLE CONSTANT! 1.4)} 

El 3 

866 Ml DECLARE U«V_GAS_CONST SCALAR SINGLE CONSTANT! 8.31432 10 1» 

867 Ml DECLARE MOLE_WT_ZERO SCALAR SINGLE CONSTANT! 28.9644 ) \ 

868 Ml DECLARE SPEED_OF_SOUNDlCONST SCALAR SINGLE CONSTANT! SORT! GAMMA UNIV_GAS_CONST / MOLE_WT_ZERO I ) | 


CURRENT SCOPE 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AEROJ>ARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

A ERO_ PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO.PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 

AERO_PARAMETERS 


HAL/S STD 360-24.20 


INTERMETRICS 


INC 


APRIL 27, 1987 


14:13:28.85 


STMT 


SOURCE 

CURRENT SCOPE 


C| 

MACH NUMBER 

I AERO PARAMETERS 


cl 

Cl 

I AERO PARAMETERS 
I AERO„PARAMETERS 

869 

Ml 

SPEED_OF_SOUND = SPEED_0F_S0UND_C0NST M_TO_FT SQRT ( ATMOS . TM ) J 

1 AERO_PARAMETERS 

•870 

Ml 

IF SPEED__OF_SOUND = 0 THEN 

| AERO_PARAMETERS 

871 

Ml 

MACH = 0> 

1 AERO_PARAMETERS 

872 

Ml 

ELSE 

I A E RO_PAR AME TE RS 

872 

Ml 

MACH = V_REL_MAG / SPEED_0F_S0UND> 

1 AERO_PARAMETERS 


C| 

REYNOLDS NUMBER 

I AERO PARAMETERS 


c| 

cl 

1 AERO PARAMETERS 



1 AERO_PARAMETERS 


El 

-6 1.5 

| 

873 

Ml 

VISCOSITY = 1.458 10 KG_TO_SLUG ATMOS. TS / U 110.4 + ATMOS . TS J M_TO_FT>J 

I AERO_PARAMETERS 

874 

Ml 

IF (VISCOSITY = 0 OR ATMOS. H > 300000 J THEN 

| AERO_PARAMETERS 

875 

Ml 

REYNOLDS_NUMBER = 0* 

I AERO_PARAMETERS 

876 

Ml 

ELSE 

| AERO_PARAMETERS 

876 

Ml 

REYNOLOS_NUMBER = ATMOS. RHO V_REL_MAG C_BAR / VISCOSITY j 

1 AERO_PARAMETERS 


Cj 

VISCOUS PARAMETER 

S AERO_PARAMETERS 


Cl 

cl 

1 AERO PARAMETERS 
1 AERO_PARAMET£RS 

877 

Ml 

IF REYNOLDS_NUMBER = 0 THEN 

1 AERO_PARAMETERS 

878 

Ml 

V_BAR = Oi 

1 AERO_PARAMETERS 

879 

Ml 

ELSE 

I AERO_PARAMETERS 



HAL/S STD 360-24.20 


INTERMETRICS 


INC. 


APRIL 27, 1987 


14:13:28.85 


Oo 


STMT 

879 Ml 

C| 

Cl 

C| 

880 Ml 

881 Ml 

882 Ml 

883 Ml 

884 Ml 

885 Ml 

Cl 

Cl 

C| 

886 Ml 

887 Ml 

888 Ml 

889 Ml 

890 Ml 

891 Ml 

892 Ml 

893 Ml 

894 Ml 

895 Ml 

C| 

Cl 

Cl 

El 

896 Ml 
El 

897 Ml 
El 

897 Ml 


SOURCE 

DOj 

T_HALL FOR VISCOUS PARAMETER 

1 IF t ATMOS. H < 249000) THEN 

1 * TJHALL = 2178.0 DEG_R_TO_OEG_ % K | 

1 ELSE IF ( ATMOS. H >= 249000) AND C ATMOS. H < ‘360000) THEN i 

1 T.HALL = (5913.0 - 0.015 ATMOS. H) DEG_R_TO_DEG_K $ 

1 ELSE IF ( ATMOS. H >= 360000) THEN 

1 T_WALL = 504.0 DEG_R_TO_DEG_K* 

GAMMA FOR VISCOUS PARAMETER 

1 IF f ATMOS. H < 100000) THEN 

1 GAMMA_VBAR = 1.4j 

1 ELSE IF ( ATMOS. H >= 100000) AND ( ATMOS. H < 170000) THEN 

1 GAMMA_VBAR a 1.7 - 3.00E-6 ATMOS. Ht 

1 ELSE IF ( ATMOS. H >= 170000) AND C ATMOS. H < 225000) THEN 

1 GAMMA_VBAR = 1.375 - 1.09E-6 ATMOS. H) 

1 ELSE IF ( ATMOS. H >= 225000) AND ( ATMOS. H < 300000) THEN 

1 GAMMA_VBAR = 1.220 - 4.00E-7 ATMOS. H) 

1 ELSE IF ( ATMOS. H >= 300000) THEN 

1 GAMMAVBAR = 1.1* 

T_PRIME AND C_PRIME FOR VISCOUS PARAMETER 


1 T_PRIME = (.468 ♦ .532 T_HALL / ATMOS. TS ♦ .195 (GAMMA_VBAR - 1) MACH / 2) ATMOS. TS» 

.5 -( 5/ATMOS. TS) 

1 C_PRIME = ( T_PRIME / ATMOS. TS) (( ATMOS. TS ♦ 122.1 10 ) / (T_PRIME ♦ 122.1 10 

( 5/T PRIME) 1 
1 )! f 


CURRENT SCOPE 

I AERO_PARAMETERS 

I AERO_PARAMETERS 
I AERO_PARAMETERS 
I AERO^PARAMETERS 

I AERO_PARAMETERS 

I AERO_PARAMETERS 

I AERO_PARAMETERS 

I AERO_PARAMETERS 

I AERO_PARAMETERS 

I AERO_PARAMETERS 

I AERO_PARAMETERS 
| AERO_PARAMETERS 
| AERO_PARAMETERS 

I AERO.PARAMETERS 

I AE RO_P AR AME TE RS 

l AERO_PARAMETERS 

I AE RO_PARAMETERS 

I AERO_PARAMETERS 

I AE RO_P AR AME TE RS 

I AERO_PARAMETERS 

I AERO^PARAMETERS 

I AERO_ PARAMETERS 

I AERO_PARAMETERS 

I AERO_PARAMETERS 
I AERO_PARAMETERS 
1 AERO_PARAMETERS 

I 

I AERO_PARAMETERS 

I 

I AERO_PARAMETE RS 

I 

I AERO_PARAMETERS 


ORIGINAL PAGE IS 
OF POOR QUALITY 



OO 

OO 


HAL/S STD 360-2^.20 INTERMETRICS, INC 

STMT SOURCE 

El .5 

898 Ml 1 V_BAR = MACH IC_PRIME / REYNOLDS_NUMBER I j 

899 Ml END | 

900 Ml CLOSE AERO_PARAMETERS } 

##** BLOCK SUMMARY **** 

COMPOOL VARIABLES USED 
M_TO_FT, KG_TO_SLUG 

OUTER STRUCTURE TEMPLATES USED 
ATMOSPROP 


APRIL 27, 1987 1*;13:28.85 

CURRENT SCOPE 

I 

I A E RO_P AR AME TE RS 
I AERO_PARAMETERS 
| AERQ_PARAMETERS 



HAL/S STD 360-24.20 


INTERMETRICS 


INC. 


APRIL 27, 1987 


14:13:28.85 


OO 

CO 


STMT 


SOURCE 


CURRENT SCOPE 


901 Ml LOOKUPS 

901 Ml PROCEDURE (ALPHA, ALT, V_BAR , MACH) ASSIGN1 CL, CD)) 

Cl 

C| FUNCTION: LOOK-UP OF CL AND CO VERSUS ALPHA, ALTITUDE, VISCOUS 

Cl PARAMETER, AND MACH NUMBER. 

C| INPUTS: ALPHA - ANGLE OF ATTACK (DEG) 

Cl ALTITUDE - ALTITUDE ABOVE FISHER ELLIPSOID (FT) 

Cl V_BAR - VISCOUS PARAMETER 

C| MACH - MACH NUMBER 

C| OUTPUT: CO - DRAG COEFFICIENT 

Cl CL - LIFT COEFFICIENT 

C| FLOWJ?EGIME: 1 = USE ALTITUDE DATA 

Cl 2 = USE V BAR DATA 

Cl 3 = USE MACH DATA 

C| REFERENCE: NOT DOCUMENTED 

C | 

Cl 

cl INPUT AND OUTPUT VARIABLES 
Cl 

902 Ml DECLARE ALPHA SCALAR SINGLE) 

903 Ml DECLARE ALT SCALAR SINGLE) 

904 Ml DECLARE CD SCALAR SINGLE) 

905 Ml DECLARE CL SCALAR SINGLE) 

906 Ml DECLARE MACH SCALAR SINGLE) 

907 Ml DECLARE V_BAR SCALAR SINGLE) 

C| 

C| LOCAL VARIABLES 
C | 

908 Ml DECLARE ALPHA_FRACT SCALAR SINGLE) 

909 Ml DECLARE ALPHA_MAX SCALAR SINGLE CONSTANT 50 ) ) 

910 Ml DECLARE ALPHA_MIN SCALAR SINGLE CONSTANT) 0 ) ) 

911 Ml DECLARE ALT_L SCALAR SINGLE) 

912 Ml OECLARE ALT_MAX SCALAR SINGLE CONSTANT* 537000 > ) 

913 Ml DECLARE ALT_MIN SCALAR SINGLE CONSTANTt 300000 >) 

914 Ml DECLARE ALTJWJN SCALAR SINGLE CONSTANT* ALT_MAX - ALT_MIN ) ) 

DECLARE CD_1 SCALAR SINGLE) 


I LOOKUP 

I LOOKUP 

I LOOKUP 
I LOOKUP 
| LOOKUP 
| LOOKUP 
I LOOKUP 
I LOOKUP 
| LOOKUP 
I LOOKUP 
j LOOKUP 
I LOOKUP 
I LOOKUP 
I LOOKUP 
| LOOKUP 
j LOOKUP 

| LOOKUP 
1 LOOKUP 
I LOOKUP 

I LOOKUP 

I LOOKUP 

I LOOKUP 

| LOOKUP 

I LOOKUP 

| LOOKUP 

j LOOKUP 
I LOOKUP 
I LOOKUP 

| LOOKUP 

I LOOKUP 

I LOOKUP 

I LOOKUP 

I LOOKUP 

I LOOKUP 

I LOOKUP 


915 Ml 


LOOKUP 


ORIGINAL PAGE IS 

OF P00 n O' 'AI.ITY 



dn^ooi | 
dnwoi | 


|H I£6 
|W X£6 


U*r*?lO' *0250* 

<£ZSO- 4 SW0* *8020 ' 4 8*X0* 4 £0X0' ‘0500* I1NVXSN03 319NIS HV1V3S (9)AV«aV 318VA _ aV8~A 3«V133a 

* (0500*0 UNV1SN03 319NIS aVIVOS NIH~MVS - A 3aV1330 |H 0£6 


dnxooi s 
dmiooi s 
dnxooi | 
dnxooi | 

dfDIOOl j 

dnxooi j 
dnxocn j 
dnxooi | 
drmooi j 

dDXOOl I 

dnxooi j 
dnxooi | 
dnxooi | 

dfWOOl j i 

3doos iN3aano 

58*82 : £T 5 M Z86T ‘Z2 IXddV * 3 N I 


*(t»*7Z0*0)IWV!SNO3 319NIS dVIVOS XVW _ «V9~A 3aV133(J |H 626 
<319NIS dVIVOS l _ aVfl~A 3dV133a |H 826 
* 319NIS MV1V3S VHd1V~NIS 3aV1330 |N Z26 
*( 2)1NV1SN03 319NIS «V1V3S NIH _ H3VH 3«V133Q |H 926 

MOT )1NV1SN03 313NIS d¥1V3S XVH~H3VH 3dV133a |H 526 

<319NIS «V1V3S *TH3VH 3aV1330 JH *26 
*319NIS d3931NI C 3aVT33d |H £26 

( 3T3NIS d3931NI X 3aV1330 |H 226 

* 319NXS a393iWI 3HI93a~M0U 3aVT330 |H 126 

«319NIS aV1V3S i3Vad 3HV1330 |W 026 

* 319NIS aV3V3S VHd“IV~S03 3aV1330 |H 616 
1313NIS a\HV3S 2“13 3avi330 |H 8T6 

( 3T3NIS aV1V3S I _ 13 3aVT33G |H ZX6 

f 3T9NIS aVlVOS 2~03 3avi330 |H 9X6 

33anos IW1S 

4 S3IU3HM 3iNI 02**2-09£ (US S/1VH 



HAL/S STO 360-2*. 20 


INTERMETRICS 


INC 


APRIL 27, 1987 


14:13:28.85 


STMT SOURCE CURRENT SCOPE 

Cj | LOOKUP 

Cl ERV AERODYNAMIC DATA TABLES NOT LISTED I LOOKUP 

C| | LOOKUP 



dnxooi 


<QN3 


|M 656 



dnxooi | nna nv / (Niw“nv - i nvi = aovbj , z |w £56 

I r 4 2 i+r «2 r*z _ _ |s 

drwotn | uovaj'vHdiv ( nv~Q3 - nv“a3) ♦ nv ao = z oo 2 |w 256 

I r‘i t+r‘t r* x |s 

dnxooi | UDvaj VHdiv < . nv ao - nv“aai ♦ nv aa = t ao z |w ts6 

I _ r ‘2 i*r 4 2 _ r‘z _ is 

dnxooi j udvmj vHdiv ( nv“io - nv 131 + nv no = z id z |w os& 

I _ r‘t i*r«i _ r«i Is 

dnxooi | uovdd^vMdiv ( nv“i3 - nv 13) + nv~io = i io z |w 666 

dnxoon | <(xvw~nv «nv «Niw~nvnvAaiw - "Tnv z |w 866 

I 3SV3 dflXOOl | *00 I |W Z66 

driMOOl I 10 

dnxooi j viva sannnv jo 

dnxooi | jo 

dnxooi | 13 

drmxn | ssisvi naauoo oniso VHdiv boj 3iviod«3iM jo 

dnxooi i 1 3 

dnXOOl | i 3WI93B~M01d 3SV3 00 |W 966 

dnwxn i *2 = 3WI93H - MOId |M 5+76 

dflXOOl | 3S13 |H 566 

dOXOOl | t£ = 3WI03H - MO1J |M 666 

dflXOOl 1 N3H1 (XVW”H3VW => HOVHI dl 3S13 |W £66 

dnxooi J it = 3WI93B - MOTd |H 266 

dOXOOl | N3H1 (NIW"nV =< 11V > dl |H 166 

dflXOOl | ' <0 = 3WI93B~M01d |M 066 

dnxooi | I o 

dflXOOT | dn XOOT aod 3Sn 01 3WI93a M013 3NIWB3130 1 3 

dnxooi j * 1 3 

dnxooi | i((VHdnV)31V3Nnai - VHdlV) = 13VBd“VHd1V |W 6£6 

dnxooi | Mos ‘it ♦ < VHdiv) 3iv3man ‘tnvAaiw = r |w 9£6 

dnxooi | 1 3 

dnxooi | X3V11V do 319NV dO X30NI 3NIWB3130 |3 

dnxooi | I o 

3doos iwaaano aoanos iw is 

59*92 :si:6t Z86T <Z2 HBdV ' 3 N I ‘S3IB13WB31HI 02“62-09£ 01S S/1VH 


CM 

CD 



ORIGINAL PAGE S3 
OF POOR QUALITY 


I 

drwocn s 
I 

dmiooi j 
I 

dnxooi j 
dn>woi | 
dnwxn | 

£ 3SV3 dnwxn | 
dflMOCH I 

dn>iocn j 

dHMOOl I 

drowxn | 

I 

dnwxn | 
I 

dn>«xn | 

I 

dnwxn | 

I 

dnxooi | 

I 

dnwxn | 
dnwxn | 
dnxooi | 


dnxooi | 

dnxooi | 

dnXOOl | 

dnxooi I 

^ 3SV3 dnxooi | 

dnxooi | 
dflMOOl | 
dnxooi I 

3d03S iNaaano 

S8* 82 5 £!*♦>! Z86T S Z2 HHdV 



CO 

CD 



HAL/S STD 560-24.20 


INTERMETRICS 


INC 


APRIL 27, 1987 


14; 13 : 28.85 


CO 


974 Ml 
S| 

975 Ml 

976 Ml 

977 Ml 

Cl 

C| 

C| 

978 Ml 

979 Ml 

980 Ml 


SOURCE 

2 CD 2 = CD MACH * (CD MACH - CD MACH ) ALPHA TRACT \ 

I+1,J I+1,J+1 I+1,J 

2 TRACT = ( MACH_L / 2 - I)| 

1 END J 

END* 


INTERPOLATE BETWEEN TABLES 

CL « CL_1 ♦ ( CL_2 - CL_I ) FRACTl 
CD - CD_1 + ( CD_2 - CD_1 ) FRACT t 
CLOSE LOOKUP l 


CURRENT SCOPE 
I LOOKUP 

I 

I LOOKUP 
| LOOKUP 

I LOOKUP [£> CASE END 

| LOOKUP 
I LOOKUP 
| LOOKUP 

I LOOKUP 

I LOOKUP 

| LOOKUP 



HAL/S STD 360-24.20 


INTERMETRICS » INC.' 


981 Ml USATMOS62: 

981 Ml PROCEDURE! R, POLE) ASSIGN! ATMOS )l 


Cl 

Cl FUNCTION: COMPUTES THE DENSITY OF THE ATMOSPHERE IO-1SOKMI 
c | AT A SPECIFIED ALTITUDE USING THE 1962 U.S. STANDARD 

Cl ATMOSPHERE MODEL. 

Cl INPUTS: R - INERTIAL POSITION VECTOR I FT ) 

C | POLE - INERTIAL NORTH POLE (NO) 

Cl OUTPUTS: ATMOS. H ~ ALTITUDE (FT) 

C| ATMOS. RHO - DENSITY < LBM/FT**3 ) 

c ATMOS. TS - STATIC TEMPERATURE OF AIR ( DEG K ) 

C| ATMOS. TM - MOLECULAR TEMPERATURE OF AIR (DEG K) 


C| NOMENCLATURE: A 


A - EARTH EQUATORIAL RADIUS (M) 

F - EARTH FLATTENING (ND) 

RO - EARTH RADIUS (PHI = 45 DEGS) (M) 

PHI - GEOGRAPHIC LATITUDE (DEGS) 

PSI - GEOCENTRIC LATITUDE (DEGS) 

GO - SEA-LEVEL GRAVITY (M/S/S) 

MO - MEAN MOLECULAR HEIGHT OF AIR (ND) 

RR - UNIVERSAL GAS CONSTANT 

(JOULES / (DEG K) / (KG-MOL)i 


982 Ml. DECLARE ATMOS ATMOSPROP -STRUCTURE ) 


983 Ml DECLARE R VECTOR! 3 i DOUBLE j 

984 Ml DECLARE POLE VECTOR! 3 ) DOUBLE) 


985 Ml DECLARE GO SCALAR SINGLE CONST ANTI 9.80665 ) ) 

986 Ml DECLARE MO SCALAR SINGLE CONSTANT! 28. 9644 ) ) 


987 Ml DECLARE RR SCALAR SINGLE CONSTANT! 8314. 32 ) I 

988 Ml DECLARE K2 SCALAR SINGLE CONSTANT! GO MO / RR>» 


989 Ml DECLARE M_TO_FT SCALAR DOUBLE CONSTANT! 1 / .3048)) 

990 Ml DECLARE KG_TO_LBM SCALAR SINGLE CONSTANT ( 1 / .45359237)) 


991 Ml DECLARE A SCALAR DOUBLE CONSTANT! 6378178 ) ) 

992 Ml DECLARE F SCALAR DOUBLE CONSTANT! X / 298.32)1 


APRIL 27, 1987 


14:13:28.85 


CURRENT SCOPE 

| USATMOS62 

| USATM0S62 

J USATM0S62 
| USATMOS62 
| USATM0S62 
| USATMOS62 
| US ATMOS 6 2 
| USATMOS62 
| USATM0S6 2 
| USATM0S62 
| USATM0S62 
| USATM0S62 
| USATM0S62 
| USATMOS62 
| USATMOS62 
| USATMOS62 
| US ATMOS 6 2 
| USATMOS62 
| US ATMOS 6 2 
| US ATMOS 6 2 
| US ATMOS 6 2 
| US ATMOS 6 2 
| US ATMOS 6 2 


| USATMOS62 


1 US ATMOS 6 2 
| US ATMOS 6 2 


| USATMOS62 
| US ATMOS 6 2 
| US ATMOS 6 2 
| USATMOS62 


| US ATMOS 6 2 
| US ATMOS 6 2 


| US ATMOS 6 2 
| USATMOS62 



HAL/S STD 360-24.20 


INTERMETRICS, INC. 


APRIL 27, 1987 


14113:28.85 


CO 

CD 


STMT 

El 

993 Ml 

994 Ml 

995 Ml 

996 Ml 

996 M! 

997 Ml 

997 Ml 

998 Ml 

998 Ml 

999 M| 
999 Ml 
999 Ml 

1000 Ml 

1000 Ml 

1001 Ml 
1001 Mi 


SOURCE 

4 2 

DECLARE R0 SCALAR SINGLE CONST ANT ( A SQRTU 1 ♦* !1 - F) ) / II ♦ (1 - F) )))* 

DECLARE I INTEGER SINGLE J 

DECLARE NSEGS INTEGER SINGLE CONST ANT( 12 )> 

DECLARE H_BASE ARRAY! NSEGS ♦ 1) SCALAR SINGLE CONSTANT! 0, 11000, 20000, 32000, 47000, 52000, 
61000, 79000, 90000, 100000, 110000, 120000, 150000 )* 

DECLARE TM_BASE ARRAY! NSEGS + 1) SCALAR SINGLE CONSTANT! 288. 15, 216.65, 216.65, 228.65, 270.65, 
270.65, 252.65, 180.65, 180.65, 210.65, 260.65, 360.65, 960. 65)* 

DECLARE T_BASE ARRAY! NSEGS ♦ 1) SCALAR SINGLE CONSTANT! 288. 15, - 216 . 65, 216.65, 228.65, 270.65, 
270.65, 252.65, 180.65, 180.65, 210.02, 257.00, 349.4^, 892.79)* 

DECLARE RHO_BASE ARRAY! NSEGS ♦ 1) SCALAR SINGLE CONST ANT(1 . 2250 , 0.36392, 0.088035, 0.013225, 
0.0014275, 0.00075943, 0.00025109, 0.00002001, 0.000003170, 0.0000004974, 0.00000009829, 
0.00000002436, 0.000000001836)1 

DECLARE DTMDH ARRAY! NSEGS) SCALAR SINGLE INITIAL! -.0065, .0000, .0010, .0028, .0000, -.0020, - 
.0040, .0000, .0030, .0050, .0100, .0200)* 


CURRENT SCOPE 

I 

i US ATMOS 6 2 
| USATMOS62 
I USATMOS62 
I USATM0S62 
| USATMOS62 
1 US ATMOS 6 2 
| US ATMOS 6 2 
I US ATMOS 6 2 
| US ATM OS 6 2 
I USATMOS62 
| US ATMOS 6 2 
I USATMOS62 
I USATMOS62 
| US ATMOS 6 2 


DECLARE DTDH ARRAY! NSEGS) SCALAR SINGLE INITIAL! 006500, .000000, .001000, .0028000, .000000000 | USATMOS62 
, -.002000, -.004000, .000000, .002937, .0046980, .009249000, .018110)* I USATMOS62 


1002 Ml DECLARE SCALAR DOUBLE, I USATMOS62 

1002 Ml R_MAG , H, SPSI * I USATMOS62 

1003 Ml DECLARE SCALAR SINGLE, ! USATMOS62 

1003 Ml DH, EXPO, ALTADJ, ALTRATIO, GRAVRATIO, TEMP RATIO* I USATMOS62 


C | 

C| DETERMINE ALTITUDE ABOVE FISHER ELLIPSOID 

Cl 

El 

1004 111 R_MAG = ABVAL1 R ) * 

E| 

1005 Ml SPSI - ! R / R_MAG ) . POLE* 


| US ATMOS 6 2 
| USATMOS62 
j USATMQS62 

I 

| US ATMOS 6 2 

I 


US ATMOS 6 2 



29SOH1VSTI | 
29S0WiVSn I 

I 

29sowivso ! 

I 

29SOW1VSO I 

I 

29SOHlVSn g 


29SOW1VSO g 

4 0N3 

I 

|H 6101 

29SOM1VSTI | 

411X3 

2 

!H eiOT 

1 

29SOW±VSn | 

T+N 

N3H1 ( 3SV8 H > HI JI 

2 

IS 

IW ZIOI 

29S0WlVSn | 

4H = I 

2 

I W 9101 

29SOklVSn | 

4S33SN 01 t = N AHVdOdW31 dOJ 00 

I 

|H STOI 

29SOWlVSn | 
1 

4 OIlVUAVaO = OIlVdAVHS 
2 

\ 

|H 6I0T 

O 

29S0WlVSn | 

4H OIlVdAVdO » H 

I 

|W £101 

29SOM1VSO I 

N3H1 (00006 > H> dl 

I 

|H 2tOT 

29sowivsn g 

4ch ♦ oa» / on = oiivHAvao 

I 

| H IIOI 

29S0RLVSH | 

400 


|H 0101 

29SOU1VSH g 

3S13 


|H 0101 

29SOHLVST1 | 

40 = OHd'SOWlV 


|H 6001 

29SOW1VSH | 

29SOW1VSO | 
29SOWIVSn | 
29SORLVSTI | 

N3H1 (OOOOST < HI JI 


[H 8001 

1 “i 

3arUVH3dN31 dVin3310W ONV 4 3dniVd3dH31 3I1V1S * A1ISN3Q 3NIWd3130 

1 J 
13 



1 J 

29S0W1VSTI ! 

413*01 W H = H-S0W1V 


I H ZOOI 

29SOWiVSn i 
1 

4 ( ( ISdS - II (J - 2) J - tUd&S / (1 - If V - I1J-01-M / 9VW~dl = H 
2 


gw 9001 
13 

3d03S lN3ddfl3 

33dfTOS 


IMIS 


S9*8Z:£1tM Z86T 4 LZ UddV ' 3NI ‘S3IM13HH31NI 02*f>2-09£ 01S S/1VH 


400 

N3H1 (00006 > H» JI 


I IW 47201 
I |W £201 


1 1 

i H01Q HO * 3SV8 1 - S1S0W1V 


IS 

I |W 2201 


I I 

4 HO WIG HQ * 3SVQ W1 = MTS0M1V 


IS 

t |H I20T 


4 3SV8 H - H = HO 


IS 

I jH 0201 




HAL/S STD 

360-24.20 INTERMETRICS, INC. 

APRIL 27, 

1987 14:13:28.85 


STMT 

SOURCE 


CURRENT SCOPE 


1025 M| 2 

si 

IF (DTMDH = 0) THEN 
I 


1 U6ATMOS62 
| 


1026 Ml 2 

s| 

ATMOS. RHO = RH0_BASE EXP(-K2 DH / TM BASE )} 
I " I 


I US ATMOS 6 2 


1027 M| 2 

ELSE 


1 USATMOS62 


1027 Ml 2 

DO} 


1 USATMOS62 


1028 M| 3 

s| 

EXPO = 1 + K2 / DTMDH } 
I 


1 USATMOS62 
1 


El 

1029 Ml 3 

s| 

EXPO 

ATMOS. RHO = RH0_BASE 1 TM BASE / ATMOS . TM l ; 

I I 


1 

1 US ATMOS 6 2 
1 


1030 MI 2 

END} 


I USATMOS62 

v 

1031 Ml 1 

END} 


1 US ATMOS 6 2 

CO 

oo 

1032 MI 1 

ELSE 


1 USATMOS62 


1032 Ml 1 

DO} 


1 US ATMOS 6 2 


1033 Ml 2 

s| 

ALTAOJ = RO + H_BASE - TM BASE / DTMDH } 
l ~ 1 I 


1 US ATMOS 6 2 
1 


1034 Ml 2 
s| 

ALTRATIO = (RO + H) / (RO + H_BASE >} 

I 


j US ATMOS 6 2 
| 


1035 Ml 2 

s| 

TEMP RATIO = TM_BASE / ATMOS. TM} 
I 


1 USATMOS62 


El 

1036 Ml 2 
SI 

2 

EXPO = (K2 / DTMDH ) ( RO / ALTADJ) } 
I 


I 

1 US ATMOS 6 2 
J 


1037 Ml 2 
si 

ATMOS. RHO = RHO_BASE TEMPRATIO EXPHK2 / DTMDH ) GRAVRATIO 
I I 

ALTRATIO (DH / 

ALTADJ > S | USATMOS62 
1 


El 

1037 Ml 2 

EXPO 

(TEMPRATIO ALTRATIO ) } 


1 

I US ATMOS 6 2 


1038 M| 1 

END} 


t US ATMOS 6 2 


El 

1039 Ml 1 

3 

ATMOS. RHO = ATMOS. RHO ( KG_TO_LBM / (M_T0_FT G_TO_FPS2))} 


1 

1 US ATMOS 6 2 


1040 Ml 

END} 


! USATM0S62 


1041 Ml CLOSE 

USATMOS62} 


1 USATMOS62 



CO 

CO 


HAL/S STD 360-24. 20 INTERMETRICS, INC. APRIL 27, 1987 14:13:28.85 

STMT SOURCE CURRENT SCOPE 

**** BLOCK SUMMARY **** 

COMPOOL VARIABLES USED 
G_TO_FPS2 

OUTER STRUCTURE TEMPLATES USED 
ATHOSPROP 



200 


HAL/S STD 360-24. 20 INTERMETRICS » INC. 

STMT SOURCE 

1042 Ml CLOSE AERO_GUID> 

***# BLOCK SUMMARY **** 

COHPOOL VARIABLES USED 

G_LOAD, G_RUN_GUI DANCE , ALT_NAV, ALT_FREE2E_GUID „ GUID_PASS_LIM 


APRIL 27 9 1987 14:13:28.85 

CURRENT SCOPE 
I AERO_GUID 



n 

i 

C<J 


ro 

o 


HAL/S STD 360-24.20 INTERMETRICS, INC. APRIL 27, 1967 14:13:28.85 

**** COMPILATION LAYOUT «**« 

FSH.POOL: EXTERNAL COMPOOL? 

IL_POOL: EXTERNAL COMPOOLj 
AERO_GUID: PROCEDURE ? 

INITI AL_GUID : PROCEDURE? 

FILTERS: PROCEDURE? 

HEAT_RATE_CONTROL : PROCEDURE ? 

ATTITUDE_COMMAND : PROCEDURE? 

CORRECTOR: PROCEDURE? 

PREDICTOR: PROCEDURE? 

INTEGRATOR: PROCEDURE? 

EARTH_FIXED_FROM_REFERENCE: FUNCTION? 

AERO_PARAMETERS; PROCEDURE? 

LOOKUP; PROCEDURE? 


USATM0S62: PROCEDURE? 



HAL/S STD 360-24.20 INTERMETRICS » INC. APRIL 27, 1987 14? 13: 28.85 


SYMBOL & CROSS REFERENCE TABLE LISTING: 

(CROSS REFERENCE FLAG KEY: 4 = ASSIGNMENT » 2 = REFERENCE, 1 = SUBSCRIPT USE, 0 = DEFINITION J 

OCL NAME TYPE ATTRIBUTES & CROSS REFERENCE 


526 

ATMOSPROP 

STRUCTURE TEMPLATE 

ALIGNED XREF: 0 0526 2 0574 

2 

0582 

2 

0583 

2 

0584 

2 

0591 




2 0592 2 0704 2 0734 2 0738 

2 

0747 

2 

0748 

2 

0764 

2 

0785 




2 0789 2 0793 2 0853 2 0869 

2 

0873 

2 

0874 

2 

0876 

2 

0880 




2 0882 2 0883 2 0884 2 0866 

2 

0888 

2 

0889 

2 

0890 

2 

0891 




2 0892 2 0893 2 0894 2 0896 

2 

0897 

2 

0982 

2 

1007 

2 

1009 




2 1021 2 1022 2 1026 2 1029 

2 

1035 

2 

1037 

2 

1039 



526 

K 

1 SCALAR 

SINGLE , ALIGNED XREF: 0 0526 

2 

0584 

2 

0738 

2 

0748 

2 

0789 




2 0793 2 0874 2 0880 2 0882 

2 

0883 

2 

0884 

2 

0886 

2 

0888 




2 0889 2 0890 2 0891 2 0892 

2 

0893 

2 

0894 

4 

1007 



526 

RHO 

1 SCALAR 

SINGLE, ALIGNED XREF: 0 0526 

2 

0591 

2 

0592 

2 

0764 

2 

0876 




4 1009 4 1026 4 1029 4 1037 

6 

1039 







526 

TS 

1 SCALAR 

SINGLE, ALIGNED XREF: 0 0526 

2 

0873 

2 

0896 

2 

0897 

4 

1022 

526 

TM 

1 SCALAR 

SINGLE, ALIGNED XREF: 0 0526 

2 

0869 

4 

1021 

2 

1029 

2 

1035 

991 

A 

SCALAR 

DOUBLE, ALIGNED, STATIC, CONSTANT 



XREF: 

0 

0991 

2 

0993 




2 1006 









572 

A_DRAG_MAG 

SCALAR 

SINGLE, ALIGNED, STATIC XREF: 

0 

0572 

4 

0586 

2 

0587 

2 

0588 




2 0589 









573 

A LIFT MAG 

SCALAR 

SINGLE, ALIGNED, STATIC XREF: 

0 

0573 

4 

0587 

2 

0588 



87 

A_NAV 

3 - VECTOR 

DOUBLE, ALIGNED, INITIAL XREF: 

0 

0087 

2 

0586 

2 

0587 



701 

A^PRED 

3 - VECTOR 

DOUBLE, ALIGNED, STATIC XREF: 

0 

0701 

4 

0735 

2 

0741 

4 

0779 




2 0819 2 0821 2 0825 2 0827 

2 

0831 

2 

0833 

2 

0837 



702 

A_PRED_MAG 

SCALAR 

SINGLE, ALIGNED, STATIC XREF: 

0 

0702 

4 

0741 

2 

0742 

2 

0743 

810 

ACCUM_ ACCEL 

3 - VECTOR 

DOUBLE, ALIGNED, STATIC XREF: 

0 

0810 

4 

0819 

6 

0825 

6 

0831 




2 0837 









811 

ACCUM_VEL 

3 - VECTOR 

DOUBLE, ALIGNED, STATIC XREF: 

0 

0811 

4 

0818 

6 

0824 

6 

0830 



* 

2 0836 









751 

AERO ACCEL 

3 - VECTOR 

SINGLE, TEMPORARY XREF: 

0 

0751 

4 

0774 

2 

0779 



515 

AEROJSUID 

PROCEDURE 

XREF: 0 0515 NOT REFERENCED 









851 

AERO_PARAMETERS 

PROCEDURE 

XREF: 2 0583 2 0747 0 0851 









901 

ALPHA 

SCALAR 

SINGLE, ALIGNED, INPUT-PARM 


XREF: 

0 

0901 

0 

0902 

2 

0938 




2 0939 









352 

ALPHA CMD 

SCALAR 

SINGLE, ALIGNED, INITIAL XREF: 

0 

0352 

4 

0553 

4 

0635 

6 

0645 

516 

ALPHA_DES 

SCALAR 

SINGLE, ALIGNED, STATIC XREF: 

0 

0516 

4 

0551 

2 

0635 

2 

0668 




2 0672 2 0676 6 0693 6 0694 









401 

ALPHA_EI 

SCALAR 

SINGLE, ALIGNED, INITIAL XREF: 

0 

0401 

2 

0551 

2 

0553 



908 

ALPHA_FRACT 

SCALAR 

SINGLE, ALIGNED, STATIC XREF: 

0 

0908 

4 

0939 

2 

0949 

2 

0950 




2 0951 2 0952 2 0962 2 0963 

2 

0964 

2 

0965 

2 

0971 

2 

0972 




2 0973 2 0974 









513 

ALPHA MAX 

SCALAR 

SINGLE, ALIGNED, INITIAL XREF: 

0 

0513 

2 

0553 

2 

0695 

2 

0694 

909 

ALPHA_MAX 

SCALAR 

SINGLE, ALIGNED, STATIC, CONSTANT 



XREF: 

0 

0909 






NOT REFERENCED 









513 

ALPHA MIN 

SCALAR 

SINGLE, ALIGNED, INITIAL XREF: 

0 

0513 

2 

0553 

2 

0645 

2 

0694 

910 

ALPHA_MIN 

SCALAR 

SINGLE, ALIGNED, STATIC, CONSTANT 



XREF: 

0 

0910 






NOT REFERENCED 









74 

ALPHA_NAV 

SCALAR 

SINGLE, ALIGNED, INITIAL XREF: 

0 

0074 

2 

0584 





703 

ALPHA PRED 

SCALAR 

SINGLE, ALIGNED, STATIC XREF: 

0 

0703 

4 

0733 

2 

0748 



649 

ALPHA_TRY 

SCALAR 

SINGLE, ALIGNED, STATIC XREF: 

0 

0649 

4 

0666 

4 

0672 

4 

0676 




2 0733 









901 

ALT 

SCALAR 

SINGLE, ALIGNED, INPUT-PARM 


XREF: 

0 

0901 

0 

0903 

2 

0941 


2 0948 



HAL/S STD 360-24.20 


INTERM E. TRICS 


INC 


APRIL 27» 1987 


14:13:28.85 


ro 

o 

CO 


DCL 

NAME 

TYPE 

513 

ALT.FREEZE GUID 

SCALAR 

911 

ALT L 

SCALAR 

912 

ALT_MAX 

SCALAR 

913 

ALTJ1IN 

SCALAR 

79 

ALT NAV 

SCALAR 

914 

ALT RUN 

SCALAR 

513 

ALT TAEM 

SCALAR 

513 

ALT.TAEM BIAS 

SCALAR 

1003 

ALT AD J 

SCALAR 

1003 

ALTRATIO 

SCALAR 

931 

ATMOS 

STRUCTURE 

851 

ATMOS 

STRUCTURE 


574 

ATMOS 

STRUCTURE 

704 

ATMOS 

STRUCTURE 

634 

ATTITUDE_COMMAND 

PROCEDURE 

863 

C BAR 

SCALAR 

856 

C PRIME 

SCALAR 

901 

CD 

SCALAR 

933 

CD..ALT 

SCALAR ARRAY 

937 

CDMACH 

SCALAR ARRAY 

575 

CD_NOM 

SCALAR 

705 

CD PRED 

SCALAR 

935 

CD_VISC 

SCALAR ARRAY 

915 

CD_1 

SCALAR 

916 

CD_2 

SCALAR 

901 

CL 

SCALAR 

932 

CL_ALT 

SCALAR ARRAY 

517 

CL_EST 

SCALAR 

936 

CL_MACH 

SCALAR ARRAY 

576 

CL NOM 

SCALAR 

706 

CL PRED 

SCALAR 

934 

CLJ/ISC 

SCALAR ARRAY 

917 

CL_1 

SCALAR 

918 

CL_2 

SCALAR 


ATTRIBUTES £ CROSS REFERENCE 

SINGLE, ALIGNED, INITIAL XREF : 0 0513 2 0535 

SINGLE, ALIGNED, STATIC XREF: 0 0911 4 0948 2 0953 

SINGLE, ALIGNED, STATIC, CONSTANT XREF: 0 0912 2 0914 

2 0948 

SINGLE, ALIGNED, STATIC, CONSTANT XREF: 0 0913 2 0914 

2 0941 2 0948 2 0953 

DOUBLE, ALIGNE0, INITIAL XREF: 0 0079 2 0535 2 0625 
SINGLE, ALIGNED, STATIC, CONSTANT XREF: 0 0914 2 0953 

SINGLE, ALIGNED, INITIAL XREF: 0 0513 2 0793 

SINGLE, ALIGNED, INITIAL XREF: 0 0513 2 0738 

SINGLE, ALIGNED, STATIC XREF: 0 1003 4 1033 2 1036 2 1037 

SINGLE, ALIGNED, STATIC XREF: 0 1003 4 1034 2 1037 

ATMOSPROP-STRUCTURE, ALIGNED, ASSIGN-PARM XREF: 0 0981 

0 0982 4 1007 4 1009 4 1021 4 1022 4 1026 6 1029 2 1035 

4 1037 6 1039 

ATHOSPROP -STRUCTURE , ALIGNED, INPUT-PARM XREF: 0 0851 0 0853 

2 0869 2 0873 2 0874 2 0876 2 0880 2 0882 2 0883 2 0884 

2 0886 2 0888 2 0889 2 0890 2 0891 2 0892 2 0893 2 0894 

2 0896 2 0897 

ATMOSPROP -STRUCTURE, ALIGNED, STATIC XREF: 0 0574 4 0582 

2 0583 2 0584 2 0591 2 0592 

ATMOS PROP -STRUCTURE, ALIGNEO, STATIC XREF: 0 0704 4 0734 

2 0738 2 0747 2 0748 2 0764 4 0785 2 0789 2 0793 

XREF: 2 0541 0 0634 

SINGLE, ALIGNED, STATIC, CONSTANT XREF: 0 0863 2 0876 

SINGLE, ALIGNED, STATIC XREF: 0 0856 4 0897 2 0898 

SINGLE, ALIGNED, ASSIGN-PARM XREF: 0 0901 0 0904 4 0979 

NOT REFERENCED 

ARRAY! 2,51), SINGLE, ALIGNED, STATIC, CONSTANT XREF: 0 0933 
2 0951 2 0952 

ARRAY! 5,51), SINGLE, ALIGNED, STATIC, CONSTANT XREF: 0 0937 
2 0973 2 0974 

SINGLE, ALIGNEO, STATIC XREF: 0 0575 4 0584 2 0585 2 0589 

SINGLE, ALIGNEO, STATIC XREF: 0 0705 4 0748 2 0749 2 0769 

ARRAY! 8,51), SINGLE, ALIGNED, STATIC, CONSTANT XREF: 0 0935 
2 0964 2 0965 

SINGLE, ALIGNED, STATIC XREF: 0 0915 4 0951 4 0964 4 0973 

2 0?79 

SINGLE, ALIGNED, STATIC XREF: 0 0916 4 0952 4 0965 4 0974 

2 0979 

SINGLE, ALIGNED, ASSIGN-PARM XREF: 0 0901 0 0905 4 0978 

NOT REFERENCED 

ARRAY! 2,51), SINGLE, ALIGNEO, STATIC, CONSTANT XREF: 0 0932 
2 0949 2 0950 

SINGLE, ALIGNED, STATIC XREF: 0 0517 4 0593 2 0625 

ARRAY! 5,51), SINGLE, ALIGNED, STATIC, CONSTANT XREF: 0 0936 
2 0971 2 0972 

SINGLE, ALIGNED, STATIC XREF: O 0576 4 0584 2 0585 2 0593 

SINGLE, ALIGNED, STATIC XREF: 0 0706 4 0748 2 0749 

ARRAY! 8, 51), SINGLE, ALIGNED, STATIC, CONSTANT XREF: 0 0934 
2 0962 2 0963 

SINGLE, ALIGNED, STATIC XREF: 0 0917 4 0949 4 0962 4 0971 

2 0978 

SINGLE, ALIGNEO, STATIC XREF: 0 0918 4 0950 4 0963 4 0972 

2 0978 



204 


HAL/S STD 360-24.20 INTERMETRICS, I 


DCL 

NAME 

TYPE 

843 

CLAMBDA 

SCALAR 

648 

CORRECTOR 

PROCEDURE 

919 

COS ALPHA 

SCALAR 

640 

COSPHI CMD 

SCALAR 

518 

COSPHI_QDQT 

SCALAR 

596 

COSPHI 1 

SCALAR 

597 

COSPHI 2 

SCALAR 

752 

CPHI 

SCALAR 

659 

CR_ERR 

SCALAR 

658 

CR_ERROR 

SCALAR ARRAY 

657 

CRE 

SCALAR 

655 

OCRE DA 

SCALAR 

656 

DCRE DP 

SCALAR 

653 

DDRE DA 

SCALAR 

654 

DDRE DP 

SCALAR 

864 

DEG_R_JO_DEG_K 

SCALAR 

5 

DEG_TO_RAD 

SCALAR 

650 

0ELTA_ ALPHA 

SCALAR 

651 

DELTA_PHI 

SCALAR 

519 

DELTA_T_PRED 

SCALAR 


513 

0ELTA_T PRED GAIN 

SCALAR 


513 

DELTA_T PRED MAX 

SCALAR 


513 

DELTA T PRED MIN 

SCALAR 


652 

DETERM 

SCALAR 


1003 

DH 

SCALAR 


707 

DOT 

SCALAR 


663 

DR_ERR 

SCALAR 


662 

DR_ERROR 

SCALAR 

ARRAY 

753 

DRAG__ACCEL 

SCALAR 


664 

DRE 

SCALAR 


513 

DT AEROGUID 

SCALAR 


1001 

DTDH 

SCALAR 

ARRAY 

1000 

OTMDH 

SCALAR 

ARRAY 

699 

EARTH FIXED FROM REFERENCE 

3X3 MATRIX FUNCTION 

463 

EARTH FLAT 

SCALAR 


464 

EARTH_J2 

SCALAR 


465 

EARTH_MU 

SCALAR 


466 

EARTH_POLE 

3 - VECTOR 

467 

EARTH_R 

SCALAR 


468 

EARTH_RATE 

SCALAR 



C . APRIL 27, 1987 14:13:28.65 

ATTRIBUTES & CROSS REFERENCE 


DOUBLE, 

ALIGNED, 

, STATIC 

XREF: 0 

0843 

4 

0845 

2 

0847 



XREF: 

2 0536 

0 0648 









SINGLE, 

ALIGNED, 

. STATIC 

XREF:. 0 

0919 

NOT USED 




SINGLE, 

TEMPORARY 

XREF: 0 

0640 

4 

0641 

6 

0642 

2 

0643 

SINGLE, 

ALIGNED, 

STATIC, 

INITIAL 

XREF: 

0 

0518 

4 

0630 

4 

0632 

2 0638 

2 0641 










SINGLE , 

ALIGNED, 

STATIC 

XREF: 0 

0596 

4 

0628 

2 

0630 



SINGLE, 

ALIGNED, 

STATIC 

XREF: 0 

0597 

4 

0629 

2 

0630 



SINGLE , 

TEMPORARY 

XREF: 0 

0752 

4 

0767 

2 

0773 



SINGLE, 

ALIGNED, 

STATIC 

XREF: 0 

0659 

2 

0682 

4 

0804 



ARRAY! 3 ) , SINGLE 

:» ALIGNED 

, STATIC 



XREF: 

0 

0658 

4 

0682 

2 0686 

2 0687 

2 0689 









SINGLE , 

ALIGNED, 

STATIC 

XREF: 0 

0657 

4 

0689 

2 

0693 

2 

0695 

SINGLE , 

ALIGNED, 

STATIC 

XREF: 0 

0655 

4 

0686 

2 

0690 

2 

0695 

SINGLE, 

ALIGNED, 

STATIC 

XREF: 0 

0656 

4 

0687 

2 

0690 

2 

0693 

SINGLE, 

ALIGNED, 

STATIC 

XREF: 0 

0653 

4 

0684 

2 

0690 

2 

0695 

SINGLE, 

ALIGNED, 

STATIC 

XREF: 0 

0654 

4 

0685 

2 

0690 

2 

0693 

SINGLE , 

ALIGNED, 

STATIC, CONSTANT 



XREF: 

0 

0864 

2 

0881 

2 0883 

2 0885 










DOUBLE , 

ALIGNED, 

CONSTANT 


XREFs 

0 

0005 

2 

0006 

2 

0007 

2 0555 

2 0556 

2 0641 2 

0766 








SINGLE, 

ALIGNED, 

STATIC, CONSTANT 



XREF: 

0 

0650 

2 

0672 

2 0684 

2 0686 










SINGLE, 

ALIGNED, 

STATIC, CONSTANT 



XREF: 

0 

0651 

2 

0677 

2 0685 

2 0687 










SINGLE , 

ALIGNEO, 

STATIC 

XREF: 0 

0519 

4 

0739 

4 

0743 

4 

0744 

6 0745 

2 0787 

2 0820 2 

0821 2 

0826 

2 

0827 

2 

0832 

2 

0833 

2 0836 

2 0837 










SINGLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0513 

2 

0743 





SINGLE , 

ALIGNEO, 

INITIAL 

XREF: 0 

0513 

2 

0744 

2 

0745 



SINGLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0513 

2 

0739 

2 

0745 



SINGLE, 

ALIGNED, 

STATIC, : 

INITIAL 

XREF: 

0 

0652 

4 

0690 

2 

0691 

2 0693 

2 0695 










SINGLE, 

ALIGNED, 

STATIC 

XREF: 0 

1003 

4 

1020 

2 

1021 

2 

1022 

2 1026 

2 1037 










DOUBLE , 

ALIGNED, 

STATIC 

XREF: 0 

0707 

4 

0801 

2 

0802 

6 

0803 

2 0804 

4 0805 

2 0806 6 

0807 2 

0808 







SINGLE, 

ALIGNED, 

STATIC 

XREF: 0 

0663 

2 

0681 

4 

0808 



ARRAY ( 3 ) , SINGLE 

, ALIGNED 

, STATIC 



XREF: 

0 

0662 

4 

0681 

2 0684 

2 0685 

2 0688 









SINGLE, 

TEMPORARY 

XREF: 0 

0753 

4 

0769 

2 

0770 

2 

0774 

SINGLE, 

ALIGNED, 

STATIC 

XREF: 0 

0664 

4 

0688 

2 

0693 

2 

0695 

SINGLE , 

ALIGNED , 

INITIAL 

XREF: 0 

0513 

2 

0620 





ARRAY! 12 ) , SINGLE, ALIGNED, STATIC, INITIAL 


XREF: 

0 

1001 

2 1022 











ARRAY! 12 ) » SINGLE, ALIGNED, STATIC, INITIAL 


XREF: 

C 

1000 

2 1021 

2 1025 

2 1028 2 

1033 2 

1036 

2 

1037 





DOUBLE 

XREF : 

0 0699 2 

0796 








DOUBLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0463 

2 

0557 





DOUBLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0464 

2 

0777 





DOUBLE , 

ALIGNED, 

INITIAL 

XREF: 0 

0465 

2 

0778 





DOUBLE, 

ALIGNED, 

INITIAL 

XREF: 0 

0466 

2 

0582 

2 

0734 

2 

0776 

2 0777 

2 0785 










DOUBLE , 

ALIGNED, 

INITIAL 

XREF: 0 

0467 

2 

0777 

2 

0804 

2 

0808 

DOUBLE , 

ALIGNED, 

INITIAL 

XREF: 0 

0468 

2 

0845 

2 

0846 





HAL/S STD 360-24.20 


INTERMETRKS » -INC 


APRIL 27, 1987 


14:13:28.85 



DCL 

NAME 

TYPE 

ATTRIBUTES & CROSS REFERENCE 




708 

EF FROM REF 

3X3 MATRIX 

DOUBLE, ALIGNED, STATIC 

XREF: 0 0708 

4 0796 2 0797 

2 0798 


520 

EF FROM REF AT EPOCH 

3X3 MATRIX 

DOUBLE, ALIGNED, STATIC 

XREF: 0.0520 

4 0549 2 0847 



308 

EF TO_REF_AT_EPOCH 

3X3 MATRIX 

DOUBLE, ALIGNED, INITIAL 

XREF: 0 0308 

2 0549 



1003 

EXPO 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 1003 

4 1028 2 1029 

4 1036 





2 1037 





992 

F 

SCALAR 

DOUBLE, ALIGNED, STATIC, CONSTANT 

XREF: 0 0992 

2 0993 





2 1006 





571 

FILTERS 

PROCEDURE 

XREF : 2 0534 0 0571 





598 

FIRST PASS 

BIT(l) 

ALIGNED, STATIC, INITIAL 

XREF: 0 0598 

2 0611 4 0613 



921 

FLOW_REGIME 

INTEGER 

SINGLE, ALIGNED, STATIC 

XREF: 0 0921 

4 0940 4 0942 

4 0944 





4 0945 2 0946 





920 

FRACT 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0920 

4 0953 4 0966 

4 0975 





2 Q978 2 0979 





11 

FT_TG_NM 

SCALAR 

DOUBLE, ALIGNED, CONSTANT 

XREF: 

0 0011 2 0012 

2 0604 





2 0808 





82 

G LOAD 

SCALAR 

DOUBLE, ALIGNED, INITIAL 

XREF: 0 0082 

2 0532 



513 

G RUN GUIDANCE 

SCALAR 

SINGLE, ALIGNED, INITIAL 

XREF: 0 0513 

2 0532 



13 

G_T0_FPS2 

SCALAR 

D0U3LE, ALIGNED, CONSTANT 

XREF: 

0 0013 2 0014 

2 0017 





2 1039 





865 

GAMMA 

SCALAR 

SINGLE, ALIGNED, STATIC, CONSTANT 

XREF: 0 0865 

2 0868 


857 

GAMMA_VBAR 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0857 

4 0887 4 0889 

4 0891 





4 0893 4 0895 2 0896 





754 

GRAV ACCEL 

3 - VECTOR 

DOUBLE, TEMPORARY 

XREF: 0 0754 

4 0778 2 0779 



1003 

GRAVRATIO 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 1003 

4 1011 2 1013 

6 1014 





2 1037 





521 

GUI0_PASS 

INTEGER 

SINGLE, ALIGNED, STATIC 

XREF: 0 0521 

2 0535 6 0537 

2 0538 





4 0539 4 0550 




ro 

513 

GUID_PASS_LIM 

SCALAR 

SINGLE, ALIGNED, INITIAL 

XREF: 0 0513 

2 0538 


O 

985 

GO 

SCALAR 

SINGLE, ALIGNED, STATIC, CONSTANT 

XREF: 0 0985 

2 0988 

Ol 

1002 

H 

SCALAR 

DOUBLE, ALIGNED, STATIC 

XREF: 0 1002 

4 1006 2 1007 

2 1008 





2 1011 2 1012 6 1013 2 

1017 2 1020 

2 1023 2 1034 



526 

H 

SCALAR 

**** SEE STRUCTURE TEMPLATE ATMOSPROP 




996 

H_BASE 

SCALAR ARRAY 

ARRAY! 13 > , SINGLE, ALIGNED, STATIC, CONSTANT XREF: 

0 0996 





2 1017 2 1020 2 1033 2 

1034 




595 

HEAT_RATE_CONTROL 

PROCEDURE 

XREF: 2 0540 0 0595 





513 

HS 

SCALAR 

SINGLE, ALIGNED, INITIAL 

XREF: 0 0513 

2 0614 



755 

HS_NORM_PRED 

SCALAR 

SINGLE, TEMPORARY 

XREF: 0 0755 

NOT USED 



599 

HS_2 

SCALAR 

SINGLE, ALIGNEO, STATIC 

XREF: 0 0599 

4 0614 Z 0615 

2 0625 


665 

I 

INTEGER 

SINGLE, TEMPORARY 

XREF: 4 0665 

2 0666 1 0681 

1 0682 


922 

I 

INTEGER 

SINGLE,' ALIGNED, STATIC 

XREF: 0 0922 

4 0958 1 0962 

1 0963 





1 0964 1 0965 1 0966 4 

0970 1 0971 

1 0972 1 0973 

1 0974 





2 0975 





994 

I 

INTEGER 

SINGLE, ALIGNED, STATIC 

XREF: 0 0994 

4 1016 I 1020 

1 1021 





1 1022 1 1025 1 1026 1 

1028 1 1029 

1 1033 1 1034 

1 1035 





1 1036 1 1037 





709 

X_INPLANE 

3 - VECTOR 

DOUBLE, ALIGNEO, STATIC 

XREF: 0 0709 

4 0800 2 0801 

2 0804 





2 0805 2 0808 





756 

I LAT 

3 - VECTOR 

SINGLE , TEMPORARY 

XREF: 0 0756 

4 0772 2 0773 



757 

I LIFT 

3 - VECTOR 

SINGLE, TEMPORARY 

XREF: 0 0757 

4 0773 2 0774 



710 

I_N0RMAL 

3 - VECTOR 

DOUBLE, ALIGNED, STATIC 

XREF: 0 0710 

4 0799 2 0800 

2 0804 





2 0808 





522 

I_TARGET_EF 

3 - VECTOR , 

DOUBLE, ALIGNED, STATIC 

XREF: 0 0522 

4 0567 2 0800 

2 0601 





2 0804 





758 

I_VEL 

3 - VECTOR 

SINGLE, TEMPORARY 

XREF: 0 0758 

4 0771 2 0772 

2 0773 





2 0774 






HAL/S STO 360-24.2© INTERMETRICS, 


DCL 

NAME 

TYPE 

542 

INITIAL GUID 

PROCEDURE 

523 

INITIALIZE GUIDANCE 

BIT(l) 

711 

IMTEG LOOP 

SCALAR 

809 

INTEGRATOR 

PROCEDURE 

712 

IR_E 

3 - VECTOR 

923 

J 

INTEGER 


957 

K 

INTEGER 

84 

K_LOD_NAV 

SCALAR 

600 

K QDOT 

SCALAR 

601 

K QDOT RATE 

SCALAR 

472 

K RHO FILTER GAIN 

SCALAR 

83 

K_RHO_NAV 

SCALAR 

990 

KG TO LBM 

SCALAR 

18 

KG TO SLUG 

SCALAR 

602 

K1 GAIN 

SCALAR 

603 

K1 QDOT 

SCALAR 

988 

K2 

SCALAR 

475 

L^OVER 0 FILTER GAIN 

SCALAR 

544 

LAT 

SCALAR 

513 

LAT_TARGET 

SCALAR 

759 

LIFT_ACCEL 

SCALAR 

577 

LOD_MEAS 

SCALAR 

578 

LOD_NOM 

SCALAR 

713 

LOD„PRED 

SCALAR 

545 

LONG 

SCALAR 

513 

LONG_TARGET 

SCALAR 

901 

LOOKUP 

PROCEDURE 

10 

h_TO_FT 

SCALAR 

989 

M_TO_FT 

SCALAR 

901 

MACH 

SCALAR 

579 

MACH 

SCALAR 

851 

MACH 

SCALAR 

924 

MACH L 

SCALAR 

925 

MACH_MAX 

SCALAR 

926 

MACH_MIN 

SCALAR 

714 

MACH_PRED 

SCALAR 

314 

MASS NAV 

SCALAR 

867 

MOLE WT ZERO 

SCALAR 

986 

MO 

SCALAR 

1015 

N 

INTEGER 

995 

NSEGS 

INTEGER 

513 

OMEGA^QDOT 

SCALAR 

604 

OMEGA_QDOT_SQUARED 

SCALAR 


N C . APRIL 27, 1987 14: 13s 28.85 

ATTRIBUTES & CROSS REFERENCE 


XREF: 

2 0529 0 0542 









ALIGNED 

, STATIC, INITIAL 

XREF: 0 

0523 

2 

0527 

4 

0530 



SINGLE, 

ALIGNED, STATIC 

XREF: 0 

0711 

4 

0750 

2 

0814 



XREF: 

2 0780 0 0809 









DOUBLE, 

ALIGNED, STATIC 

XREF: 0 

0712 

4 

0797 

2 

0799 

2 

0805 

2 0808 










SINGLE , 

ALIGNED, STATIC 

XREF: 0 

0923 

4 

0938 

1 

0949 

1 

0950 

1 0951 

1 0952 1 0962 

1 0963 1 

0964 

1 

0965 

1 

0971 

1 

0972 

1 0973 

1 0974 









SINGLE, 

TEMPORARY 

XREF: 4 

0957 

2 

0958 

1 

0959 



SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0084 

4 

0569 

6 

0590 

2 

0593 

2 0749 










SINGLE , 

ALIGNED, STATIC 

XREF: 0 

0600 

4 

0626 

2 

0629 



SINGLE, 

ALIGNED, STATIC 

XREF: 0 

0601 

4 

0627 

2 

0628 



SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0472 

2 

0591 





SINGLE , 

ALIGNED, INITIAL 

XREF: 0 

0083 

4 

0568 

6 

0591 

2 

0592 

2 0764 










SINGLE, 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 

0990 

2 

1039 

DOUBLE , 

ALIGNED, CONSTANT 

XREF; 

0 

0018 

2 

0873 



SINGLE , 

ALIGNED, STATIC 

XREF: 0 

0602 

4 

0615 

2 

0625 



SINGLE , 

ALIGNED, STATIC 

XREF; 0 

0603 

4 

0625 

Z 

0626 

2 

0627 

SINGLE , 

ALIGNED, STATIC, 

CONSTANT 



XREF; 

0 

0988 

2 

1026 

2 1028 

2 1036 2 1037 









SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0475 

2 

0590 





SINGLE, 

ALIGNED, STATIC 

XREF: 0 

0544 

4 

0556 

6 

0557 

2 

0558 

SINGLE, 

ALIGNED, INITIAL 

XREF; 0 

0513 

2 

0556 





SINGLE, 

TEMPORARY 

XREF; 0 

0759 

4 

0770 

2 

0774 



SINGLE , 

ALIGNED, STATIC 

XREF: 0 

0577 

4 

0588 

2 

0590 



SINGLE , 

ALIGNED, STATIC 

XREF: 0 

0578 

4 

0585 

2 

0590 



SINGLE , 

ALIGNED, STATIC 

XREF: 0 

0713 

4 

0749 

2 

0770 



SINGLE , 

ALIGNED, STATIC 

XREF: 0 

0545 

4 

0555 

2 

0562 

2 

0566 

SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0513 

2 

0555 





XREF: 

2 0584 2 0748 i 

0 0901 








DOUBLE, 

ALIGNED, CONSTANT 

XREF; 

0 

0010 

2 

0013 

2 

0869 

2 0873 










DOUBLE, 

ALIGNED, STATIC, 

CONSTANT 



XREF; 

0 

0989 

2 

1006 

2 1007 

2 1039 









SINGLE, 

ALIGNED, INPUT-PARM 

XREF; 

0 

0901 

0 

0906 

2 

0943 

2 0969 










SINGLE , 

ALIGNED, STATIC 

XREF: 0 

0579 

4 

0583 

2 

0564 



SINGLE , 

ALIGNED, ASSIGN- 

PARM 

XREF: 

0 

0851 

0 

0854 

4 

0871 

4 0872 

2 0896 2 0898 









SINGLE , 

ALIGNED, STATIC 

XREF: 0 

0924 

4 

0969 

2 

0970 

2 

0975 

SINGLE, 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 

0925 

2 

0943 

2 0969 










SINGLE , 

ALIGNEO, STATIC, 

CONSTANT 



XREF: 

0 

0926 

2 

0969 

SINGLE, 

ALIGNED, STATIC 

XREF: 0 

0714 

4 

0747 

2 

0748 



SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0314 

2 

0589 

2 

0615 

2 

0769 

SINGLE, 

ALIGNED, STATIC, 

CONSTANT 



XREF; 

0 

0867 

2 

0868 

SINGLE, 

ALIGNED, STATIC, 

CONSTANT 



XREF; 

0 

0986 

2 

0988 

SINGLE, 

TEMPORARY 

XREF: 4 

1015 

2 

1016 

1 

1017 



SINGLE, 

ALIGNED, STATIC, 

CONSTANT 



XREF; 

0 

0995 

2 

0996 

2 0997 

2 0998 2 0999 

2 1000 2 

1001 

2 

1015 





SINGLE, 

ALIGNED, INITIAL 

XREF: 0 

0513 

2 

0616 

2 

0617 



SINGLE, 

ALIGNED, STATIC 

XREF: 0 

0604 

4 

0616 

2 

0626 
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OCL 

NAME 

TYPE 

ATTRIBUTES S CROSS REFERENCE 






812 

ORIG_POS 

3 - VECTOR 

DOUBLE, ALIGNED, STATIC 

XREF: 0 0812 

4 

0816 

2 0820 

2 

0826 




2 0832 2 0836 







813 

ORIG_VEL 

3 - VECTOR 

DOUBLE, ALIGNED, STATIC 

XREF: 0 0813 

4 

0817 

2 0821 

2 

0827 




2 0833 2 0837 







25 

PHI_CMD 

SCALAR 

SINGLE, ALIGNED, INITIAL 

XREF: 0 0025 

4 

0554 

4 0636 

6 

0637 




2 0641 6 0643 6 0646 







524 

PHI_DES 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0524 

4 

0552 

2 0636 

2 

0669 




2 0673 2 0677 6 0695 6 

0696 






513 

PHI DES MAX 

SCALAR 

SINGLE, ALIGNED, INITIAL 

XREF: 0 0513 

2 

0696 




402 

PHI El 

SCALAR 

SINGLE, ALIGNED, INITIAL 

XREF: 0 0402 

2 

0552 

2 0554 



513 

PHI_MAX 

SCALAR 

SINGLE, ALIGNED, INITIAL 

XREF: 0 0513 

2 

0554 

2 0637 

2 

0646 




2 0766 







715 

PHI_PRED 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0715 

4 

0765 

6 0766 

2 

0767 




2 0768 







660 

PHI.TRY 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0660 

4 

0669 

4 0673 

4 

0677 




2 0765 







981 

POLE 

3 - VECTOR 

DOUBLE, ALIGNED, XNPUT-PARM XREF 5 

0 

0981 

0 0984 

2 

1005 

661 

PRED EXIT 

BIT!1) 

ALIGNED, STATIC XREF: 0 

0661 4 0790 

2 

0791 




697 

PREDICTOR 

PROCEDURE 

XREF: 2 0680 0 0697 







605 

QBAR 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0605 

4 

0624 

2 0628 

2 

0629 

606 

QDOT 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0606 

4 

0610 

2 0620 

2 

0621 




2 0622 2 0629 







513 

QDOT LIMIT 

SCALAR 

SINGLE, ALIGNED, INITIAL 

XREF: 0 0513 

2 

0622 

2 0629 



607 

QDOT_PAST 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0607 

2 

0620 

4 0621 



608 

QDOT RATE 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0608 

4 

0616 

4 0620 

2 

0628 

981 

R 

3 - VECTOR 

DOUBLE, ALIGNED, INPUT-PARM XREF: 

0 

0981 

0 0983 

2 

1004 




2 1005 







1002 

R„MAG 

SCALAR 

DOUBLE, ALIGNED, STATIC 

XREF: 0 1002 

4 

1004 

2 1005 

2 

1006 

717 

R_MAG_PRED 

SCALAR 

DOUBLE, ALIGNED, STATIC 

XREF: 0 0717 

4 

0728 

2 0775 

2 

0777 




2 0778 4 0781 2 0788 







66 

R NAV 

3 - VECTOR 

DOUBLE, ALIGNED, INITIAL 

XREF: 0 0086 

2 

0582 

2 0727 



716 

R_PRED 

3 - VECTOR 

DOUBLE, ALIGNED, STATIC 

XREF: 0 0716 

4 

0727 

2 0728 

2 

0731 




£ 0734 2 0772 2 0775 2 

0781 2 0783 

2 

0785 

2 0768 

2 

0797 




2 0798 2 0816 4 0820 4 

0826 4 0832 

4 

0836 




6 

RAD TO_DEG 

SCALAR 

DOUBLE, ALIGNED, CONSTANT 

XREF: 

0 

0006 

2 0643 



718 

RDOT_PRED 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF; 0 0718 

4 

0788 

2 0789 



858 

REYNOLDS_NUMBER 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0858 

4 

0875 

4 0876 

2 

0877 




2 0898 







526 

RHO 

SCALAR 

***** SEE STRUCTURE TEMPLATE ATMOSPROP 






999 

RHO_BASE 

SCALAR ARRAY 

ARRAY! 13), SINGLE, ALIGNED, STATIC, CONSTANT 

XREF: 

0 

0999 




2 1026 2 1029 2 1037 







580 

RHO_MEAS 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0580 

4 

0589 

2 0591 



525 

RHO_NAV 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0525 

4 

0592 

2 0610 

2 

0624 

760 

RHO_PRED 

SCALAR 

SINGLE, TEMPORARY 

XREF: 0 0760 

4 

0764 

2 0769 



513 

RHO_SL 

SCALAR 

SINGLE, ALIGNED, INITIAL 

XREF: 0 0513 

2 

0615 




987 

RR 

SCALAR 

SINGLE, ALIGNED, STATIC, 

CONSTANT 


XREF: 

0 0987 

2 

0988 

993 

RO 

SCALAR 

SINGLE, ALIGNED, STATIC, 

CONSTANT 


XREF: 

0 0993 

2 

1011 




2 1033 2 1034 2 1036 







478 

S REF 

SCALAR 

SINGLE, ALIGNED, INITIAL 

XREF: 0 0478 

2 

0589 

2 0615 

2 

0769 

927 

SIN ALPHA 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0927 

NOT USED 



844 

SLAMBDA 

SCALAR 

DOUBLE, ALIGNED, STATIC 

XREF: 0 0844 

4 

0846 

2 0847 



859 

SPEED OF SOUND 

SCALAR 

SINGLE, ALIGNED, STATIC 

XREF: 0 0859 

4 

0869 

2 0870 

2 

0872 

868 

SP E E D_OF_SOUND_CONST 

SCALAR 

SINGLE, ALIGNED, STATIC, 

CONSTANT 


XREF: 

0 0868 

2 

0869 

761 

SPHI 

SCALAR 

SINGLE, TEMPORARY 

XREF: 0 0761 

4 

0768 

2 0773 



1002 

SPSI 

SCALAR 

4 DOUBLE, ALIGNED, STATIC 

XREF: 0 1002 

4 

1005 

2 1006 
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DCL 

NAME 

TYPE 

841 

T 

SCALAR 

998 

T_BASE 

SCALAR ARRAY 

310 

T_EPOCH 

SCALAR 

22 

T GMT 

SCALAR 

719 

T PRED 

SCALAR 

860 

T PRIME 

SCALAR 

861 

T_WALL 

SCALAR 

1003 

TEMPRATIO 

SCALAR 

736 

TIME INCREMENT 

INTEGER 

526 

TM 

SCALAR 

997 

TM_BASE 

SCALAR ARRAY 

700 

TOTAL TIME STEPS 

INTEGER 

526 

TS 

SCALAR 

609 

TWO ZETA OMEGA 

SCALAR 

762 

U_PRED 

3 - VECTOR 

866 

UNIV_GAS_CONST 

SCALAR 

981 

USATMOS62 

PROCEDURE 

901 

V_BAR 

SCALAR 

651 

V_BAR 

SCALAR 

581 

V_BAR 

SCALAR 

928 

V__BAR L 

SCALAR 

929 

V BAR MAX 

SCALAR 

930 

V_BAR_MIN 

SCALAR 

720 

V_BAR_PRED 

SCALAR 

931 

V_BAR_TABLE 

SCALAR ARRAY 

513 

V_FINAL_MAG 

SCALAR 

513 

V__MAG_CHANGE 

SCALAR 

721 

V_MAG_PRED 

SCALAR 

93 

V_NAV 

3 - VECTOR 

94 

V NAV MAG 

SCALAR 

722 

V_PRED 

3 - VECTOR 


851 

V_REL_MAG 

SCALAR 

95 

V_REL_MAG 

SCALAR 

723 

V_RE L_MAG_PRED 

SCALAR 

96 

V_REL_NAV 

3 - VECTOR 

724 

V_REL_PR£D 

3 - VECTOR 

662 

VISCOSITY 

SCALAR 

725 

VR_E 

3 - VECTOR 

480 

WE_NAV 

3 - VECTOR 

546 

X 

SCALAR 


£47 Y 


SCALAR 
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ATTRIBUTES & CROSS REFERENCE 

DOUBLE, ALIGNED, INPUT-PARM XREF: 0 0841 0 0842 2 0845 

2 0846 

ARRAY ( 13 SINGLE, ALIGNED, STATIC, CONSTANT XREF : 0 0998 

2 1022 


DOUBLE , 

ALIGNED, 

INITIAL 

XREF: 

0 

0310 

2 

0845 

2 

0846 



DOUBLE, 

ALIGNED, 

INITIAL 

XREF: 

0 

0022 

2 

0726 





DOUBLE, 

ALIGNED, 

STATIC 

XREF: 

0 

0719 

4 

0726 

6 

0787 

2 

0796 

SINGLE , 

ALIGNED, 

STATIC 

XREF: 

0 

0860 

4 

0896 

2 

0897 



SINGLE, 

ALIGNED, 

STATIC 

XREF: 

0 

0861 

4 

0881 

4 

0883 

4 

0885 

2 0896 
SINGLE, 

ALIGNED, 

STATIC 

XREF: 

0 

1003 

4 

1035 

2 

1037 



SINGLE, 

TEMPORARY 

XREF: 

4 

0736 

2 

0737 






**** SEE STRUCTURE TEMPLATE ATMOSPROP 
ARRAY! 13), SINGLE, ALIGNED, STATIC, CONSTANT XREF: 0 0997 

2 1021 2 1026 2 1029 2 1033 2 1035 

SINGLE, ALIGNED, STATIC XREF: 0 0700 4 0737 NOT REFERENCED 

**** SEE STRUCTURE TEMPLATE ATMOSPROP 


SINGLE, 

ALIGNED, STATIC 

XREF: 

0 

0609 

4 

0617 

2 

0627 



OOUBLE , 

TEMPORARY 

XREF: 

0 

0762 

4 

0775 

2 

0776 

6 

0777 

2 0778 











SINGLE, 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 

0866 

2 

0868 

XREF: 

2 0582 2 0734 2 

0785 

0 

0981 







SINGLE, 

ALIGNED, INPUT-PARM 


XREF: 

0 

0901 

0 

0907 

2 

0956 

SINGLE, 

ALIGNED, ASSIGN-PARM 


XREF: 

0 

0851 

0 

0855 

4 

0876 

4 0898 

NOT REFERENCED 










SINGLE , 

ALIGNED, STATIC 

XREF: 

0 

0581 

4 

0583 

2 

0534 



SINGLE, 

ALIGNED, STATIC 

XREF: 

0 

0928 

4 

0956 

2 

0959 

2 

0966 

SINGLE , 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 

0929 

2 

0956 

SINGLE, 

ALIGNED, STATIC, 

CONSTANT 



XREF: 

0 

0930 

2 

0956 

SINGLE, 

ALIGNED, STATIC 

XREF: 

0 

0720 

4 

0747 

2 

0748 



ARRAY! 8), SINGLE, ALXGNEO 

, STATIC^ 

, CONSTANT 


XREF: 

0 

0931 

2 0959 

2 0966 










SINGLE, 

ALIGNED, INITIAL 

XREF: 

0 

0513 

2 

0636 

2 

0765 



SINGLE , 

ALIGNED, INITIAL 

XREF: 

0 

0513 

2 

0636 

2 

0765 



DOUBLE, 

ALIGNED, STATIC 

XREF: 

0 

0721 

4 

0730 

2 

0765 

4 

0782 

DOUBLE, 

ALIGNED, INITIAL 

XREF: 

0 

0093 

2 

0729 





DOUBLE , 

ALIGNEO, INITIAL 

XREF: 

0 

0094 

2 

0636 





DOUBLE, 

ALIGNED, STATIC 

XREF: 

0 

0722 

4 

0729 

2 

0730 

2 

0731 

2 0782 

2 0783 2 0788 2 

0798 

2 

0817 

2 

0818 

2 

0820 

4 

0821 

2 0824 

2 0826 4 0827 2 

0830 

2 

0832 

4 

0833 

2 

0836 

4 

0837 

SINGLE , 

ALIGNED, INPUT-PARM 


XREF: 

0 

0851 

0 

Q852 

2 

0872 

2 0876 











DOUBLE , 

ALIGNED, INITIAL 

XREF: 

0 

0095 

2 

0583 

2 

0586 

2 

0589 

2 0610 

2 0624 2 0625 










SINGLE , 

ALIGNED, STATIC 

XREF: 

0 

0723 

4 

0732 

2 

0747 

2 

0769 

2 0771 

4 0784 










DOUBLE, 

ALIGNED, INITIAL 

XREF: 

0 

0096 

2 

0586 





SINGLE, 

ALIGNED, STATIC 

XREF: 

0 

0724 

4 

0731 

2 

0732 

2 

0771 

4 0783 

2 0784 










SINGLE, 

ALIGNED, STATIC 

XREF: 

0 

0862 

4 

0873 

2 

0874 

2 

0876 

DOUBLE, 

ALIGNED, STATIC 

XREF: 

0 

0725 

4 

0798 

2 

0799 



SINGLE, 

ALIGNED, INITIAL 

XREF: 

0 

0480 

2 

0731 

2 

0783 

2 

0798 

SINGLE, 

ALIGNED, STATIC 

XREF: 

0 

0546 

4 

0559 

2 

0560 

4 

0561 

6 0562 

2 0563 2 0567 










SINGLE , 

ALIGNED, STATIC 

XREF: 

0 

0547 

4 

0563 

2 

0564 

4 

0565 
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DCL 

NAME 

TYPE 

ATTRIBUTES & CROSS REFERENCE 







548 

Z 

SCALAR 

6 0566 
SINGLE, 
2 0567 
DOUBLE , 

2 0567 

ALIGNED, STATIC 

XREF: 

0 

0548 

4 

0558 

2 

0559 

763 

513 

Z_PRED 

SCALAR 

TEMPORARY 

XREF: 

0 

0763 

4 

0776 

2 

0777 

ZETA_QDQT 

SCALAR 

SINGLE , 

ALIGNED, INITIAL 

XREF: 

0 

0513 

2 

0617 



2 OS63 



j 
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